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It is with great humility that I stand here 
today in order to attempt an evaluation of 
the contributions to the advance of meteorol- 
ogy and oceanography which my late teacher 
and friend, Vilhelm Bjerknes, made during his 
long and active life. He laid an ambitious plan 
for his work, he followed his plan with great 
consistency — and he lived to see many of 
his ideas bear fruit. 

To Bjerknes the prediction of the weather 
soon became the central problem towards 
which he directed his attention. In his en- 
trance lecture at the University of Leipzig, 
which in 1913 had called him to the chair of 
professor of geophysics and director of the 
new Geophysical Institute, he formulated his 
programme by saying: 

“From the observed state of the atmosphere 
at one moment we must compute the state 
at a later moment. The problem of the exact 
computation of a future state, which long ago 
has been solved in astronomy, must now be 
tackled in all seriousness in meteorology. 

“The problem is of immense dimensions. It 
can be solved only as result of a long devel- 
opment. One single research worker cannot 
get far, no matter how great his efforts are. 
However, I am convinced that in spite of 
all this it is not too early to point at this 
problem as the goal of the research. One does 
not always select a goal in the hope of reaching 
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it soon. To steer towards a distant goal deter- 
mines the course. In this case the far-away 
goal will give us an indispensable plan for 
our work.” 

Fifteen years carlier the name of Bjerknes 
had been unknown in meteorology. During 
his early years he had mainly been working 
in the field of classical hydrodynamics and 
had brought part of the remarkable work of 
his father, C. A. Bjerknes, to a conclusion. 
At this time he had frequently felt that the 
assumptions defining an ideal fluid placed too 
great limitation upon the character of the 
problems to which the results could be applied. 
In his search for an expansion he discovered in 
1898 his famous theorem of circulation that 
later on became the core of what he called 
physical hydrodynamics. Unknown to Bjerk- 
nes the same theorem had been discovered 
already in 1886 by Silberstein who had pub- 
lished his findings but had failed to draw any 
further conclusions. The very great merit of 
Bjerknes is that he immediately realised the 
importance of the theorem to the dynamics 
of the atmosphere and the ocean, and drew 
attention to the applications in his first paper 
on the subject. In this first paper he did not 
consider circulation relative to the earth, but 
this matter was taken up a few years later, in 
1902. 

From the very beginning Bjerknes’ papers 
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aroused great interest among meteorologists 
and oceanographers. One reason why Bjerk- 
nes’ work became rapidly known was that 
by a happy coincidence his ideas, which open- 
ed up new possibilities for dealing with the 
three-dimensional problems of the atmosphere, 
were advanced just at the time when the explo- 
ration of the free atmosphere started. Around 
year 1900 this fascinating field, for which 
Köppen proposed the name aerology, was 
suddenly being developed. In the 1890-5 
systematic meteorological observations from 
manned balloons had been organized in 
Germany by R. Assmann. From 1892 paper 
balloons with recording instruments ha 

been sent up by the French meteorologists, 
Hermite, Besangon and Teisserenc de Bort, 
and in 1894 A. Lawrence Rotch had started 
using kites at the Blue Hill observatory, Mass. 
In 1901 Assmann introduced the rubber 
balloon, which Teisserenc de Bort immediately 
took into use, and in the following year, 1902, 
he and Assmann announced independently 
the discovery of the stratosphere. 

In the midst of all this activity Bjerknes’ 
ideas of tackling the problem of the compu- 
tation of the future weather as a three-dimen- 
sional problem naturally aroused great hopes. 
The meteorologists who had turned with 
enthusiasm to the study of the free atmosphere 
foresaw the possibility that the observations 
gathered by their efforts might revolutionize 
the forecasting procedure. 

In oceanography the situation was simi- 
larly favourable. Around 1900 the accuracy of 
oceanographic observations was greatly im- 
proved and the technique of deep sea measure- 
ments perfected. Numerous observations were 
being collected, and now Bjerknes’ theoretical 
contribution promised to render the needed 
tools for the rational use of the observations. 

Still another fortunate circumstance con- 
tributed to the furthering of Bjerknes’ work. 
In this case I do not refer to a coincidence, but 
to a development for which Bjerknes himself 
has a great deal of the credit because he under- 
stood how to impart his own enthusiasm to 
others. In 1904 he visited the at that time re- 
cently established Carnegie Institution of 
Washington where he gave a number of 
lectures which created such an interest that 
from 1906 and until shortly before his death 
Bjerknes received an annual grant from the 
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Carnegie Institution in order to continue the 
development of his theories and apply them to 
the atmosphere and the ocean. 

I believe it will be difficult to point at any 
modest grant which has rendered greater 
returns than the one given Bjerknes. There are 
two reasons for this. One is that thanks to 
the grant Bjerknes could employ the assistants 
he needed in his work. He did need assistants, 
because his own work centered completely 
around the further development of the theoret- 
ical tools. He has to my knowledge never 
attempted to draw a weather map, nor has 
he ever discussed actual meteorological obser- 
vations. In the course of the years many of 
his assistants worked with the data, they tried 
to interpret the observed conditions, and, 
step by step, to gain better understanding 
of the physical processes in the atmosphere or 
the ocean. Bjerknes gave them complete 
freedom in their work. He was no hard 
taskmaster, but he laid the course. He was 
himself a very independent thinker, who in his 
writings rarely quoted results obtained by 
other authors. His papers are characterized by 
their elegant mathematics and clear style, 
because he never published a paper before he 
had corrected every ambiguous sentence and 
removed every superfluous word. He expected 
the same clarity from his assistants, and to 
some degree he wanted his assistants to be as 
independent in their research as he himself 
was. There were instances when he was 
openly displeased by finding his assistants 
spending too much time studying meteorolog- 
ical literature. In all events, Bjerknes’ assi- 
stants helped to advance his ideas, partly when 
they were working with him, and partly in 
later years, when many of them had students 
who worked under their direction. In this 
manner the Bjerknes influence has spread like 
rings in the water, reaching to ever greater 
distances. 

The other reason for the importance of the 
grant from the Carnegie Institution was that 
his grant was given during years when it was 
difficult for a young man to start on a scien- 
tific career. Bjerknes was free to offer assist- 
antships to any interested student and could 
select among the more promising ones. Several 
of the young men who for some periods 
became his “Carnegie assistants” did not 
continue in meteorology or oceanography, 
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but most of them did, and could later on, 
as already mentioned, bring Bjerknes’ view 
points to their own students. 

I have dealt to such an extent with Bjerknes’ 
assistants and his relation to them, partly 
because Bjerknes’ attitude to them was charac- 
teristic of him, and partly because speaking of 
Bjerknes’ influence on the development of 
meteorology and oceanography the contri- 
butions of his assistants cannot be overlooked, 
although Bjerknes was the dominating leader. 
He laid out the course, and his thinking ex- 
erted the greatest influence upon those who 
came in contact with him. 


After Bjerknes turned to the application of 
hydrodynamics to the problems of geophysics 
his activity falls in four different periods. The 
first was the period from 1898 to 1912 when 
he was in Stockholm until 1907 and in Oslo 
from 1907 until 1912. During these years the 
first two volumes of the broadly planned work 
‘Dynamic Meteorology and Hydrography”, 
the “Static” and the ‘‘Kinematic’’, were pub- 
lished by the Carnegie Institution, and during 


these years the impact of Bjerknes’ thoughts 
started to become apparent in the meteorolog- 
ical and oceanographic literature. Thus, in 
a paper of 1909 Bjerknes lists a dozen papers 
in which his theorem of circulation has been 
applied. 

He conducted an intensive campaign for 
the use of rational units in aerological research, 
and after he had urged his case at several 
meetings of the International Aerological Com- 
mission this body recommended in 1912 to the 
Directors’ Conference that in the aerological 
publications pressure should be reported in 
millibar and that geopotential should be used 
as the vertical co-ordinate. 

The second phase is covered by the years in 
Leipzig, from 1913 to 1917. Bjerknes accepted 
the ofter from Leipzig for two reasons. In the 
first place because he could expect to obtain 
more students and, therefore, might hope to 
develop a school in meteorology; in the second 
place because facilities were promised him 
for the working up of the aerological data 
that were being accumulated on agreed-upon 
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international days. Bjerknes would use these 
for a three-dimensional analysis of the state of 
the free atmosphere partly to develop methods 
for representing the state of the atmosphere 
by means of a series of charts, partly to gain 
experience as to the extent to which that partic- 
ular material was suited for a study of the 
changes in the state of the free atmosphere. 
His stay in Leipzig became much shorter 
than he had anticipated, and the results were 
less because only one and one half year after 
his appointment the first world war started, 
he lost his German students and also the needed 
assistance for preparing the aerological charts. 
As an interesting side-light on the role of 
meteorology during the first war, as con- 
trasted to its part during the last war, I may 
mention that no demands for participation 
in meteorological research related to military 
problems were made on Bjerknes or on his 
two or three Norwegian assistants, who con- 
tinued their strictly academic work and were 
paid by the grant from the Carnegie Institution 
in Washington even after the United States 
had entered the war. 

The Leipzig period was terminated in 1917 
when Bjerknes was asked to return to Norway 
and to accept the chair of professor of theoret- 
ical meteorology at the Geophysical Institute 
which recently had been established in Bergen. 
There he entered the most exciting period 
in his life, the period in which the importance 
of the fronts was recognized by his young 
collaborators, when the structure of the cy- 
clones was discovered and the term polar 
front was coined. The credit for the actual 
discovery of the structure of the cyclones 
and of the polar front goes to his collaborators, 
but Bjerknes himself had long been interested 
in the lines of convergence and divergence 
which appear on charts on which the wind 
direction is shown by means of stream lines. 
He had carly conceived the idea that a cyclone 
starts as a wave-like disturbance at a discon- 
tinuity surface, and one of his most important 
papers from the Leipzig period represents 
his first attempt at approaching that problem, 
the paper “Über Wellenbewegung in kom- 
pressiblen, schweren Flüssigkeiten”. 

In Bergen he had, however, to neglect his 
theoretical work because he and his assistants 
started immediately the organisation of a 
unique weather forecasting service, basing it 
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on inadequate observations, but using the 
newly gained insight in the importance of the 
fronts. This development exercised a consider- 
able influence on Bjerknes’ approach to the 
whole problem of forecasting, as evident from 
an address he gave at the meeting of Scandina- 
vian geophysicists in Gothenburg in August 
1918. He started by saying: 

“It is possible that different roads may lead 
to the solution of the problem of a satisfactory 
and practical method for weather forecasting. 
It is conceivable that among these there may be 
convenient roads which may lead around the 
strict requirement as to complete under- 
standing of the events whose development shall 
be predicted. Personally I am not interested 
in such roads. I am only interested in the road 
that leads straight through the inexorable 
demand of complete understanding.” 


But then he went on to say: 


“If we formulate the task (of predicting) 
as a mathematical problem, we may have no 
hope of finding a solution in a foreseeable 
future. But if we transfer it into a more 
practical form we find that we may advance 
simply by following consistently the road 
that the meteorologists started on when they 
first made the synoptic charts the basis of the 
study of the weather processes and by means 
of these found Bar’s law and so on. We do not 
need to doubt that along this road we may 
reach results which will give us increased 
understanding and thus will influence the 
meteorological practice.” 


In his Leipzig address in 1913 he had stressed 
the direct application of the laws of hydro- 
dynamics to given initial conditions. Now, 
after having been more in contact with the 
practical aspects of forecasting, he recognized 
the value of the step-wise approach, but he 
never relaxed his demand for complete under- 
standing of the physical processes. He held 
his banner high and under this banner meteorol- 
ogists and oceanographers from all parts of 
the world gathered. 


During the Leipzig and Bergen periods the 
applications to meteorology came into the 
foreground, and Bjerknes himself and his 
direct collaborators gave little thought to 
oceanography. Still, his influence on oceano- 
graphy has been equally great because for 
nearly fifty years Bjerknes demand as to 
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striving towards understanding of the physics 
of the events has directed the progress. 

The fourth and last period in Bjerknes’ 
life started in 1926, when he returned to the 
University of Oslo as professor of mechanics 
and applied mathematics. He retired in 1932, 
but remained active until shortly before his 
death. During these years he returned to his 
theoretical work. For a long time the third 
volume, ”Dynamics’, of his great work, 
“Dynamic Meteorology and Hydrography” 
had been in preparation. Drafts had been 
written to parts of the volume, but had been 
discarded because new ideas were advanced. 
Many of these appeared in the “Physikalische 
Hydrodynamik” which was published in 
1933 by himself and three of his collaborators, 
but the work with the “Dynamics” continued. 

I said that Bjerknes lived long enough to 
see many results of his life’s work. He saw 
aerology being developed to a degree which 
even in his most optimistic moments he had 
hardly considered possible, he saw the three- 
dimensional analysis of the atmosphere be 
made part of the routine work in forecasting, 
he experienced an expansion beyond expecta- 
tion of meteorological research and, best of 
all, he saw that most of this research was 
carried out along the lines he had advocated, 
aiming towards interpretation of events in 
terms of physics. Shortly before his death 
he had also the great satisfaction of seeing the 
manuscript of the “Dynamics” completed by 
himself and some of his early collaborators. 
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I cannot close a brief appreciation of Bjerk- 
nes’ activities without mentioning his very 
broad cultural interests. In Norway he was an 
untiring advocate for an expansion of the 
facilities of the university both in order to 
train better teachers and in order to create 
better conditions for scientific research. Along 
these lines he did not achieve as much as in 
his special field, but during the last years of 
his life he had at least the satisfaction that in 
Norway the possibilities for young men to 
use their abilities were increased greatly as 
compared to the years when he himself was 
neatly the only professor at the university 
who, thanks to the Carnegie grant, could 
offer a student an opportunity to attempt 
entering upon a scientific career. 


He was fortunate in retaining his clear 
mind to the last. Less than eight months ago 
he took active part in scientific meetings, and 
on his 89th birthday, only three weeks before 
his death, he was weak and had to stay in 
bed, but he insisted on seeing the many of the 
friends who called upon him to pay their 
respects. 


I am tempted to close with a poor transla- 
tion from the old Norse ‘‘Haavamaal” in 
which much of the wisdom of our forebears 
is condensed: 1 


‘“Cattle dies, man dies, 
you yourself shall die. 
Your name never dies.” 


Age Determination of Pacific Chalk Ooze 


by Radiocarbon and Titanium Content 


By G. ARRHENIUS!, G. KJELLBERG?, and W. F. LIBBY? 


(Manuscript received 1 november 1951) 


Abstract 


Measurements have been carried out on the average radiocarbon activity of a sediment 
core of eupelagic chalk ooze and the age of the lower boundary was calculated by integra- 


tion of the radiocarbon decay function. 


In such cores where constancy of the rate of minerogenous accumulation is indicated and 
can be checked by lateral control, the titanium content, calibrated in terms of absolute age 
by means of radiocarbon measurement, can be used for dating as far back in time as the 
constancy of the rate of minerogenous sedimentation can be proved. 


1. Principles of age measurement in core 61 B 
of the Swedish Deep Sea Expedition 


In the works of Jory (1908), PETTERSSON 
(1930), Piccor (1944) and Urry (references in 
Urry 1948) the possibilities have been de- 
monstrated for dating samples of pelagic 
clays on the basis of their content of ionium 
supported radium. The radiocarbon method 
(Lippy, ANDERSON, ARNOLD, ENGELKEMEIR, 
GROSSE, HAMILL, INGHRAM, KIRSHENBAUM, 
Reıp and WEINHOUSE 1946—51) has opened 
new possibilities especially with regard to 
the dating of highly calcareous deep sea 
deposits. 

If cores of deep sea sediments are to be used 
for geochronological purposes their facies 
relationships and stratigraphy should be well 
known. Furthermore such cores ought to be 
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correlated members of coherent profiles 
through the deep sea bottom, as this makes 
it possible to recognize layers deposited un- 
der abnormal conditions. Finally, if the 
half-life of the radioelement used is com- 
paratively short, an area with a compara- 
tively high rate of accumulation has to be 
chosen for sampling. 

The sediment cores raised by the Swedish 
Deep Sea Expedition under the leadership of 
Professor H. Pettersson provide an excellent 
material for studies of deep sea sedimentation. 
The closely spaced samples from the vast, 
stable Central Pacific area (cf. GUTENBERG and 
RICHTER 1949, p. 90) have been investigated 
and the results have made it possible to find 
a series of sequences which correspond in 
detail to the demands put forth above (ARR- 
HENIUS 1950 a and b). For the present 
purpose a member of this series was chosen, 
which had an exceptionally high rate of 
calcium carbonate accumulation, viz. the 
core 61 (10.3 m long, raised with the Kullen- 
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berg piston corer) together with 61 B (0.5 m 
ong, raised with a short open tube corer 
according to PHLEGER 1951), both from lat. 
0° 6° S, long. 135° 58’ W and 4437 m depth. 
In core 61 the uppermost layer is missing, 
owing to the sampling technique, but in 61 B 
the topmost layer is also present. This material 
was prepared for radiocarbon analysis in the 
spring of 1949. 

The cores 61—61 B have been closely cor- 
related with other cores in a profile extending 
from lat. 6° 44’ N to 3° 00’ S and it has 
been possible to prove that at the locality 61 
the changes in the rate of minerogenous 
sedimentation were smaller than the analytical 
error during the time of deposition of 
core 61B. This conclusion is based upon 
measurements of the total amount of biogenous 
material on one hand and of minerogenous 
material on the other hand, accumulated at 
different localities along the section and be- 
tween definite relative time limits. 

A closer description of the sedimentation in 
the East and Central Pacific will be published 
in the Reports of the Swedish Deep Sea Ex- 
pedition, Vol. V. 

As a measure of the minerogenous com- 
ponent the concentration of titanium recal- 
culated for salt free sediment has been used. 
Titanium is especially useful for this purpose 
as its concentration in the minerogenous 
component appears to be practically constant 
over the area in question and high enough to 
allow an easy and accurate determination. 
Furthermore titanium does not appear to be 
affected by secondary dissolution in this 
environment. 

The close proportionality between the con- 
centration of titanium and minerogenous 
matter in eupelagic deposits was first demon- 
strated in the fundamental works of CORRENS 
(1935) and Revere (1944) and was later 
stressed by Koczy (1950 and 1951) who used 
Correns’ measurements. With regard to the 
possibilities of estimating the rate of minero- 
genous sedimentation, which were conspicuous 
from these works, the analysis of titanium was 
included in the 1948 working program for 
treating the sedimentary material of the 
Swedish Deep Sea Expedition (G. ARRHENIUS, 
1950 a). 

In connection with his work on the core 
241 from the Equatorial Atlantic Ocean, 
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WISEMAN (1950) has also pointed out the pos- 
sibility of using the concentration of titanium 
for estimating the clay component. He ap- 
parently based an interpretation of Late 
Pleistocene chronology (1951) on the vertical 
distribution of minerogenous matter assuming 
a constant rate of deposition of such matter. 

It should be emphasized in this connection 
that a constant rate of sedimentation and a 
constant titanium content of the minerogenous 
component should not be looked upon as a 
general rule for deep sea sediments but must 
be proved in all cases. A section through the 
East Pacific from lat. 0° to 17° N and at an 
approximate distance of 800 nautical miles 
from the American coast did not reveal the 
high regularity of minerogenous sedimenta- 
tion found in the Central Equatorial area. 
Frequent examples of irregularities in the 
rate of minerogenous deposition are found 
on rises and ridges and in throughs and pits 
of the Atlantic, Pacific and Indian Oceans. 

It should also be pointed out that a smooth 
bottom topography or freedom from vol- 
canism, as in parts of the stable Central 
Pacific area or in parts of the Indian and 
Atlantic oceans, might considerably reduce the 
probability for tectonic disturbances, sediment 
slumping and abnormal dilution with pyro- 
clastic material but is by no means a guarantee 
for a constant rate of minerogenous sedimen- 
tation. Redepositional changes in this rate by 
as much as 50 % appear in the eupelagic envi- 
ronment generally not to be accompanied by 
changes in grain size or by qualitative changes 
in the composition of the clay material. The 
only possibility of checking changes in the 
rate of minerogenous sedimentation in a 
pelagic environment therefore seems to be 
a quantitative comparison of corresponding 
layers in several sediment cores along a profile 
of wide geographical range. 

The biogenous components of the sediment 
consist mainly of skeletal remains of calcium 
carbonate from foraminifera and coccoli- 
thophoridae and of silica from radiolaria and 
diatoms. The aim of the present investigation 
was to use the radiocarbon of the calcium 
carbonate for dating purposes. The possibility, 
however, that the calcium carbonate had 
exchanged ions with the bottom water after 
deposition could not be ignored. In such a 
case it would not be possible to use the radio- 
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carbon content of the carbonate for age 
calculations. An important means for control- 
ling secondary exchange is furnished by the 
palaeothermometric method of UREY based 
upon the temperature sensibility of the ratio 
between the oxygen isotopes O16 and O3 
(Urey, LOWENSTAM, Epstein and McKInney, 
1951). If the calcite had exchanged all its 
carbonate with carbonate ions from the deep 
sea bottom water or from the interstitial 
sediment water, a temperature of a few degrees 
centigrade ought to be indicated by the oxygen 
isotope distribution of the calcite. If, on the 
other hand the shells have not been subjected 
to any exchange of carbonate after the death 
of the organisms, the temperature indicated 
should be that of the surface layer of the ocean. 
If partial replacement of the original calcite 
had taken place some intermediate temperature 
would be found. 

Measurements of the oxygenisotope distribu- 
tion carried out by Dr S. Epstein gave the 
following results in core 61: 0—2 cm 21°C, 
TOI 57cm at, 20 Cm 10.5 48200 em 173°; 
all measurements with an estimated standard 
deviation of 5°C. All temperatures found 
thus coincide within the limits of the experi- 
mental error with an assumed average tempera- 
ture for the biotope of tropical planktonic 
foraminifera which are most abundant in 
the uppermost 100 m layer of the ocean 
(PHLEGER 1951). It is thus very probable that 
no considerable exchange of the shell carbonate 
took place after deposition on the deep sea 
bottom and the radiocarbon measurements 
might safely be used as a basis for age calcula- 
tions. This conclusion refers to core 61 only 
and the results may not directly be applied to 
other types of deep sea deposits or other 
localities. It is, however, probable that the 
stability of the carbonate characterizes all 
bathyal pelagic deposits with a very high 
content of calcite, poor in humus and not 
affected by rising juvenile solutions.* 

The stratigraphy of core 61 and 61 B in- 
dicated 3 cm/1,000 years as the order of magni- 
tude of the average rate of Post-Glacial sedi- 


1 It might be pointed out in this connection that 
radiocarbon measurements in autochtonous eupelagic 
carbonate deposits are not affected by the source of 
error recently indicated by BARTIETT (1951) as in- 
fluencing the isotopic composition of tufa and lake 
marl, 
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Fig. 1. Principles of age calculation. T is the age of 
the lower boundary of an ideal sediment column, which, 
subjected as a whole to radiocarbon analysis, yields the 
activity N,,. The curve represents the radiocarbon decay 


function. As is easily seen, the surface of the rectangle 
N: T is equal to the surface limited by the coordi- 


nate axes, the line y— T and the decay curve. 


mentation at this locality. The half life level 
of C!4 would under such conditions be found 
at approx. 17 cm depth. For the radiocarbon 
analysis approx. 100 g of calcium carbonate 
are required. Due to the small dimensions of 
the core (3 cm diameter) and to the fact 
that material had to be used also for other 
analyses it was, in this case, impossible to 
extract and analyse samples from definite time 
levels and to construct a decay curve from the 
results. Instead a coherent sample from o to 
45 cm depth in core 61 B was subjected as a 
whole to radiocarbon analysis and the result 
was used as a basis for the age calculation. The 
principles of this calculation appear from fig. 1. 
If N,, is the average activity value obtained 
from analysis of a homogenous sediment 
column, deposited with a constant rate of 
accumulation and if T is the age of the lower 
boundary of the column then the product 


a. 
Nu T = f No: e—* + dt if N, = No: e—# 
O0 


is the disintegration function. 

In the actual case, however, the rate of 
sedimentation has not been constant, the 
changes being caused by changes in the rate 
of biogenous sedimentation. As the miner- 
ogenous rate of sedimentation appears to have 


AGE DETERMINATION OF CHALK OOZE 


been constant in the case of core 61 B and its 
neighbours and as the concentration of titanium 
is proportional to the concentration of miner- 
ogenous matter (see above pp 223—224 and 
Vol. V, Rep. Swed. Deep Sea Exp.), the in- 
tegral over the titanium concentration as a 
function of the depth in the sediment will be 
proportional to the time interval between the 
integration limits. The relative age of any 
level in core 61 B is thus expressed by the 
titanium integral from the present sediment 
surface to the level in question. 

A factor which must be taken into con- 
sideration in the age calculation is, that due 
to the changes in the balance between the 
biogenous production of calcium carbonate 
and of silica the amount of calcium carbonate 
and therefore of radiocarbon deposited per 
time unit varies with time. As the calcium 
carbonate concentration has been measured in 
closely spaced samples from core 61 B, cor- 
rections for the varying carbonate deposition 
can easily be applied. 

The calculation must also take into account 
the fact that at different levels of the sediment 
column samples of varying size had been ex- 
tracted for other chemical and palaeontolog- 
ical analyses previous to the radicarbon anal- 
ysis. The corrections for this stepwise variation 
in size of the samples were applied in the 
same way as the corrections for carbonate 
content. 

The variations in composition of the samples 
investigated were not big enough to require a 
correction for the varying specific gravity 


of the deposit. 


2. Constants and experimental values 


Assay of modern calcium carbonate shell 
No = 16.2 + 0.5 disintegrations/min/g of C. 

Half life of C14 = 5 568 + 30 years. 

Assay of column 61B N, = 6.18 + 0.24 
disintegr/min/g of C. 

The results of the measurements of carbonate, 
titanium and mass of sediment are shown in 
table 1 together with the calculated cumula- 
tive fractions of the total time T. The core 
was divided into 22 elementary samples, or 
strips, for each of which the same measure- 
ments were made. The concentrations of car- 
bonate and titanium and the mass values all 
refer to salt free sediment. 
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TT Cumul. 
Sample | [CaCO;] | [TiO:] le fraction 
depth +108 =O? SP | of total 
in dg ; 
cm (qi) On time T 
; (fi) 
o—2 762 33 72) 0.054 
Te 764 31 55 0.104 
Ay 766 30 13 0.141 
ST, 780 24 42 0.170 
DAS 795 18 RO 0.192 
8.$—11.5 795 18 83 0.236 
11.5713.5 794 19 24 0.267 
1955 15.5 805 18 55 0.296 
LS? 817 17 21 0.316 
I Fas 840 22 I4 0.325 
17.5—19 862 27 27 0.358 
19 —2I 808 44 55 0.430 
2I —22.5 755 61 14 0.504 
22.5—23 VS Di 13 0.527 
2 25% 760 38 138 0.604 
28.5 2725 763 22) 79 0.640 
27.5—32 768 22 248 0.738 
3:2) 05S 773 23 71 0.766 
83.35-38.45 778 26 275 0.872 
38.5—40 783 29 73 0.908 
40 —43 799 25 166 0.969 
43 —45 815 20 $5 I.000 


3. Calculation of the results 


Notations: 
mass of material in the i:th strip p; 
fraction of carbonate in ,, a qi 
time limits of , Re 
; 3? Pag i 


T is the age of the lowest level of the column. 
The index i varies from 1 through 22, and 
fo =O; ja2 = 

To derive an equation that will permit 
the determination of T, it is simplest to 
start with the assumption that the quantities 
needed are known as continuous functions of 
the time ¢ along the column. 

Suppose that c(t)dt is the mass of the car- 
bonate contained between the levels of the 
core that correspond to times t and t + dt. 
With No = the radioactivity of modern car- 
bonate per unit mass, the total activity of the 
core becomes: 


72 
Na No S c(t)e—* dt (1) 
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and the activity per mass unit carbonate: 
m 
S c(t) e—* de 
Nm = No : — (2) 

J c(t) dt 

oO 
But c(t)dt = q(t) - m(t)dt, where m(t)dt is the 
mass of all material between the levels t and 
t + dt, and q(t) the fraction of carbonate in 
the same portion. The variation of m(t) thus 
accounts both for the variations in the rate of 
sedimentation and for the variations in the 
area of the cross-section which occur because 
material was extracted from the core for 
other experimental purposes, prior to this in- 
vestigation. q(t) depends also on the variation 
of the rate of accumulation of silica. 

Inserting this in (2) and removing the 
denominators we get: 


Nm S ad ms) de = No f q(t)m(t)e—*dt (3) 


[0] 


or, by writing the integrals as sums of in- 
tegrals over the elementary strips: 


j=22 4 


= No 22 Ji q(t) m(t)e—** dt (4) 


Supposing now that q(t) and m(t) are con- 
stant over each strip, we get: 


1— 22 


Nin >> qimi(ti — t;—1) = 


ı=I 


i—=22 
No’ 


n z qim; [exp (— At;_,) — exp (— At;)] 
(5) 


As m;(t;—t—1) = pi, and.t; = fi - T, we 


finally get: 


=) ee The 
— exp (— ATfi)| (6) 
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Putting: 
P22, N, 
= na | Xt 
i=I 


_ Ah) = exp(— | 
aS | 


the upper age limit T is the positive root of the 
equation u(t) = o. 

Although this equation might conceivably 
have several positive roots, it can be shown 
that this is not the case by study of the exact 
equation (3), which can have only one positive 
root. This conclusion is valid also for discon- 
tinuous q(t), m(t), such as we have assumed 
here, and the equation u(t) = o does therefore 
determine T unambiguously. 

The value of T was found to be 14 162 
years. The accuracy of this result is discussed 
in the next section. 

For the numerical work, the digital relay 
computer BARK was used.t A program for the 
evaluation of # was so arranged that, with 
only slight modifications, it could also be used 
to calculate the different partial derivatives 
necessary for the error estimation. (See next 
section). Once this program had been set up, 
each new value of # was computed in about 
5/2 minutes. It was therefore easy to com- 
pute both u and the derivatives for several 
neighbouring values of t, which gave a con- 
vincing guarantee against accidental errors of 
computation. 


4. Sources of error and their influence 


If, for the purposes of the present section, 
we introduce the notation: 

dx = the standard deviation of x, we get 
the following formula for dT. 


ine (Ge) de + (Gx) ane + 
OT [7 


ST 3 | i= 22 IT D 
+ (>) AN, Ir D: (&) dp;? Ar 


+ /0T\: |. IT\: 
+ YG) @ ty G) 4 


* For a description of the BARK and its performance, 
see G. KJELLBERG and G. Ngovius, (1951). 
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i= 22 Dy \2 t= 22 
PACE 


IT 


The derivatives were computed according 
to the formulae: 


Draw 7 

Pr) 6CO«CS 

ou Æ — _ NT 
ON z BS No? 


| FROM MIE RD oe 
fe — 
Nm 


oe ae aT 
2 N, ART 
PAT 0 ep) 


i fr 


ou u ou 
oa Gi dpi 
ou Pidi 


9 — fi—fi-s 
(ea) opti) 
hi 
WIEN: 
— AT exp(—aTf | er; 
{ exp(— ATf) — exp(— ATfi+1) 
tr 


— AT exp(— ATfi) | 
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(For the computation of the derivatives 


Où du Ou du 
Ot? Op; dar fi 


the Bark was used.) 


The following numerical results were ob- 
tained. The units employed are: + and T in 
years, N,, and No in disintegrations per min 
per g C, p; in hectograms, q; and f; dimen- 
sionless fractions.! 


oT h 
Or 2x 54 
ce = 0.48712 - 10-4 
Die bos ae, : 
2 = 0.050 
ON en 
UE A 
ON Zu 333 
and 
du du du 
Opi 9; dfi 
I — 0.7689 — 0.0121 + 0.0716 
> — 0.6577 — 0.0474 + 0.0675 
3 — 0.5743 — 0.0097 + 0.0526 
4 — 0.5204 — 0.0280 + 0.0660 
5 — 0.4844 — 0.0195 + 0.0996 
6 — 0.4270 — 0.0446 + 0.0900 
9) — 0.3637 — 0.0110 + 0.0613 
8 — 0.3223 — 0.0220 + 0.0554 
9 — 0.2913 — 0.0075 + 0.0242 
10 — 0.2785 — 0.0046 + 0.0330 
II — 0.2530 — 0.0079 + 0.0552 
12 — 0.1680 — 0.0II4 + 0.0410 
13 — 0.0766 — 0.0014 + 0.0128 
IA — 0.0294 — 0.0005 + 0.0683 
IS + 0.0165 + 0.0030 + 0.0896 
16 + 0.0635 + 0.0066 + 0.1219 
197) + 0.1141 + 0.0368 + 0.1078 
18 + 0.1576 + 0.0145 + 0.1045 
19 + 0.1992 + 0.0704. + 0.0942 
20 + 0.2389 + 0.0223 + 0.0590 
DT + 0.2680 + 0.0557 + 0.0504 
22 + 0.2946 + 0.0199 


1 In section 3 No and Ny were defined as disin- 
tegrations per min. per g CaCO,. It is, however, 
obvious that the choise of unit does not matter, as 
long as No and Ny are measured with the same unit. 
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The value of the standard deviation of each 
variable is known from the measurements, 
except for p; and qj, where one should also 
take into account the error involved in the 
assumption that q(t) and m(t) are constant over 
each strip. This leads to estimated s.d.:s for 
the q;, larger than would result from_ the 
measurements alone, while the errors of the 
pi are not affected. The argument proving this, 
and the procedure by which the standard devia- 
tions of the g; were estimated (with the aid 
of differences from table 1), will be omitted, 
as it will soon appear that neither group of 
errors contributes appreciably to the total error 


dT. 


After these corrections had been made, the 
standard deviations of the p;, q; and f; were 
taken to be, respectively: 


1° 107-3, between 3- 1073 and 9: 1073, I: 1073 
The standard deviations of the other variables 
are given in section 2. Inserting the values of 


derivatives and s.d.:s in equation (7), one 
obtains: 


dT = 843 
and we should therefore write: 
T = 14 200 + 900 years. 


The following list gives the influence of 
each separate source of error. 


oT 

le # =” 
OG 

| IN, dNo — 522 
PE 


EF 
SLY 
en 

+ 
= 
we 

A 
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It is evident that the errors dN,, and dNo 
dominate completely, and in particular that the 
errors in p, q and f are negligible. 


5. Possibilities for extrapolation of the 
results and further use of the integration 
method for age determination 


The possibilities of extrapolating laterally 
the results gained from core 61 B depend upon 
whether isochronous key levels can be found 
and traced through the sedimentary sequence 
of the oceans. The investigations of the Swed- 
ish Deep Sea Expedition has shown that this 
condition is fulfilled in certain types of calcium 
carbonate lithofacies where detailed correlation 
can be carried out over distances exceeding 
1 700 nautic miles. The world encircling net 
of stations of the Swedish Deep Sea Expe- 
dition provides a material which is highly 
promising for the establishment of a global 
marine geochronology (cf. PETTERSSON, 1950). 

The possibilities for extrapolating backwards 
in time from the radiocarbon age determina- 
tion of core 61 B depend upon the ability to 
evaluate the rate and changes of rate of 
biogenous and minerogenous sedimentation 
in the same way as was done with core 61 B. 
This work has not yet been completed but 
as far as can be judged from present results, the 
vertical extrapolation can be carried down to 
the lower Pleistocene or upper Pliocene in the 
East Pacific chalk ooze profile. This mostly 
corresponds to the entire depth of the se- 
quence sampled with the Kullenberg piston 
corer. 

The easiest way of extending the radiocarbon 
time scale backwards in time appears to be 
to control the constancy of minerogenous 
sedimentation and, in such cases where con- 
stancy is found, to calibrate the titanium in- 
tegral in terms of absolute time. The age of 
any depth level can thereafter be calculated 
from the corresponding titanium integral. 
Constancy of minerogenous sedimentation can 
be expected and proved only in exceptional 
cases. In regions where no such constancy is 
found the dating can be acquired through 
correlation with key localities such as number 
61 of the Swedish Deep Sea Expedition. 
Correlation methods for deep sea sequences 
will be described in detail in Vol. V of the 
Reports of the Expedition. 
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Multiple Currents in the Gulf Stream System” 


By F. C. FUGLISTER, Woods Hole Oceanographic Institution 


(Manuscript received 20 September 1951) 


Abstract 


A new interpretation of the accumulated temperature and salinity data from the Gulf 
Stream Area indicates that the System is made up of a series of overlapping currents. 
These currents are separated by relatively weak countercurrents. Data from a recent survey 


are presented as supporting this hypothesis. 


In nearly all discussions of the Gulf Stream 
and of the North Atlantic Current the terms 
countercurrent and eddy are frequently em- 
ployed, (ISELIN 1936, 1940, and SVERDRUP 
et al. 1942). Many temperature and salinity 
profiles across the currents indicate the pres- 
ence of countercurrents and/or eddies and 
certain theoretical considerations require that 
they exist (RossBy 1936). On the other hand 
charts showing the average distribution of 
temperature and salinity (Wüst and DEFANT 
1936 and FUGLISTER 1947, etc.) and charts 
showing the currents in the western North 
Atlantic do not indicate that countercurrents 
or eddies are permanent features of the sys- 
term. 

It is the purpose of this paper to suggest that, 
by eliminating the averaging processes, coun- 
tercurrents can be introduced as, at least, semi- 
permanent features of the circulation and a 
new picture of the current system will evolve. 
The existing current charts fall into two 
general classes, those that show a relatively 
narrow current flowing from Cape Hatteras 


* Contribution No. 584 from the Woods Hole Oce- 
anographic Institution. 
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to the Grand Banks area where the current is 
shown to fork or branch, and those charts 
that depict a current at Cape Hatteras that 
gradually broadens as it moves toward the 
mid-Atlantic. The data used for these charts 
interpreted in a different manner shows that, 
in the Gulf Stream and North Atlantic Current 
areas, there is no one continuous current but a 
series of overlapping easterly currents separated 
by countercurrents. 

The averaging processes tend to obscure 
certain details that are known to exist. Per- 
haps the most important of these details, in 
the present connection, is the very narrow 
zone of minimum temperatures or low pressure 
generally found on the left-hand side, facing 
downstream, of a strong current. In an attempt 
to show this zone and to illustrate what seems 
a more logical interpretation of the accumulat- 
ed sub-surface data, a schematic chart of 
temperatures at a depth of 200 meters in the 
Gulf Stream area has been drawn. The data 
used are from various sources but primarily 
from the multiple ship survey of 1950, (Fuc- 
LISTER and WORTHINGTON 1951). This chart 
(see fig. 1) shows three easterly currents 
between Cape Hatteras and the longitude of 
Flemish Cap. There is a minimum temperature 
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zone inshore of the sharp gradient off Cape 
Hatteras. The temperatures in this zone in- 
crease to the east. This increase in temperature 
is balanced by a corresponding decrease in 
temperature in the warm zone to the right of 
the gradient. If we assume for the moment 
that the currents flow along the isotherms 
then the major current off Cape Hatteras, 
without becoming wider or branching, simply 
becomes weaker as it flows toward the east. On 
both sides of this current there are weak 
countercurrents. A second easterly current 
forms just off the continental shelf north of 
Cape Hatteras, increases in strength as it 
flows toward the Grand Banks area and then 
in its turn gradually weakens in the same 
manner as the first current. Finally a third 
stream is depicted forming in the Grand Banks 
area and moving toward the northeast. 

In order to check on the above interpretation 
of data obtained from various sources a survey 
was made of part of the area by the AIBA- 
TROSS III in June 1951. The data collected on 
this cruise consisted of half-hourly bathy- 
thermograph observations to a depth of 250 


Schematic chart of temperatures (°C.) at a depth of 200 meters in the Gulf Stream Area. 


meters and (over many sections) hourly sur- 
face current measurements with the GEK (von 
ARX 1950). As was the case with the multiple 
ship survey of 1950 the majority of surface 
current vectors showed the same current 
directions as would be deduced from a plot of 
the 200-meter isotherms. A plot of these 
isotherms and the track of the ALBATROSS III 
are shown in figure 2. The entire area of 
interest could not be covered in the 26 days 
the research vessel was at sea but the results 
obtained strongly support the multiple current 
hypothesis. In fact the data indicated the 
existence of four major easterly currents separat- 
ed by countercurrents of varying magni- 
tudes. 

Aside from the extravagant amount of 
extrapolation used in both fig. ı and fig. 2 
the most serious criticism of the multiple 
current hypothesis as presented here may be 
that the temperature structure at 200 meters 
does not reflect the pressure field. It may be 
considered as representing only a relatively 
shallow surface layer. Actually a similar pic- 
ture of the current regime in this area could 
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Fig. 2. Temperatures (°F.) at a depth of 200 meters in the Gulf Stream Area. An interpretation of the data 
obtained on the ALBATROSS III cruise of June 1951. 


be obtained by using surface dynamic height Of the available sections to the east of th 
anomalies referred to the 2000 decibar sur- Montauk—Bermuda line and west of the soth 
face. The difficulty would be that observations meridian only 7 reach from the continenta 
taken over a period of more than ro years shelf out to the Gulf Stream. These 7 all sho 
would have to be used and even more extra- an easterly current off the continental shel 
polation would be required. and a countercurrent between it and the Gul 
However, an examination of all the available Stream. None of these sections reached fa 
deep temperature and salinity sections normal enough south to check on the existence of 
to the coast between Cape Hatteras and the countercurrent to the right of the Stream 
Grand Banks gives the following interesting When the data from the three most easterl 
results. Out of fifteen sections between Mon- sections are averaged, the depth of the sigma- 
tauk, Long Island, and Bermuda, variations 27.2 surface at each maximum and minimu 
in the depth of the sigma-t 27.2 surface show: point from the shelf outward is: 
7 sections with a slight westerly current 
very. close in by the continental shelf: metre 7732 = 387 383 89 
13 sections with a stronger easterly current difference —58 HAS = COMTESSE 
just off the shelf; 
13 sections with a countercurrent to the left 
of the Gulf Stream. 
12 sections with a countercurrent to the 


Assuming that the differences in depth o 
the sigma-t surface between these points in. 
dicate the relative strength of the current 


: then it appears that the easterly current o 
right of the Gulf Stream. : 
8 ult Stream the continental shelf increases and the Gul 


The average depth in meters of the 27.2 SI8~ Stream decreases in strength as they mov 
ma-t surface at each maximum and minimum toward the Grand Banks area 


point from the shelf outward is as follows: If this multiple current system is accepte 


average 255 231 336 214 949 833 then there is a problem of nomenclature. As 
difference —24 +105 —122 +735 —116 group the currents discussed here form 
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portion of the Gulf Stream System. This 
portion of the System fits the description of 
the North Atlantic Current, as given in THE 
OCEANS (Sverprup et al. 1942), except 
that now the name should apply to the group 
of currents as far west as Cape Hatteras in- 
stead of only to those currents to the east of 
the Grand Banks. If the term Gulf Stream is 
to be retained it should be understood that 
it applies only to the direct continuation of 
the Florida Current and it does not include 
all easterly-flowing currents between Cape 
Hatteras and the Grand Banks. Whether or 
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not any of the other currents in the group 
should be specifically named is a question that 
might wait until more is known about their 
permanency. 

Since June of this year the Arsatross III 
has been engaged in surveying the area to the 
east of Flemish Cap. Aboard ship the various 
currents are referred to as No. 1 (for the Gulf 
Stream), No. 2, No. 3, etc., for each successive 
current encountered. No. 4, the current off 
Flemish Cap, has been followed as far east 
as the 23rd meridian and the possible exi- 
stence of a No. 5 has been indicated. 
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Height Variations in the Concentration of Ions Near 
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By HARALD NORINDER and REINHARDS SIKSNA 
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Abstract 


Results of measurements of the concentration of small- and large-ions at different heights 
near the ground during quiet summer nights at Uppsala are given. Characteristic shapes 
for the variation of the concentration of ions are illustrated. At the ground the concentra- 
tion of small ions increases more than at higher levels, where also an enhanced content of 
small ions can be observed during the stated weather conditions. The number of large ions 
descreases after sunset reaching a minimum at sunrise, after sunrise a comparatively rapid 
increase follows. It seems that the concentration of large ions near the ground is actually 
a little below that at the higher levels. Preliminary explanations of the variations are 


outlined. 


Introduction 


We have developed the concept that some 
of the variations of the ionic density observed 
by us in the summer of 1949 (17) were caused 
by the accumulation near the ground of emana- 
tion exhaled from the soil, and consequently 
also by its decay products, under certain 
meteorological conditions in which a minimal 
mixing of air takes place (18), (19). A detailed 
investigation of the variations of the ionic 
density, under meteorological conditions which 
might favour the accumulation of emanation 
near the ground, was carried out at the In- 
stitute for High Tension Research at Uppsala 
in the summer of 1950; the ionic density was 
measured at different levels above the ground 
in order to obtain an insight into the mixing 
conditions prevailing. Measurements of ele- 
ments characteristic of the electrical properties 
of atmospheric air at different levels are not 


very numerous, though they may be of great 
importance. Investigations of the potential- 
gradient at three different levels have been 
carried out by DAUNDERER (6) and by Norin- 
DER (15). BROWN (4) has measured the space- 
charge in atmospheric air at different levels 
near the ground. BECKER (1) and Dey (7) have 
measured the content of emanation at two 
levels. Conductivity of the air, mainly due 
to small ions, in the lowest levels of the 
atmosphere has been measured at different 
heights by HoGG (10). The vertical distribu- 
tion of atmospheric condensation nuclei has 
been measured by YUNKER (24) at two levels. 

We have further investigated the electrical 
elements of atmospheric air near the ground 
at different levels by measurements of the 
ionic density. In this manner we hope to gain an 
insight into the matter of ionization processes 
near the ground and their correlation with 
the conditions of mixing of air there. 
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Arrangement of the measurements 


To be able to measure the ionic density at 
different levels above the ground a special 
device was built by which air was sucked from 
three levels above the ground (h, = o.1, 
0.05 m, ha = 0.5, 0.25 m, h, = 2.1 m). The 
ionic density of the air sucked from these 
different levels was measured with the ion- 
counters used during the preceding measure- 
ments at the Institute (16), (17), (23). Small 
ions of both polarities were measured with 
two reconstructed EBERT ion-counters (16) 
and the positive large ions with an ISRAËL 
large-ion counter (23). Some difficulties arose 
during the measurements. At sunset insects 
were sucked into the EBERT counters and 
remained on the insulators of the electrometers. 
The apparatus could be cleaned by compress- 
ed air or by dissembly. Another kind of 
difficulty was encountered with the large-ion 
counter during this series of measurements; the 
insulation losses of the counter were variable, 
increasing abnormally during the night and 
diminishing again during the early morning 
hours. It seemed at first that measurements 
of the large ions were impossible because of 
the above mentioned abnormal irregularity 
in the deterioration of the insulation. It may 
be supposed that the most probable cause for 
the observed diminution of the insulation was 
a radioactive deposit, precipitated on the 
condensers of the counter during the aspiration 
of air through it. If the number of large ions 
must be measured with a large-ion counter 
in the manner used by us during selected 
summer nights when an enhanced quantity 
of active deposit may be expected, especial 
attention must be paid to the determination 
of the fictive insulation “losses” of the ap- 
paratus as a function of the time of aspiration 
and of the time lapse after the interruption of 
aspiration. 


Results of measurements 


The measurements of the ionic density were 
carried out during selected nights when 
conditions could be expected to be favourable 
for the accumulation of emanation because 
of stratification of the air near the ground. 
The nights selected were without clouds and 
windless. Variations of the measured concen- 
tration of small and large ions during some of 
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Fig. 1. Variation of concentration of ions at two levels 
on July 16—17, 1950; positive small ions, - ---- 
negative small ions, positive large ions. 


the nights selected are shown in figs. 1, 2, 3, 4, 
and 5. The positive ionic density is denoted 
by the full line ) and the negative 
ionic density by the dashed line (----- 
Values corresponding to different levels are 
denoted by different designations as shown 
in the figures. The noted time is Central- 
European. The temperature given in the 
figures may be regarded as a factor, character- 
istic of the meteorological conditions pre- 
vailing during the observation periods. 

Small ions. From an examination of the 
variations shown in figs. 1, 2, and 4 two 
superposed kinds of variation of the concentra- 
tion of small ions can be identified: 

(1) Variations of a longer time scale which 
extend over the entire observation period, 
and (2) variations of a shorter time scale with 
a duration of one to several hours. Variations 
of still shorter period may exist, but are not 
detected with the arrangement of measure- 
ments used. For a general insight into the 
basic effect variations of the first kind are of 
fundamental interest. The variations of a 
shorter time scale can be regarded as the 
result of variations in the environmental con- 
ditions. The increase of the ionic density at 
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Fig. 2. Variation of concentration of small ions at three 
levels on July 20—21, 1950; positive small ions, 


Fig. 4. Variation of concentration of small ions at three 
levels on August 23—24, 1950; positive small 
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the two higher levels (h, = 0.25 m and 
hs = 2.10 m) at 22/ in fig. 4 can be regarded 
as an extreme example of such short-period 
variations. This increase was perhaps caused 
by air, with a higher content of ions or ion 
producers, rising locally from the ground. A 
stronger decrease can be observed at all three 
levels following the increase of ionic density 
at the higher levels h, and hg. A similar effect 
can be observed during the same night be- 
tween 2" 30" and 3" 40". In this way the 
shorter-period variations may be regarded as 
local disturbances which deform the general 


made to smooth the curves of the measured 
trend of the small ions detected, neglecting 
variations of a shorter time scale. These 
smoothed curves which may be regarded 
as idealized trends of the effect are shown in 
fig. 6. 

The shape of the variations of the small ion 
concentration. The following features of the 
variation in concentration of small ions during 
the effect may be pointed out: 

(1) At the ground an increase is observed 
which proceeds very rapidly at sunset; during 
the night the ionic density remains compara- 
tively high with only a slow decrease until 1—2 
hours after sunrise when it decreases rapidly. 
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(2) At the higher levels an increase also 
begins at sunset and continues until 2—3 hours 
after sunrise; after that time a rapid decrease 
follows. 

(3) At higher levels the number of small 
ions is smaller. 

(4) The concentration of positive ions is 
higher than that of negative ions. 

(s) The sudden increase in the concentration 
of negative ions over the positive ions at the 
beginning of the increase at sunset (fig. 1) 
deserves attention. A similar phenomenon was 
observed also during other series of measure- 
ments, especially at higher levels (for example 
see fig. 2). 

(6) In addition to variations of a relatively 
long time-scale (extending over the entire 
observation period), variations of a shorter 
time scale (with a duration of one to several 
hours) may be observed. 

Large ions. As mentioned before, it seemed 
impossible at first to obtain measurements of 
the large ions because of the abnormal de- 
terioration of the insulation of the apparatus 
during the measurements. For the same reason 
we were unable to determine the spectral 
composition of large ions with the method 
used earlier (23). The measurements of large 
ions during these series were carried out at 
an adjusted limit mobility ky = 2.14 X 10-1 
cm sec-l: V cm-!. Taking into account the 
trend of the insulation it was attempted to 
calculate the variations of the concentration 
of large ions during the effect. Surprisingly, 
the trends obtained for the variation of the 
concentration of large ions were of a charac- 
teristic type. Two families of smoothed curves 
of these variations are shown in fig. I using 
different combinations of the observed curves. 
The general features of the variations of 
large ions during the effect are as follows: 

(t) A decrease may be noted from sunset 
to sunrise, when a minimum is reached; after 
sunrise a comparatively rapid increase follows. 

(2) It appears that the concentration of 
large ions near the ground is really a little 
lower than that at the higher levels. 

The trend of the concentration of large 
ions is given in fig. 7 unexpectedly exactly; 
only the height of the maximum after sun- 
rise could not be determined exactly with the 
data at hand. It must be noted here that by 
evaluating the trend of the variations of large 
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Fig. 6. Idealized trend of the variations of the concentra- 
ion of small ions at the ground (h,), and at two higher 
levels (hy < hg) during quiet summer nights at Upp- 
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Fig. 7. Idealized trend of the variations of the concen- 

tration of positive large ions at the ground (h,), and 

at two higher levels (hg < h,) during quiet summer 
nights at Uppsala. 


ions during individual nights and also by the 
smoothing of the curves shown in fig. 7 every 
effort was made to suppress the variations; 
nevertheless a characteristic trend remained. 


Points to be kept in mind in explaining 
the variations 


We have previously given a preliminary 
general explanation of the observed variations 
of the concentration of ions during certain 
specified weather condition (18), (19). Indica- 
tions will be given here, of a way by which a 
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more detailed explanation of the whole effect 
may be obtained. À survey of these subjects 
will be given later on. 

Accumulation of emanation near the ground. 
The additional ions or the producers of ions 
necessary to produce the increase of the ionic 
density could not have been conveyed to 
the place of observation from another region 
because the horizontal motion of the air was 
feeble (in some cases calm prevailed). At the 
same time the vertical air motion must have 
been negligibly small, since with a minimum 
temperature at the ground turbulence must 
have been almost completely suppressed. The 
cause of the increased ionic density must be 
found at the place of observation itself. In 
our case the most probable prime cause of 
the increased ionic density might be expected 
to be the accumulation of emanation exhaled 
from the soil into the air layers nearest the 
ground, as shown more detailed in (18), (19). 
This effect may occur only under certain 
meteorological conditions under which mixing 
of the air is suppressed. 

Austausch-processes and electrical phenomena 
in the atmospheric air. The mixing processes in 
atmospheric air may be considered as Austausch- 
processes — an idea introduced by SCHMIDT 
(21). Therefore the process of the accumulation 
of emanation in the layers of air near the 
ground must be considered as a process stip- 
ulated by Austausch. Attempts have been 
made to apply the theory of Austausch to 
several problems of atmospheric electricity. 
The distribution of radioactive substances in 
the free atmosphere has been considered in 
terms of Austausch by Hess and Scumipr (9), 
SCHMIDT (22), PRIEBSCH (20), and Lerrau (13). 
Variations of the potential-gradient have been 
considered in terms of Austausch-theory by 
Brown (5), FRANKE (8), IsRAËL (11), and of 
the electrode-effect by Lerrau (13). 

Our problem of the increase of the concen- 
tration of small ions near the ground under 
weather conditions, during which the accu- 
mulation of emanation exhaled from the soil 
and stratification of ions near the ground may 
be expected, is more complicated, but it is 
evident that this problem must be considered 
as one stipulated by Austausch. Consideration 
of the individual processes which may be of 
importance to the whole problem could not 
give us results of full certainty because of the 
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lack of several data necessary for the calcula- 
tions (for example, the combination processes 
between small- and large ions and nuclei, the 
recombination processes in the track of an 
«-particle etc.). Even so, conclusions of in- 
terest could be obtained by such considerations. 

Some considerations concerning the trend of the 
number of large ions. In order to show the per- 
spectives which may appear through a more 
detailed consideration of the phenomena 
described in this article we will here mention 
some events which may have a connection 
with the observed variations of large ions. 

The minimum prior to sunrise and the 
enormous increase of the concentration of 
large ions after sunrise was at first a surprise 
to us. But we have taken into consideration 
facts known about the diurnal variation of 
condensation nuclei; large ions may be con- 
sidered as a fraction of the condensation nuclei. 
Unfortunately, “in spite of the numerous 
observations of nuclei, there is not yet a 
definite relationship established between time 
of day and number of nuclei” as stated LAnDs- 
BERG (12). BRADBURY and MEURON (3) have 
stated that “The diurnal variation in nuclei- 
density shows a definite maximum in the 
forenoon and in the evening with a pronounced 
minimum prior to sunrise and a midday minimum 
of lesser extent. Such a variation might be ex- 
pected in the normal daily convective cycle.” 
(Underlined by us.) Characteristics similar to 
those found by BRADBURY and Meuron for 
the ground level have been established by 
YUNKER (24) for a level so cm and another 
20 m above the ground. It is seen that the 
nuclei-density at the 20-metre level is on the 
whole higher than that at the ground. That 
the concentration of nuclei at night may be 
much smaller than the day average is also 
shown by Noran and TEDDE (14). A theoret- 
ical consideration of the trend of the concentra- 
tion of nuclei is given by FRANKE (8) in terms 
of Austauschprocesses. The course of the 
concentration of nuclei during night-hours 
calculated by Franke with a minimum and a 
higher maximum shows a similarity to the 
trend of large ions observed by us. We cannot 
accept in whole the mechanism used by 
FRANKE, but it appears that a description of 
the effect observed by us could be given 
with some modification. 

It seems that a detailed investigation of 
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the associate phenomena may give us important 
information about the events considered here. 
In addition it could be mentioned here that 
also variations of sight may contribute to the 
understanding of our problems. BrettcH and 
RıtschL (2) have considered the diminution 
of sight after sunset and a probable explana- 
tion of this effect is given in terms similar to 
those by which the observed marked diminu- 
tion of sight 1—2 hours after sunrise during 
radiation nights is explained. 


A more detailed account of the investigation 
presented here will appear in Arkiv för Geofysik 
published by Kungl. Svenska Vetenskaps- 
akademien, Stockholm. 
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Die Feldkomponentenmiihle 


Ein Gerät zur Messung der drei Komponenten des luftelektrischen 


Feldes und der Flugzeugeigenladung bei Flugzeugaufstiegen 
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Abstract 


Recording the atmosperic-electric potential gradient from an aeroplane it is very useful to 
get all the three field-components. This is the case especially in investigations in greater altitudes 
and in thunderclouds. In this paper an apparatus is described working on the principle of the 
‘Feldmiihle’ (field strength meter) and giving simultaneous recording of the atmospheric- 
electric potential gradient separated in the three rectangular components of the field. Furthermore 
a method is pointed out to separate the field of the own charge of the aeroplane and to record 
this one separately. This enables the ‘Feldkomponentenmiihle’ (field component meter) to 


be used in motor driven aeroplanes. 


1. Einleitung 


Alle bisher durchgeführten Messungen des 
luftelektrischen Feldes in der freien Atmos- 
phäre angefangen von den Freiballonaufstiegen 

is zu den modernen Messungen im Segel- 
oder Motorflugzeug (1, 2, 3, 4, 5,9) beschränken 
sich auf die Messung oder Registrierung der 
vertikalen Feldkomponente. Solange man bei 
luftelektrisch ungestörtem Wetter nur bis zu 
Höhen von einigen Kilometern über dem 
Erdboden vordringt, kann man erwarten, 
daß das luftelektrische Feld vertikal gerichtet 
ist, und somit die Horizontalkomponenten 
fehlen. Bei Untersuchungen von Wolken und 
insbesondere von Gewitterfeldern aber treten 
neben der Vertikalkomponente auch erhebliche 
Horizontalkomponenten auf: So ist man z. B. 
bei der Registrierung der Vertikalkomponente 
in Gewitterfeldern gezwungen, über die 
geometrische Verteilung der Wolkenladungen 
gewisse Annahmen zu machen, um diese 
Messungen überhaupt auswerten zu können. 


SIMPSON und ROBINSON (3) nehmen deshalb 
eine Ladungsverteilung in übereinanderliegen- 
den Kugeln an, während E. Watt (6) die 
Ladungsanordnung in Schichten vorzieht, um 
für die Auswertung der Messungen einfach 
zu übersehende theoretische Potential- und 
Feldverhältnisse zu haben. H. WICHMAnN (7) 
leitet auf Grund der meteorologischen Ver- 
hältnisse in der Gewitterwolke gekoppelt mit 
den  elektrizitätsbildenden Vorgängen eine 
andere Ladungsverteilung ab, die sich nicht 
mehr in so einfache geometrische Formen 
pressen läßt. Es verbleibt also bei der Auswer- 
tung der Messungen ein gewisser Grad vou 
Unsicherheit, der sich durch die gleichzeitige 
Registrierung aller drei Feldkomponenten be- 
seitigen lassen würde. 

Auch bei Schönwetterflügen wird die Hori- 
zontalkomponente des luftelektrischen Feldes 
dann von wesentlicher Bedeutung sein, wenn 
man mit dem Registriergerät Höhen erreicht, 
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die größer sind als die Gipfelhöhe der Gewitter- 
wolken. Hier muß sich die Elektrizitätsströ- 
mung nach den Seiten ausbreiten und somit 
auch eine horizontale Feldkomponente vor- 
herrschen. Dieser seitliche Ausgleich wird nicht 
erst in der Ionosphäre stattfinden (8), sondern 
bereits in tiefer gelegenen Luftschichten, so 
daß man etwa in dem Raum zwischen 15 
und so km Höhe mit dem Auftreten von hori- 
zontalen Feldkomponenten rechnen darf. 

Aus diesen und anderen Gründen wurde 
vom Verfasser bereits 1944 mit der Entwick- 
lung eines Meßgerätes für Flugzeugaufstiege 
begonnen, das die drei Komponenten des luft- 
elektrischen Feldes gleichzeitig registriert. 
Durch dieses Instrument sollte außerdem die 
Möglichkeit gegeben werden, das von der 
Eigenladung des Flugzeuges herrührende Feld 
von dem luftelektrischen Feld abzutrennen 
und gesondert zu registrieren. Denn es ist 
bekannt, daß Motorflugzeuge durch die Aus- 
puffgase auf Spannungen von $—10000 V 
aufgeladen werden, so daß das Eigenfeld des 
Flugzeuges die schwachen luftelektrischen Fel- 
der in der Höhe überdeckt. 

Die praktische Erprobung der Anlage im 
Flugzeug konnte leider nicht mehr durchge- 
führt werden. Wenn hier trotzdem dieses Gerät 
beschrieben werden soll, so geschieht es aus 
dem Grunde, weil auch heute noch die prak- 
tische Flugerprobung in absehbarer Zeit in 
Deutschland nicht gesichert erscheint. 


2. Theoretische Überlegungen zur Feld- 
komponentenmühle 


Das Meßgerät wurde nach dem Prinzip der 
automatisch rotierenden Wilsonplatte gebaut, 
das unter dem Namen Feldmeßmaschine oder 
speziell in der Luftelektrizität als Feldmühle 
bekannt ist. In Verbindung mit der modernen 
Verstärkertechnik hat sich dieses Meßverfahren 
wegen seiner Unempfindlichkeit gegen mecha- 
nische Beanspruchungen und wegen der 
geringen Isolationsanforderung zur Registrie- 
rung luftelektrischer Felder im Flugzeug aus- 
gezeichnet bewährt. 

Das Gerät mißt letzthin die durch das 
äußere Feld influenzierte Ladung auf einem 
Teilstück der Flugzeugoberfläche, bzw. die 
dieser Ladungsdichte proportionale Feldstärke. 
Zum Verständnis der 3 Komponentenmessung 


241 


Abb. ı. Leitender Körper im elektrischen Feld beliebiger 
Richtung. 


soll eine kurze theoretische Überlegung ein- 
geschaltet werden. 

Wir setzen einen beliebig geformten leiten- 
den Körper in ein beliebig gerichtetes elek- 
trisches Feld. Dieses Feld mit der Feldstärke F 
zerlegen wir entsprechend unserem karte- 
tesischen Koordinatensystem x, y, Z in seine 
Komponenten X, Z, Z. (Abb. 1.) 

Weiterhin schneiden wir uns 3 beliebige 
Oberflächenstücke 1, 2, 3 aus dem Körper 
aus, die aber durch einen Draht mit der 
übrigen Oberfläche des Körpers leitend ver- 
bunden bleiben sollen. Auf jedem Oberflächen- 
stück wird von jeder Komponente des äußeren 
Feldes eine Ladung q influenziert, die der je- 
weils betrachteten Feldkomponente propor- 
tional ist. Die den einzelnen Oberflächenstüc- 
ken zugehörigen Ladungen kennzeichnen wir 
durch die Indizees 1, 2, 3 und unterteilen die 
Gesamtladung in jedem Oberflächensegment 
in die Anteile, die durch die Feldkomponenten 
X, Y, Z induziert werden. Diese Teilladungen 
kennzeichnen wir durch weitere Indizees 
X, Y, Z. Es setzt sich also die Gesamtladung qr 
auf dem Oberflächenstück ı zusammen aus 
den Teilladungen qx, qiy, gız. Dasselbe gilt 
entsprechend für die Oberflächenstücke 2 
und 3. Damit erhalten wir folgendes Glei- 
chungssystem. 


Bi = x + qır + nz 
ga = qex + gay + gz (1) 
93 = gsx + gay + gaz. 


Die Teilladungen dx, qiy,.....--- sind 
nun proportional den Feldkomponenten X, 


Ag 


Y, Z. Diese führen wir mit Hilfe der Propor- 
tionalitatsfaktoren ar, br, &,... ein und er- 


halten aus Gl. (1.) 


qi == di X + Dir GZ 
QE RX + b2Y + fb, (2) 
B = a3X at ne b3Y oe GZ. 


Da die a, bc, konstante Fak- 
toren darstellen, die wir etwa als Teilkapazitäts- 
koeffizienten für die entsprechende Feldrich- 
tung definieren und messen können, dürfen 
wir Gl. (2) nach den Feldkomponenten X, 
Y, Z auflösen und erhalten mit neuen Kon- 
Stanten. Ar Dr. Car veuskane en 


x = Aıgı =e Bxqz air Ciq3 
Ye — Asgı t B2q2 r C2q3 (3) 
ie == Asg: + B;q2 SE C3q3. 


Die Konstanten Ar, Br, Ci,....... können 
wir nach bekanntem Rechenverfahren aus 
den a, bi, Gr, gewinnen. Wir wollen 
jedoch diese Gleichungen übergehen, da diese 
Konstanten bei der praktischen Ausführung 
des Gerätes zweckmäßig durch Eichung ge- 
wonnen werden. 

Als Ergebnis halten wir fest, daß es nach Gl. 3 
theoretisch möglich ist, durch Messung der 
Oberflächenladung an drei verschiedenen Stel- 
len der Oberfläche des Körpers, der sich in 
einem beliebig gerichteten homogenen Feld 


Abb. 2. Photographie der Feldkomponentenmiihle. 
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befindet, die Komponenten des äußeren Feldes 
zu bestimmen. Die praktische Ausführung zeigt 
darüber hinaus, daß es keine Schwierigkeiten 
bereitet, die erforderliche Rechenoperation 
durch das Gerät automatisch a zu 
lassen. 


3. Praktische Ausführung der Feld- 
komponentenmühle 


Die Messung der Oberflächenladung bzw. 
der Feldstärke nach dem Feldmühlenprinzip 
soll hier als bekannt vorausgesetzt werden. 
Wir wenden uns gleich der technischen Kon- 
struktion der Feldkomponentenmiihle zu. Abb. 
2 zeigt eine Photographie solch eines Gerätes. 
Die Feldmiihle besteht also aus einem kurzen 
Zylinder, dessen Kopf in 8 voneinander isolierte 
Segmente aufgespalten ist. Hierbei entfallen 
auf den Zylinderdeckel 4 Segmente, wahrend 
das obere Stiick des Zylindermantels ebenfalls 
in 4 symmetrische Teilstiicke unterteilt ist. 
Uber diesen Segmenten rotiert ein etwas 
größerer geerdeter zylindrischer Flügel, in 
dem Ausschnitte derart angebracht sind, daß 
bei der Rotation je 2 gegenüberliegende Seg- 
mente im gleichen Rythmus frei gegeben und 
abgeschirmt werden. Jedes Segmentepaar des 
Zylindermantels ist durch einen Widerstand 
miteinander verbunden, der durch einen Mittel- 
abgriff mit der Flugzeugoberfläche verbunden 
ist. Die Wirkungsweise dieser Anordnung 
wollen wir uns an der Schemafigur Abb. 3 a— 
3d klar machen. 

Es sind hier die 4 Segmente 1, 3 und 2, 4 des 
Zylindermantels im Schnitt aufgezeichnet. Sie 
werden von einem horizontalen Feld influen- 
ziert, dessen Richtung parallel zur Verbin- 
dungslinie der Segmente 1, 3 gegeben ist, und 
das senkrecht zur Zylinderachse steht. Wir 
schen aus Abb. 3b, daß das Segment ı nur 
negative und das Segment 3 nur positive In- 
fluenzladung trägt. Die Segmente 2 und 4 
dagegen tragen nach Abb. 3a auf ihrer linken 
Hälfte negative und auf der rechten Hälfte 
positive Ladung. 

Wir wollen nun den Ladungsfluß zwischen 
den Segmenten ı und 3 bei einer Rotation 
des Flügels verfolgen. Als Ausgangsstellung 
wählen wir die in Abb. 3a abgebildete Lage 
des Flügels, in der die Segmente ı und 3 von 
dem äußeren Feld vollkommen abgeschirmt 
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Abb. 3a und 3 b. Mantelsegmente unter dem Einfluss eines horizontalen Feldes. 


Abb. 3 c. Von den Segmenten 1 und 3 gelieferte Spannung. 


Abb. 3 d. Von den Segmenten 2 und 4 gelieferte Spannung. 
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sind. Dreht sich der Flügel im Uhrzeigersinn 
bis in die Stellung 3b, so fließen unter der 
Einwirkung des äußeren Feldes Elektronen 
über den Widerstand Rı von Segment 3 auf 
Segment 1. Dreht sich der Flügel weiter, bis 
er wieder die Stellung 3a erreicht, so gleichen 
sich die Influenzladungen über den Wider- 
stand Rı aus, es fließen also Elektronen von 
Segment ı auf Segment 3. Bei weiterer 
Drehung wiederholt sich das Spiel von neuem. 
Dieser dauernde Ladungsfluß durch den Wider- 
stand Rı erzeugt über diesem eine Spannung, 
deren Verlauf in Abb. 3c schematisch dar- 
gestellt ist. Diese Wechselspannung wird dann 
wie üblich durch einen Verstärker verstärkt 
und durch einen zwangsgesteuerten Gleich- 
richter gleichgerichtet dem Schreiber zuge- 
führt. Die Amplitude dieser Spannung ist der 


einfallenden Feldstärke proportional. Der 
zwangsgesteuerte Gleichrichter sorgt außer- 


dem dafür, daß auch das Vorzeichen des 
Feldes nicht verloren geht. 


Wir betrachten nun die Wirkung des Feldes 
auf die Segmente 2 und 4. Gehen wir ebenfalls 
von der Stellung 3a aus, so werden sich beim 
Drehen des Flügels in die Stellung 3b in der 
ersten Halbzeit die negativen Ladungen von 
Segment 2 mit den positiven von Segment 4 
ausgleichen, während in der zweiten Halbzeit 
sich die positiven Ladungen von Segment 2 
mit den negativen von 4 ausgleichen. Dreht 
sich der Flügel weiter in Stellung 3a, so er- 
gibt sich derselbe Ladungsfluß nur mit umge- 
kehrtem Vorzeichen. Die Kurve 3d zeigt 
ungefähr den resultierenden Spannungsverlauf 
über dem Widerstand R.. Dieser ist gegenüber 
dem Spannungsverlauf über dem Widerstand 
Rı so in der Phase verschoben, daß durch 
eine mit dem ersten Stromkreis synchrone 
Zwangsgleichrichtung die resultierende Gleich- 


“sea ed 


Abb. 4. Influenz des Segmentepaares 1 und 3 durch ein 
senkrecht einfallendes Feld. 
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spannung gleich o wird. (Die gleichgerichtete 
Spannung ist in Abb. 3c und 3d gestrichelt 
eingezeichnet.) 


Drehen wir jetzt-die Feldrichtung so, daB 
sie mit der Verbindungslinie der Segmente 2 
und 4 zusammenfällt, dann liefern sinngemäß 
diese die Spannung, die wir eben für die Seg- 
mente ı und 3 abgeleitet haben, während 
letztere die gleichgerichtete Spannung o geben. 
Durch den Kunstgriff der zwangsgesteuerten 
Gleichrichtung sprechen also sowohl die Seg- 
mente ı und 3 als auch die Segmente 2 und 4 
nur auf die Feldkomponente an, die parallel 
zu der Verbindungslinie der entsprechenden 
Segmente gerichtet ist. 


Wir müssen uns jetzt noch davon über- 
zeugen, daß diese Segmentpaare auch unemp- 
findlich sind gegenüber der Feldkomponente, 
die in Richtung der Zylinderachse einfällt. 
Zu diesem Zweck betrachten wir Abb. 4. In 
dieser sind die beiden Segmente ı und 3 im 
Aufriß dargestellt. 


Durch das in Achsrichtung einfallende Feld 
wird auf beiden Segmenten die gleich große 
negative Ladung influenziert. Diese Ladung 
fließt dem Segment ı über den linken halben 
Widerstand zu und dem Segment 3 über den 
rechten halben Widerstand. Da in den beiden 
Widerstandshälften die Stromrichtung gegen- 
sinnig ist, kompensiert sich der dadurch ent- 
stehende Spannungsabfall über dem ganzen 
Widerstand zu o. Dasselbe gilt für das Seg- 
mentepaar 2 und 4. 


Voraussetzung für diese Kompensation ist 
die Symmetrie des Feldlinienbildes. Montiert 
man die Feldmühle z. B. auf der Oberseite 
des Flugzeuges, so ist diese Symmetrie für das 
in Richtung der Flugzeugflügel weisende Seg- 
mentepaar immer dann gegeben, wenn man 
die Mühle über der Mittellinie des Flugzeugs 
anbringt. Für das nach Kopf- und Schwanzende 
zeigende Segmentenpaar müßte die elektrisch 
neutrale Linie erst bestimmt werden. Läßt 
sich diese Linie aus irgend einem Grund für 
die Montage nicht verwenden, so kann man 
durch einen anderen Erdungsabgriff am Wi- 
derstand die vollständige Kompensation er- 
zwingen. 


Abschließend können wir feststellen. Die 
Segmentepaare des Zylindermantels geben 
jedes nur für jeweils eine Komponente eines 
äußeren Feldes eine gleichgerichtete Spannung. 
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Abb. 5a. Deckelsegmente unter der Einwirkung eines senkrecht einfallenden Feldes. 
Abb. 5b. Deckelsegmente unter der Einwirkung eines horizontal einfallenden Feldes. 


Gegen die jeweils anderen beiden Feldkompo- 
nenten ist das Segmentenpaar unempfindlich, 
teils durch die zwangsgesteuerte Gleichrich- 
tung, teils durch eine entsprechende Erdung 
des überbrückenden Widerstandes. 


Etwas einfacher zu übersehen sind die Ver- 
hältnisse für die Segmente im Zylinderdeckel. 
Hier genügt es, wenn von den zwei Segment- 
paaren überhaupt nur eines benutzt wird. 
Dieses wird nach Abb. sa kurzgeschlossen 
und über einen Widerstand mit dem Flugzeug 
verbunden. 


Die Wirkung für parallel zur Zylinderachse 
einfallende Felder ist dieselbe, wie bei der 
normalen Einkomponenten-Mühle. Bei sen- 
krecht zur Zylinderachse einfallenden Feldern 
hebt sich die Wirkung entweder schon inner- 
halb ein und desselben Segments oder durch 
das Kurzschließen beider Segment auf, wie 
aus Abb. sb ohne weiteres ersichtlich ist. 

Wenn man nicht Wert darauf legt, daß man 
für jeden Komponentenkanal des Verstärkers 
einen getrennten Eingangswiderstand erhält, 
könnte man auf die Segmente im Zylinder- 
deckel überhaupt verzichten. Denn parallel 
zur Zylinderachse einfallende Felder kann man 
auch mit einer Schaltung nach Abb. 6 durch die 
Zylindermantelsegmente registrieren. 

Hier fließt über den eingeführten Wider- 
stand R; die Influenzladung auf die Segmente 1 
und 3. Die Spannung für die senkrecht ein- 
fallende Feldkomponente wird über dem Wi- 


derstand R; abgenommen. Der Querwider- 
stand bleibt dabei ohne Wirkung. Die Ähn- 
lichkeit dieser Schaltung mit den Ausführungen 
nach Abb. sa und sb ist so evident, daß nach 
den oben gemachten Ausführungen eine ein- 
gehende Diskussion nur eine Widerholung 
bedeuten würde. 

Mit den bisherigen Ausführungen ist ge- 
zeigt, daß bei einer der angegebenen Kon- 
struktionen der Feldkomponentenmühle je- 
weils ein Segmentepaar immer nur auf eine 
Feldkomponente anspricht, während es durch 
die anderen beiden Komponenten nicht beein- 
flußt wird. 


Abb. 6. Registrierung vertikaler Felder mit Hilfe der 
Zylindermantelsegmente. 
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4. Die Registrierung der Flugzeug- 
eigenladung 


Eine Komplikation entsteht bei der Feldregi- 
stricrung im Flugzeug dadurch, daß Motor- 
flugzeuge durch die Auspuffgase so hoch auf- 
geladen werden, daß ihre Eigenfelder die 
schwachen luftelektrischen Felder in der Höhe 
überdecken. Man kann diese Schwierigkeit 
natürlich dadurch umgehen, daß man ein 
Segelflugzeug verwendet, von dem man eine 
etwa vorhandene Eigenladung mit Hilfe 
eines in der neutralen Linie angebrachten Aus- 
gleichers wegbringt. Andererseits wird grade 
bei Messungen in Regen- und Gewitterwolken 
auch das Segelflugzeug durch Einfangen von 
geladenem Niederschlag eine derartige Eigen- 
ladung annehmen, die durch einen Ausgleicher 
nicht restlos beseitigt werden kann. 


Bei der Komponentenmühle bietet sich zur 
Abtrennung dieser Eigenladungsfelder folgen- 
der Weg. Es wird hierzu für die Zylinder- 
mantelsegmente eine Schaltung nach Abb. 6 
und für die Zylinderdeckelsegmente eine 
Schaltung nach Abb. sa verwendet. Wir 
bezeichnen die parallel zur Mühlenachse ein- 
fallene Komponente des äußeren Feldes mit 
Fx und die von der Eigenladung herrührende 
Feldstärke mit E. Für das Eigenladungsfeld 
gelten bei den einzelnen Segmenten andere 
Teilkapazitätskoeffizienten als für das äußere 
Feld. Die ersteren bezeichnen wir mit bı, b2 
und die letzteren mit ar, 42. 


Mit der Schaltung nach Abb. 7 können wir 
jetzt das Eigenladungsfeld von dem äußeren 
Feld auf folgende Weise separieren. 


eb, a fate, 


7 | | | “ = 
ai Q ly 
te. 


fo, Ux (oy 


Abb. 7. Schaltung zur Trennung der Wirkung des 
Eigenladungsfeldes von der des äusseren luftelektrischen 
Feldes. 
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Die über dem Widerstand Ry entstehende 
Spannung U; ist gegeben durch 


U; = ar1Fy + bıE, (4) 
während für die über R; entstehende Spannung 
U; gilt 

U; = aFx + b2E. (5) 

Durch passende Wahl der Widerstände R; 
und R, können wir erreichen, daß ar = a, wird. 
Dann erhalten wir durch Gegeneinanderschal- 
ten der Widerstände R; und Ry an der Punkten 


3 und 4 die Eingangsspannung U, für den 
Verstärkerkanal der Eigenladung. Denn es ist 


bE. (6) 


62 Ue 


Es ist also U. nur proportional der Eigen- 
ladungsfeldstarke E. Dasselbe Schaltbild ver- 
wenden wir fiir die andere Kombination der 
Zylindermantel- und Deckelsegmente. Hier 
gleichen wir die R; und Ry entsprechenden 
Widerstände aber so ab, daß in den Gleichungen 
(4) und (s) by = b: wird. Dann fällt der Eigen- 


ladungsanteil heraus, und wir erhalten 
Ur U; = Us = (ar — a2) Fx. (7) 


Ux ist also eine Spannung, die nur propor- 
tional der parallel zur Mühlenachse einfallen- 
den Feldkomponente Fx ist. 

Eine Stôrung der Anzeige der horizontalen 
Feldkomponenten durch die Eigenladungs- 
feldstirke erfolgt nicht, wenn die Mantelseg- 
mente nach Abb. 3a bzw. 7 geschaltet werden. 
Es gilt dann für das Eigenladungsfeld dasselbe, 
was an Hand der Abb. 4 fiir die vertikale 
Feldkomponente dargelegt wurde. Durch den 
Mittelaberiff des Querwiderstandes fällt die 
Wirkung des Eigenladungsfeldes heraus. 

Eine praktische Schwierigkeit könnte aber 
dadurch entstehen, daß sowohl die Werte für 
a und a; als auch für by und b2 fast gleich 
groß sein werden. Durch die Differenzbildung 
in den Klammern von Gl. 6 und 7 erhalten 
wir also für U, und Ux kleine Werte. Nach 
dem Stand der heutigen Verstärkertechnik be- 
reitet es zwar keine Schwierigkeit auch kleine 
Wechselspannungen genügend hoch zu ver- 
stärken. Aber schwerer wird es sein, den 
genauen Abgleich der Widerstände zu errei- 
chen. Außerdem muß die von der Kontakt- 
spannung und die von der mehr oder minder 
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guten Erdung des rotierenden Flügels herrüh- 
rende Störspannung der Mühle klein gehalten 
werden. 


Sollte sich also herausstellen, daß dieser Weg 
nur mit großem Aufwand beschritten werden 
kann, so bleibt immer noch die Möglichkeit 
offen eine einfache Hilfsmühle mit hinzuzu- 
ziehen. Diese kann dann an einem geeigneten 
Ort und zwar an der Unterseite des Flug- 
zeuges montiert werden. Dadurch erhalten die 
entsprechenden Koeffizienten &r, a2 und br, bz 
verschiedenes Vorzeichen, verstärken also 
einander. Für die Konstruktion der Kom- 
ponentenmühle gelten sonst sinngemäß die- 
selben Erfahrungen, die man beim Bau der 
einfachen Feldmühle gesammelt hat, auf 
die aber hier nicht mehr eingegangen wer- 
den soll. 


Zusammenfassung 


Bei der Registrierung des Luftelektrischen 
Feldes im Flugzeug ist es äusserst wünschens- 
wert, alle drei Feldkomponenten aufzuneh- 
men. Dies gilt besonders bei Untersuchungen 
in grossen Höhen und bei Gewitterfeldern. 
In dieser Arbeit wird eine Apparatur be- 
schrieben, die nach dem Feldmühlenprinzip 
arbeitet und die gleichzeitige Registrierung 
der nach kartesischen Koordinaten zerlegten 
Komponenten des luftelektrischen Feldes ge- 
stattet. Ferner wird ein Weg aufgezeigt mit 
dieser Apparatur das Feld der Eigenladung 
des Flugzeuges von dem luftelektrischen Feld 
abzutrennen und für sich gesondert zu re- 
gistrieren. Dadurch ist es möglich die Feld- 
komponentenmühle auch im Motorflugzeug 
zu verwenden. 
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Numerical Tendency Computations from the Barotropic 


Vorticity Equation 


B. BOLIN: and J. CHARNEY, Institute for Advanced Study, Princeton 


(Manuscript received ro October 1951) 


A bs tract 


As a preliminary test of the usefulmess of the barotropic vorticity equation for numer- 
ical forecasting a series of tendency computations have been made for a period of tn 
days in February ı951. Despite the difficulty of comparing Instantaneous tendencies with 
observed changes over a finite time interval (r2 hours) some positive information was 
obtained concerning computational errors as well as errors im the atmospheric model. 
The computations support the conclusion reached by CRaRNEY and Errassen (1049) that 
the Rocky Mountains are of importance in their influence on the large wale upper level 


flow pattern over the United States. 


I. Introduction 


In a recent article by CHARNEY, FJORTOFT 
and von NEUMANN (1950): a series of numerical 
forecasts for the soo mb surface was made 
using the barotropic vorticity equation. Un- 
fortunately, the number of forecasts was too 
small to permit general conclusions about the 
accuracy of the computations, nor was it 
possible to separate truncation errors from 
errors due to defects in the model. Final con- 
clusions of this kind cannot be drawn until a 
large number of complete 24-hour forecasts 
have been made. Such computations are plan- 
ned for the high-s electronic computer 
being built at the Institute for Advanced 
Study. 

In the four complete forecasts that were 
described in (1), it was found that the initi- 
ally computed tendency gives a fairly good 
idea of = 24-hour change obtained by a 
complete integration. This is obviously still 


ı Now at University of Stockholm. 
2 In the following we shall refer to this article by (1). 


more true for a 12-hour forecast. Thus one 
might expect to obtain some useful results 
from a comparison between initial tendencies 
and the subsequent observed 12-hour change 
provided some proper way of co i 


can be devised. With this idea in mind a 


mom ean 


series of tendency computations was carried | 


through for a period in February 1951. It was 
possible to obtain some information from 
these computations regarding computational 
errors as well as defects in the model. It 
was also instructive to study the computed 
tendencies in some detail in that they furnish 
direct illustrations of the quasi-barotropic be- 
haviour of the atmosphere. 

It should be stressed that a series of com 
forecasts has to be made as a final test of the 
computational scheme as well as of the model. 
We can obviously not obtain any information 
from these computations about errors due to 
the finite extrapolation in time and the extent 


to which computational errors put a limit on| 


the time over which a forecast can be made. 


COMPUTATIONS FROM THE BAROTROPIC VORTICITY EQUATION 


This study should rather be looked upon as 
a series of experiments intended to serve as a 
guide for further work. 


2. The solution of the vorticity equation 


The reader is referred to (I) for a detailed 
discussion of the problems that arise when 
trying to solve the barotropic vorticity equa- 
tion for a certain initial set of data. Here only 
a general outline is given of the procedure 
followed. 

Both from empirical and theoretical in- 
vestigations we know that there exists a level 
in the mid-troposphere which is approximately 
non-divergent. As in (I) we shall assume that 
the soo mb level coincides with this level. (It 
will be seen that some of the major errors in 
the tendency computations are introduced by 
this assumption.) Thus the vorticity equation 
may be written 


ae A 
NEN ET fi): (1) 


€ is the relative vorticity, v is the horizontal 
wind velocity and f the Coriolis parameter. 
It was shown in (I) how this equation may 
be transformed into 


v2 2 =5( Vz +62) =] (7,2), 
(2) 


where the geostrophic wind equation has been 
used for the evaluation of the vorticity. Here 
z denotes the height of the soo mb surface, 
m is the scale factor in the transformation from 
spherical coordinates to the particular map 
used in the computations, 7 is the absolute 
vorticity and J (&:, %2) is defined as J (a1, #2) = 
dax |Ox + d&2/0y — Oa2/Ox + Oxx/dy. It is our in- 
tention to solve this Poisson equation in order 
to obtain the initial value of 9z/9f and compare 
that tendency with the observed 12-hour 
change of z. For hand computation equation (2) 
is most rapidly solved by the relaxation meth- 
od described by SOUTHWELL (1946). However, 
it is necessary to specify 92/dt on the bound- 
ary. Since these quantities are not known, we 
are forced to introduce artificial boundary 
conditions. When carrying out a complete 
24-hour forecast, this question deserves special 
attention, but as we only are going to compute 
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the initial tendency we may assume that 
dz/dt = o on the boundary. This will undoubt- 
edly introduce errors close to the boundary, 
but the magnitude of these errors decreases 
rapidly in going away from the boundary. 
The reason is that the effect of the correct 
boundary conditions at a fairly small distance 
from the boundary becomes equivalent to 
dz/t = o on the boundary because of the 
interference of positive and negative boundary 
values. In the verification, no points have been 
used at a distance less than 2 AS from the 
boundary, AS being the grid interval used 
in solving the finite difference equation cor- 
responding to (2). 

For the solution of equation (2) we choose 
a rectangular area on the map with sides 
p : AS and q- AS, where AS is the distance 
between the gridpoints in a square grid and 
p and q are integers. The coordinates of the 
grid points are given by i- AS and j: AS 
(65, Fgh P fOr lyn 2). eG) SAS. Come 
puted for every other point (i +j=2n) 
according to the formula 


I 
(22); = (asp [RE eee 
epee epee | G) 
3 
i+j=2n, 


Di pie I Le eat 


Ji; is then obtained for points where i + j = 
= 2n + I according to 


I 


Ji (n, 2) cs 4 (45): (eee) x 
Nr ee N ee 
len) (4) 
iFj=2n+7 


eier, ee) 
gm? 
fü 
the non-homogeneous term of equation (2) 
is obtained for points on the map forming 
a square grid with the grid interval AS via! 
The equation is solved by relaxation methods. 


(7?2); + fj. Thus 


and where nj = 


250 


This procedure differs from (I) in that the 
vorticity is computed using a gridsize ers 
while equation (2) is solved using a grid in- 
terval of AS V2. The scheme outlined here 
was chosen for the following reasons: In (I) 
it was found that the value AS & 736 km 
was too large, in particular for the evaluation 
of the vorticity. Here AS was chosen consider- 
ably smaller for the computation of the 
vorticity (A S = 4° long. at lat. 45° = 315 km). 
However, the time required for solving equa- 
tion (2) by relaxation methods increases very 
rapidly with increasing number of grid points. 
In order to keep this time within reasonable 
limits, the relaxation was carried out for the lar- 


ger grid size (AS V2). The method cannot 
be used for a complete forecast, because tend- 
encies are only obtained for every other 
point and new vorticities cannot be computed 
over the smaller grid after the first time ex- 
trapolation has been carried out. Finally it 
might be of some interest to mention that a 
complete calculation of dz/ot takes four to 
five hours for an experienced computer for 
a grid with p = 18 and q = 15. This includes 
the evaluation of the initial field of z from 
the map. 


3. Period under investigation and methods 
of verification 


In all about forty computations were made 
in which tendencies were computed for 
every twelve hours for the period February 
3, 03 GMT—February 13, 03 GMT over an 
area covering the United States, the major 
part of Canada and the western part of the 
Atlantic Ocean. For these computations the 
conventional analysis from the Weather Bureau 
—Air Force—Navy Analysis Center (WBAN) 
was used. It was pointed out in (I) that an ob- 
jective analysis would have been preferable and 
the same conclusion may be drawn from the 
present computations. It was also realized that 
the data coverage over the oceans is too poor 
to permit an analysis that is accurate enough 
for this type of computation. For this reason 
it was decided only to verify the computations 
over the United States and southern Canada. 
In order to see to what extent the analysis may 
have influenced the results over these areas, 
the maps February 4, 15 GMT—February 12, 
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15 GMT were re-analysed without reference 
to the WBAN analyses using smoothness in 
the vorticity field as a basic guiding principle. 
The same computations then were carried 
out for these re-analysed maps. 

Some simple statistical methods have been 
used for the estimation of the errors in the 
computations. From the two independent 
analyses the probable error, due to uncertain- 
ties in the analysis and the replacement of the 
differential equation (2) by a finite difference 
equation was computed. This involved the 
assumption of a Gaussian distribution of errors. 
For the comparison of the computed tendencies 
with actual observed changes over twelve hours 
several methods can be used. For the present 
study the computed value of 1/2 [(2z/d6),=;, + 
+ (dz/At)t=1,4+4| * At = x1, was compared 
with the observed value of (2; 1,4 4 — Zt=t,) = 
= x2 (At being 12 hours). This method has the 
advantage that no assumptions concerning the 
movement of the fall and rise areas have 
to be made, such estimates being always 
more or less subjective. On the other hand, 
it has the disadvantage that two compu- 
tations are involved in each verification. 
Furthermore, this verification scheme cannot 
be used if the changes are quasi-periodic with 
a period less than 36 to 48 hours. In such cases 
all changes will be underestimated. During 
the period under consideration, this has not 
affected the computations to any considerable 
extent. 


The comparison between xı and x: was 
made in a few different ways. The following 
quantities were computed: 

t) The correlation coefficient r 

2) The regression equation x1 = ax2 +b 


3) A quantity K defined by 


VZ (xx = X2)? 
VE Xi? a Vd x2? 


K= 


K has the advantage of being a function both 
of r and a. K=0 if xy; = x; for all i, and 
is approximately 0.7 if r = 0, the exact value 
depending upon ox and 62. (o1 and o being 
the root mean square of x1 and x2 respectively.) 
Quite generally K < 1. 

ets should be stressed that this study is not 
intended to be a complete statistical treatment 
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of the computed data. It contains merely 
some order of magnitude estimates of errors 
in order to permit the discussion of physical 
aspects of the problem. 


4. Computational errors 


Apart from errors due to defects in the 
barotropic model the following errors arise 
in the computations: 


1) Errors due to uncertainties in the analysis 
of the initial data. 


2) Truncation errors, e.g. errors due to the 
replacement of the differential equation by a 
difference equation. 


3) Errors depending upon the particu- 


lar method of verification in replacing 
to + At 


J (dz/dt) dt by xx. Furthermore, there are 
to 
inaccuracies in the determination of the ob- 


served change, x2. 


By comparing the tendencies that were 
obtained from the two analyses of the same 
synoptic situations we can obtain an estimate of 
the combined effect of 1) and 2). It was found 
that the probable error of the tendency com- 
putation introduced by these two factors was 
100 feet/12 hours compared with a root- 
mean-square of the tendencies themselves of 
370 feet/12 hours during this period. These 
errors may be considered to be distributed at 
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random, thus giving rise to a probable error 

of 70 feet in the determination of xı. From 

the same two sets of maps we find that the prob- 

able error in x2 is 30 feet. We cannot obtain a 

direct estimate of the errors that are introduced 
to+ At 


by replacing / (dz/at) dt by xx. In spite 
to 


of this it seems safe to conclude that the 
observed probable difference of 160 feet 
between xı and x» cannot be explained by 
computational errors. This is supported by 
other evidence as shown below. 


It is now of some interest to examine the 
agreement between xı and x: somewhat more 
closely. The quantities r and K were computed 
for an area covering the United States and 
southern Canada for every tendency calcula- 
tion that was made. In fig. 1 these values 
are given as functions of time. The solid 
curve represents the results that were obtained 
by using WBAN-analysis, the dashed curve 
is based on computations from the re-analysed 
maps. The two curves follow each other fairly 
closely except for the last few points of the 
diagram. This again indicates that factors other 
than analysis and the method of finite dif- 
ferences influence the results, which after all is 
very likely as we know that the atmosphere 
is not barotropic. The differences that are 
observed at the end of the period are mainly 
due to changes in the analysis on Feb. 11, 15 
GMT and Feb. 12, 03 GMT. These were con- 
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Graphical representation of the variation of the correlation coefficient, r (a), and the verification co- 
The solid lines represent the results obtained from the WBAN-analyses, the dashed lines those 


obtained from the reanalysed maps. 
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Fig. 2: Computed mean tendency, x,, in tens of feet/r2 hours, for the period 8 Feb. 03 GMT—15 GMT, 
1951 (a), and the corresponding observed change, x, (b). The situation was chosen to picture an average case 
(= 072, IS SO) 


siderable, in particular on the later occasion. 
It should be pointed out that when the number 
of observations is as large as it is over the 
United States, changes of the analysis to such 
an extent can be made only if some observa- 
tions are neglected in one analysis but not in 
the other. 

The mean values of r and K for the whole 
period are 0.74 and 0.43 for the WBAN-maps 
and 0.77 and 0.41 for the re-analysed maps, a 
scarcely significant difference. Itisnotimmedia- 
tely clear what these values mean, and we can 
not say what r and K would be for complete 
24-hour forecasts. The iterative method used 
for that computation certainly introduces 
errors. Tentatively we may say that the values 
of r and K obtained here are fairly close to the 
values that would be obtained in a 12-hour 
forecast but somewhat larger than one might 
expect for a 24-hour integration. In order to 
give an idea of what the values of r and K mean 
the result of a computation is shown in fig. 2. 
The situation was chosen to picture an average 
case (r = 0.73, K = 0.41). 

Returing to the curves in fig. 1 we next 
ask ourselves for the reason for the systematic 
variation in time of r and K that we observe. 
A closer analysis quite clearly shows that the 
relative errors of the computed tendencies 
decrease when the tendencies themselves be- 
come larger, in other words the results are bet- 
ter when well-defined rise and fall areas exist. 
Thus the computations were most successful 
on Feb. 6 and 7 when waves in the westerlies 


moved along the jet in a regular way. Be- 
sides this effect certain physical factors prob- 
ably have given rise to these variations. The 
material is, however, insufficient for a closer 
examination of these aspects. 


5. The barotropic model 


In this whole approach to the problem of 
forecasting large-scale atmospheric flow patterns 
the two most important questions are: 


1) To what extent does the atmosphere 
behave as a two-dimensional barotropic fluid? 


2) How is the theory for a barotropic 
atmosphere best applied to the real atmosphere? 


These questions have been discussed intensive- 
ly during the last ten years since Rossby and 
collaborators presented the first ideas on long 
waves in the atmosphere (1939). From the 
results obtained in (I) and in the computations 
carried out here we can make some further 
comments on the importance of barotropic 
processes in the atmosphere. 

It is well-known that the rotation of a fluid 
has a pronounced influence upon its motion. 
This was already noticed by Taytor (1921) 
who was able to demonstrate this effect in 
simple model experiments. The stronger the 
rate of rotation the more pronounced is the 
tendency+for the motion to become quasi- 
similar in planes perpendicular to the axis of 
rotation. It seems reasonable to assume that an 
analogous effect is present in the atmosphere. 
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Here the motion is already kinematically 
restricted to take place in a hemispherical shell 
and for large-scale flow patterns we may thus 
expect that it is the vertical component of the 
rotation of the earth that is of importance. 
Consequently this tendency for similarity from 
level to level should be found at high and 
middle latitudes but not in equatorial regions. 
By and large this agrees with what is observed 
in the atmosphere. In middle latitudes the 
similarity from layer to layer is a striking 
feature of the flow (above the friction layer), 
while recent investigations of the tropics show 
large differences between the lower and upper 
troposphere (RIEHL, 1950). This gives some 
justification for treating the motion of the 
atmosphere in the middle latitudes as two- 
dimensional. 

The atmosphere is a thermally active fluid 
and at first it may seem strange that worth- 
while results can be obtained by treating it as a 
barotropic fluid. However, we know that it 
would take approximately one weck for fric- 
tion to reduce the kinetic energy of the atmos- 
phere by so per cent if this loss were not 
replaced by baroclinic processes. The total 
kinetic energy is approximately constant, thus 
the baroclinic processes in the atmosphere do 
not replace more than about 10 per cent of 
the total kinetic energy in twenty-four hours. 
There are several reasons to believe that this 
change of potential energy into kinetic energy 
does not take place simultaneously all over the 
hemisphere but mainly in connection with 


Fig. 3a. 


soo-mb chart at 03 GMT, 6 Feb. 1951. Thin lines are the contours labelled in tens of feet. The heavy 
lines in fig. a give the distribution of absolute vorticity in units of 10-5 sec? and the dashed lines in fig. b 
show the computed tendency in tens of feet/12 hours. 
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what is known as real deepening or intensifica- 
tion of atmospheric systems. Thus outside such 
restricted areas the change in the atmospheric 
flow patterns are primarily the result of a 
redistribution of kinetic energy. Such a redistribu- 
tion essentially is a barotropic process. 

Even if the total kinetic energy of a baro- 
tropic fluid is conserved (apart from frictional 
dissipation) local intensifications of the flow 
may occur in a barotropic fluid. It has been 
shown (Kuo, 1949; FJORTOFT, 1950) that 
waves in parallel flow are unstable if certain 
criteria are fulfilled and Yen (1949) has 
pointed out that dispersion of kinetic energy 
may give rise to intensification of disturbances 
in the flow. Some examples have also been 
given indicating that these processes are of 
importance in the atmosphere. In most cases, 
however, it is difficult to tell whether baro- 
clinic processes have not also affected the de- 
velopment. Having been able to compute the 
changes in the fluid that would take place if 
the atmosphere were barotropic, we may ob- 
tain an idea of where and when the baro- 
clinicity plays a significant role. From this 
point of view it is of some interest to study 
the series of computations from Feb. 6, 
03 GMT, to Feb. 8, 15 GMT. During this 
period a wave developed from a fairly small 
disturbance in the westerlies to a very intense 
trough. The soo mb maps for Feb. 6, 03 GMT, 
Feb. 7, 03 GMT and Feb. 8, 15 GMT are shown 
in figs. 3a to sa. The heavy lines give the 
distribution of absolute vorticity. The com- 
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Fig. 4a. 


Fig. 4b. 


Fig. 4: soo-mb chart at 03 GMT, 7 Feb. 1951. Legend as in fig. 3. 


uted tendencies are given in figs. 3b—sb 
(dashed lines). We observe that the changes 
during the first 24 hours seem mainly to 
be the result of a redistribution of the vor- 
tex lines, and there is only a slight increase 
of the intensity of the vorticity maxi- 
mum in spite of the fact that the wave is 
amplifying. This indicates that the changes are 
quasi-barotropic (the diagram in fig. 1 shows 
that the most successful tendency computa- 
tions for the whole to-day period were made 
during these 24 hours). It is now interesting to 
see that the development of the wave, at least 
qualitatively, is indicated by the barotropic 
tendencies. Both on Feb. 6, 03 GMT and Feb. 
7, 03 GMT (in fact also on Feb. 6, 15 GMT, 
not shown here) the zero-line of the tendency 
field is located well to the west of the trough 
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line as computed barotropically and the fall 
area in front of the trough is more intense than 
the rise area behind (figs. 3 b, 4b). On Feb. 8, 
on the other hand, we observe a pronounced 
increase of the vorticity in the trough. It is 
also seen from fig. 1 that the tendency com- 
putations become worse during this period 
and a complete 24-hour forecast probably 
would give still larger discrepancies between 
the computed and observed changes of the flow 
pattern. Here baroclinic effects seem to have 
played an important role in the development. 
From this qualitative discussion one therefore 
gets the impression that the initiation of the 
development of the trough was barotropic 
while baroclinic effects became more and more 
important for the further intensification. 

In fig. 3 another illustration of the quasi- 
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soo-mb chart at 15 GMT, 8 Feb. 1951. Legend as in HONTE 
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barotropic behaviour of the atmosphere is 
shown (northern US around long. 100° W). 
We see clearly how the contour lines of the 
soo mb surface are ajusted when a strong jet 
pushes into a region where the flow is less 
intense. The intensification of the gradient 
is clearly indicated in the tendency computa- 
tions and agrees by and large with the observed 
changes. 

By inspecting a series of barotropic tendency 
computations of this kind it is possible to 
formulate qualitative rules for the changes 
that may be expected from different flow 
patterns. This is not the place for a detailed 
discussion of such rules; it should merely be 
stated that many of them are in accordance 
with empirical rules that have proven to be 
useful for the daily prognosis of the flow 
pattern at upper levels. 

Up to this point very little has been said 
about the best methods for applying these 
barotropic tendency computations to the real 
atmosphere. We have simply assumed that the 
flow at soo mb is governed by the barotropic 
vorticity equation, in other words that the 
soo mb surface is the level of non-divergence. 
This is not in accordance with reality, and 
considerable errors are introduced by the 
assumption. Some idea of these errors were 
obtained in the following way: For cach 
gridpoint in the area of verification the regres- 
sion coefficient a was computed using all 
tendency computations that were made. Sys- 
tematic variations of a were found as shown 
in fig. 6. During this period a had a maximum 
value of approx. 1.4 in the Great Lake Region 


Fig. 6: Geographical distribution of the regression co- 
efficient, a, computed from the whole material. 
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and was less than 0.6 over parts of the Rocky 
Mountains. This indicates that the changes as 
computed for an equivalent barotropic atmos- 
phere were stronger than the observed changes 
at soo mb in the Middle West, while the reverse 
was true over the mountains. To some extent 
this may have been caused by the use ofthe geo- 
strophic assumption which results in an over- 
estimate of the vorticity advection in troughs 
and an underestimate in ridges. During this 
period a mean trough was located over the 
middle west where the computed tendencies 
were too strong and a ridge is found over 
the Rocky Mountains. However, the dif- 
ferences are too large to be explained merely as 
an effect of the geostrophic approximation. In- 
stead we may look upon them as indications of 
a systematic variation in the height of the level 
of non-divergence, the level lying below the 
soomb surface where a > 1 and above where 
a< 1. Thus the lines of constant a as a first ap- 
proximation may be considered as contour lines 
of the level of non-divergence. It has been 
pointed out by CHARNEY (1947) and others that 
such variations of the height of that level may 
be expected in the long waves in the westerlies. 
In this particular area the influence of the 
Rocky Mountains also seems very important. 
As the level of non-divergence is located 
approximately in the middle of the atmos- 
phere (with respect to pressure) one might 
expect to find it at higher elevation over the 
mountains. As the computations were fairly 
restricted both in time and space we were not 
able to examine more closely the cause of the 
varation of the height of the level of non- 
divergence. The values of a as computed here, 
however, strongly indicate that the assumption 
of the soo mb level as the level of non-diver- 
gence introduces considerable errors. In (I) it 
was indicated that integrating through the ver- 
tical and thus treating some kind of a mean 
atmosphere tends to remove this difficulty. 
However, the interpretation of the flow pat- 
tern with respect to weather then becomes 
more difficult. 

Another possible effect of the mountains 
is indicated in these computations. It was 
found that the constant factor b in the re- 
gression equation has a very characteristic 
distribution over the area under consideration. 
(cf. fig. 7). This may be explained as a result 


of neglecting the divergence and convergence 


Fig. 7: Geographical distribution of the constant fac- 
tor, b, in the regression equation, computed from the 
whole material (tens of feet/12 hours). 


that is caused by the mountains when the flow 
is forced over them. We have made an estimate 
of the mean 12 hour change of the soo mb 
surface caused by this effect during the period 
by assuming that the divergence term in the 
vorticity equation may be computed by putting 
div: V = Vo : Vh, where vo is the surface wind 
and his the height of the mountains. It was as- 
sumed that vo = 1/2 v and the tendencies thus 
obtained are shown in fig. 8 expressed in 10 
feet/12 hours. The resemblance between fig. 7 
and fig. 8 is quite good both with respect to 
the structure of the pattern and the magnitude 
of the changes, in particular as the assumption 
that vo = '/2 visa fairly crude approximation. 
By and large the computations support an 
idea already put forward by CHarney and 
ELIASSEN (1949). 

The disagreement between the observed 
changes of the flow pattern and those de- 
rived from the tendency calculations here 
have been interpreted as a result of the varia- 
tion of the level of non-divergence and as an 
effect of the mountains. 

A few attempts were made to correct the 
barotropic computations in a manner indicated 
in (I) and in CHARNEY (1951) by assuming that 
the thermal wind rather than the actual wind 
does not vary in direction with height and that 
temperature is advected horizontally. The 
changes in the soo mb tendency were small 
and gave no systematic improvement. The 
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Fig. 8: Computed tendency field (in tens of feet/12 
hours) as a result of horizontal divergence and con- 


vergence caused by the mountain using !/, v as re- 
presentative for the mean surface wind during the 
period. 


surface tendencies, which were also calculated 
by the method, were for the most part of 
the right sign. We were, however, led to 
the conclusion that the advective hypothesis 
is not satisfactory as the correlation between 
the computed and observed patterns was 
poor. 

Summarizing the results of these computa- 
tions the following should be noticed: 

1) The data coverage over the United 
States seems sufficient to give the initial condi- 
tions with the accuracy that is required for the 
present type of computation. Over the oceans 
a much denser net of upper observations must 
be established before successful forecasts can 
be made by the present methods. 


2) A severe defect of the model used in 
these computations seems to be the assump- 
tion that the soo mb level coincides with the 
level of non-divergence. 


3) Computations of this kind are important 
in that they give us a method of separating 
purely barotropic effects from baroclinic pro- 
cesses. This will be of great help when trying 
to extend the present model of the atmosphere 
to incorporate the baroclinicity. 
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The Structure of Mean Circulations 


By R. W. JAMES, Edinburgh 


(Manuscript received 10 August 1951) 


Abstract 


A discussion of the best mode of expressing the facts of the mean westerly flow in temperate 
latitudes parametrically leads to the result that over North America the shape of the vertical geo- 
strophic westerly wind profile is nearly invariant with respect to time and place. The vertical 
wind profile can be expressed in terms either of the change pressure-gradient with height, or of 
contour-height gradient with pressure-level. In either case the profile shape is characteristic, 
and varies little with season, longitude or latitude over North America. 

The intensity of the mean circulation varies rythmically with season, but that for any one 
month may vary markedly from year to year. 

The profile shape for an oceanic region differs from that of a land region such as eastern 
North America in exibiting a lower ratio of high-level to surface wind, and in being much 
more variable from month to month. 

It is shown that the pressure-difference at 7 km between two latitudes, or the difference in 
mean height of the 200 mb surface gives a measure of the total geostrophic momentum between 
those latitudes. The square of these quantities also gives a measure of the total kinetic energy 
of the westerlies between these two latitudes. 


1. Introduction 


The accumulating mass of aerological mate- 
rial over the past ten years has permitted the 
construction of mean meridional cross-sec- 
tions of geostrophic westerly winds. This is a 
graphic method of conveying the salient char- 
acteristics of the mean circulation pattern 
in some areas of the globe, but it should not 
be regarded as the ultimate in descriptive 
aerology because it happens to be the mode 
best adapted for textbook presentation. 

The present paper will be devoted to an 
alternative, parametric representation of ge- 
neral geostrophic westerly flow in selected 
areas. The scientific advantage of parametric 
representation is that it describes mean flow in 
quantitative terms. This is obviously necessary 
if we are ever to have a quantitative theory 
of the general circulation, for a quantitative 
theory requires quantitative facts to fit it. 


The two methods are essentially alternative 
presentations of the same facts. From a com- 
plete parametric presentations, one may read- 
ily derive a corresponding cross-section, and 
alternatively, a complete set of parameters 
can be inferred from an accurate cross-sec- 
tion. 

We shall find, in passing, that an investiga- 
tion of parameters apt for describing mean 
flows leads to considerably simplified tech- 
niques which may be employed in deriving 


cross-sections. 


2. Mean geostrophic circulation as a 
function of latitude 


If we plot the mean pressure at a fixed 
height at points along a meridian against the 
cosine of the latitude, we get a curve the slope 
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of which is at each point proportional to the 
Westerly geostrophic momentum per unit 
volume, for, if 


p =f (cos y), 
1 ape TI op 
2Rasnp Ip 22a d cos p 


eu, = 


where Q is the earth’s rotational speed, a the 
earth’s radius, @U, the density of westerly 
geostrophic momentum. 

When o, the density, is known, the value 
of U,, the geostrophic wind-speed may be 
inferred. 

Hess (1948) has published summer and win- 
ter mean ascents for stations close to meridian 
80° W, ranging from Arctic Bay, latitude 
73° 16’ N, to Salinas, 2° 12’ S. From these 
data the winter pressure-profiles at heights 3 
km (lower, curve) and 10 km, are plotted as a 
function of cos @ in fig. 1. 

It will be noticed that along much of their 
length, say between the latitudes of Detroit 
and Monsonee (42° and 56° respectively), the 
two curves are almost parallel. We find the 
same shape of curve at other levels between 3 
and 10 km. Below 3 km the profile slope 
falls off. We infer from this parallelism 
between about 40° and 50° that momentum is 
approximately invariant with height in the middle 
and upper troposphere in winter over Eastern 
America. 

At latitudes lower than about 40° the two 
curves begin to diverge, the slope of the upper 
curve being the greater. Hence below about 
latitude 40° in winter momentum increases 
with height in the troposphere. 

We remark that the slope of the 3 km 
curve increases steadily with decreasing 
latitude up to between Huntingdon and 
Charleston, say in latitude 36°. This consti- 
tutes a latitude of maximum momentum and 
is termed the jet latitude. For lower latitudes 
still the momentum decreases, and becomes 
zero at the maximum of the pressure profile 
curve, in latitude 20° approximately. This 
point we shall term the doldrum latitude. 

Examining the 10 km profile we find a 
main jet latitude, again at 36°, and a subsid- 
iary, weak one between Monsonee and Port 
Harrison, say at latitude 55°. The doldrum 
latitude of this level is 14°. 

The locus of the doldrum latitudes at all 
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Fig. 1. Pressure at 3 k. m. (starred points) and at 10 km 
(circles) as a function of Cos y, meridian 80° W. approx., 
winter. 


levels is termed the doldrum axis. It will be 
observed that the doldrum axis slopes equator- 
ward with height, being in latitude 30° at the 
surface, 20° at 3 km, and 14° at 10 km. 

The jet axis is the locus of the jet latitudes 
at all levels. The main jet axis is practically 
vertical in latitude 36°. The secondary jet, 
in latitude 55° does not reach down to the 
3 km level. 

As we shall see in the next section there is a 
tendency for momentum to increase with 
height up to a certain level, and then to 
decrease. The height of maximum momen- 
tum is termed the momentum jet height, and 
the locus of momentum jet heights for all 
latitudes may be called the momentum jet 
profile. The momentum jet height is found 
to be approximatly 7 km. The intersection 
of the jet axis with the momentum jet pro- 
file gives the point at which momentum is 
a maximum. It is termed the momentum jet 
centre. The coordinates of the main momen- 
tum jet centré are p= 36°, 2 = 7 km. For the 
subsidiary jet we have, p = 55°, z=7 km. 

The momentum jet height is not the same 
as the velocity jet height; on account of the 
variation of density with height the maximum 
of velocity is found on the jet axis at a height 
of about 12 km. In describing the velocity 
field we define a velocity jet profile and a velo- 
city jet centre in a way exactly analogous to 
the corresponding momentum parameters. 

In order to characterize the westerly geo- 
strophic flow completely we must finally 
specify an intensity factor. A number of 
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possibilities exist. We could take the momen- 
tum at the main momentum jet centre (or 
velocity at the velocity jet centre). Alterna- 
tively a mean momentum between specified 
latitudes might be used, the best level to take 
being the height of the momentum jet profile, 
7 km approximately. This second alternative 
is more in accord with current usage, which 
accords a unique position to the zonal index 
or pressure difference at fixed level (usually 
the surface) between latitudes 30° and 60°. 

The zonal index at the jet level, 7 km, 
will be shown later to represent the mean 
momentum of the westerlies at all levels 
between the assigned latitudes. 

Once the zonal index at 7 km (or the mo- 
mentum at the momentum jet centre), the 
momentum jet profile, the jet axes and the 
doldrum axis have all been specified, together 
with the variation of momentum with height, 
the structure of the momentum field is com- 
pletely specified. 

Fig. 2 shows the corresponding profiles 
at 3 and 10 km for Hess’ summer data. The 
salient feature is the reduction of the zonal 
index as compared with winter. The pressure- 
drop between Charleston and Port Harrison 
at 7 km is 23.5 mb in summer, as compared 
with 42.1 mb in winter. Mean westerly 
summer momentum is just over onc-half 
the mean winter momentum. 

The doldrum axis has migrated to about 
28° latitude, and is almost vertical, while the 
momentum jet centre is now at latitude 55°. 


3. The vertical profile of momentum 


It has been shown (JAMES unpublished) 
that the latitudinal pressure-difference between 
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Fig. 2. Meridional summer profiles of pressure at 3 km 
and to km, 
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Fig. 3. Pressure at a fixed height between different lati- 

tudes. Top curve, Charleston—Sault-Ste. Marie, inter- 

mediate curve, Sault-Ste. Marie—Port Harrison, lowest 
curve, Charleston—Miami. 


two points in nearly the same longitude varies 
with height according to the relation 


Ap = (Ap)m exp[—k(=—Zm)?] () 


where Z,, is the momentum jet height, and 
(A p)m is the pressure difference at this height. 
We may term (Ap), the jet intensity between 
the pair of stations concerned. 

Fig. 3 shows the relation between Ap and 
height plotted on a logarithmic scale for station 
pairs in the Hess winter data. The upper 
curve gives the pressure difference between 
Charleston and Sault-Ste-Marie, latitudes 33° 
and 46%° respectively. The flat maximum 
occurs at about 7.3 km, the mean momentum 
jet height for this latitude range, and is prac- 
tically symmetrical for a range of s—6 km 
on either side of this central point. The profile 
is almost identical with one for Charleston- 
Caribou as yet unpublished, being closly 
represented by a profile of the form (r), as 
may be verified by plotting the cumulated 
value of Ap against z on probability paper. 

The middle curve represents the mean 
vertical profile of westerly momentum for the 
latitude range 46%°—58%° (Sault-Ste-Marie 
and Port Harrison). The mean momentum 
jet level is at 6.6 km. 

The lowest curve, Miami, latitude 26°, and 
Charleston, is again of the form (1). The mo- 
mentum jet level is at 8.4 km. The same type 
of profile is found for the pressure-difference 
between Port Harrison and Arctic Bay (lati- 
tude 73°), the momentum jet height being 
6.3 km. 

We see, therefore, that the shape of the 
vertical momentum profile is characteristic 
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over a latitude range, 25° N—73° N of eastern 
North America. It can be shown that the 
profile shape (1) holds for West and Central 
U.S., and in different seasons. It is thus a 
quasi-invariant of the mean westerly flow 
over the U.S. 

Since the momentum profile shows a very 
flat maximum the Clayton-Egnell rule, mo- 
mentum constant with height, gives a good 
approximation over quite a considerable 
range of height. For about 2 km on either 
side of the jet level the reduction in momen- 
tum is only about 6 %, while for a range of 8 
km centred at the jet level the maximum 
departure is 20%. 

A parameter specifying the properties of 
the vertical profile is the equivalent thickness, 


Zr, defined by 
ZT Apm = fe Ap dz 
re) 


Zr is the thickness of an equivalent stratum of 
air having everywhere a momentum equal 
to the momentum at jet level, with the same 
total momentum as the actual atmosphere. 
The average winter equivalent thickness is 
12.5 km between latitudes 26° and 33°, 
11.4 km between 33° and 46%°, 11.5 km 
between 46%° and 5814°, and 10.0 km for 
latitudes 5812° to 73°. These are actually un- 
derestimates, for the momentum below 1 km 
and above 16 km. has been disregarded. It 
can be shown that the equivalent thickness 
does not vary significantly with season; it is, 
in fact, a quasi-invariant of the mean circulation. 
The total westerly momentum between 
latitudes gr and g2 in a meridional slab of 
unit thickness extending vertically to infinity 1s 


oo 


P2 
af dp f oUdz 


fi 0 


When one of the latitudes lies in the region of 
low level easterlies, the lower limit of the z 
integral becomes the height of the doldrum 
axis in that latitude. 

For the total momentum between, say, lati- 
tudes 30° and 60° we may thus write 


60° 


I ZT 4) 
22 it sin @ (2 a ay 


30° 


me 
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dt 1 
where (2 is the pressure gradient along 
m 
the jet profile. Hence to an approximation, 
since the equivalent thickness varies only 


slightly with latitude, 


where @ is some mean latitude. 

The total mean westerly momentum 
between two latitudes is thus proportional 
to the pressure-difference between 30° and 
60° at the height Zr, or 7 km. 

The total momentum over 
westerlies belt is given by 


the whole 


U’ = na Zr Cot y (Ap)m/2 (3) 


where (Ap)m must now be interpreted at the 
zonal index at 7 km between latitudes 30° 
and 60°. 

The zonal index at other levels gives a 
(more approximate) measure of the total 
atmospheric momentum, because of the quasi- 
invariance of the vertical profile-shape. Howe- 
ver, the 7 km level is to be preferred as the 
momentum changes with height least here, 
and is less subject to random fluctuations. The 
worst level to take is the surface, as here the 
profile is subject to fluctuations connected 
with the boundary conditions which may not 
be appreciably reflected in the mean flow at 


height. 


4. The use of contours 


Hess’ data are given in terms of pressure 
at a fixed level. It is now standard practice 
in most weather services to quote contour 
heights of standard pressure surfaces rather 
than pressures at fixed height levels. Fig. 4 
shows the mean contour heights for the year 
1946—49, inclusive, as a function of cos y 
for U.S. stations around 80° W longitude for 
the 300 and 700 mb levels. 

These profiles show the same salient features 
as those drawn for pressure at fixed height 
levels. They have the advantage that the slope 
determines velocity directly, and not momen- 
tum. By reading off the slopes of the contour 
profiles for a number of different pressure- 
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Fig. 4. 700 mb contours (starred) and 300 mb con- 
tours (circles) as a function of latitude for meridian 80° W. 
approx, for 1946—49. The upper circles indicate 300 mb, 
contour height in the year of greatest height (1949) and 
the lower circles apply to 1947, the year of least profiles. 
The intermediate points, through which the smooth curve 

has been drawn, refer to the mean four year values. 


levels, it is possible to compute cross-sections 
with much greater facility and speed. 

Three points (circled) are shown for each 
station in the 300 mb profile in fig. 4. The 
middle point is the four year mean contour 
height, and through these points the smooth 
curve of the mean profile is drawn. The upper 
set of points gives the mean contour heights 
in the year of greatest general contour height 
(1949), while the lowest set corresponds to 
1947, the year of least general contour height. 
It will be seen that the annual mean height of 
the 300 mb surface may change by more 
than so metres from year to year, but that the 
surface rises or falls by much the same amount 
in all latitudes. If the fall or rise were uniform 
in all latitudes the slope relations would 
remain unchanged and the West-wind pattern 
would be unaltered from year to year. 

However, we should expect a greater varia- 
bility of mean contour height in high than in 
low latitudes. This is borne out in fig. 4, 
where it will be seen that the ‘49 and 47 mean 
contours differ by 30 m at Miami and by 
80 m at Greensboro. We should consequently 
expect years of high average contour height 
to be weak jet years, and years of low average 
height to display strong jets. 

The main jet latitude is approximately 36° 
for both 1947 and 1949, but the westerly 
speed at the 300 mb level is 33 % higher in 
1947 than in 1949. There is a subsidiary maxi- 
mum between the latitudes of Pittsburg and 
Buffalo in 1949, and between Buffalo and 
Caribou, or in higher latitude, in 1947, this 
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subsidiary jet being stronger in 1949. Fig. 4 
thus reveals appreciable variations in the 
strength and structure of the annual mean 
westerlies from year to year, although the 
overall pattern is not radically altered. It 
should be borne in mind, also, that systematic 
errors, changing from year to year cannot be 
completely excluded from the estimates of 
mean 300 mb contour heights. 


5. The vertical wind profile 


The difference in the mean contour height 
of a given pressure-surface between two 
stations is a measure of the mean geostrophic 
wind component between the stations in a 
direction transverse to the line joining them. 
If one station is due North of the other the 
westerly component of mean wind is measured. 

The way in which this height-difference 
varies with the pressure-surface gives a meas- 
ure of the change in geostrophic wind-speed 
with height. In an application such as this it 
is natural to take pressure as the measure of 
height. Whether we take as our argument the 
actual pressure, or its logarithm, as being 
more nearly proportional to height is purely a 
matter of convenience. In fig. 5 we plot 
mean contour height difference between 
Charleston and Caribou against pressure for 
the four-year period, 1946—49. The curve 
shown is the reduced profile, normalized so 
that (12 

It will be noticed that the points lie evenly 
on a smooth curve which is almost a linear 
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Fig. 5. Mean vertical westerly wind profile between two 

latitudes. The full curve is drawn through the four year 

mean contour-differences, Charleston—Caribou. The 

circled points represent the same relationship for Pitts- 

burg—Caribou, and the starred points are for Greens- 

boro—Pittsburg. The triangles apply to Charleston— 
Miami. 
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function of pressure up to jet level. The 
linear relationship 


U = Uo (p— pm) | (po — Pm) + 
=P Uns: (po —p) / (po — Pm) (4) 


where po is the surface pressure and p,, the 
pressure at jet level, Uo and U,, being the 
corresponding westerly geostrophic speeds, 
is indicated by the dotted line in fig. 5. If a 
closer relationship is desired a parabolic curve 
could be fitted, but for many purposes the 
linear relation will be found quite accurate 
enough. The jet level is estimated to be 
220 mb. 

The profile shown is that of mean wester- 
lies over a considerable latitude range. It is 
therefore of interest to see whether the shape 
is substantially unchanged for a narrower 
range of latitude. The starred points in fig. 5 
give the vertical profile between Greensboro 
and Pittsburg, while the circled points show 
the relationship for Pittsburg-Caribou. The 
points lie very close to the curve drawn for 
Charleston-Caribou, so that the profile shape 
is very nearly independent of the latitude 
between 33° and 47°. Above 200 mb there is 
some considerable scatter indicating that ob- 
servations at the highest levels may not be 
absolutely representative. 

Even below latitude 33° the vertical pro- 
file shape is very similar to that for higher 
latitudes. The points enclosed in triangles give 
the profile-shape for Miami-Charleston. The 
scatter is more noticeable, probably due to the 
small difference in contour height between 
these stations, but the linear relationship is 
still a reasonable approximation. The jet 
level is somewhat higher than in higher lati- 
tudes, at 190 mb. 

It is also possible to show that the vertical 
profile-shape does not change significantly 
with season. In fig. 6 the mean reduced 
vertical profile, Charleston-Caribou, 1946—49 
is again drawn, together with points for the 
months January, April, July and October, 
1949, for the same pair of stations. The points 
for the first three months cluster very closely 
around the four-year mean profile, with no 
significant seasonal change. The most dis- 
cordant month, October 1949, is indicated 
by the lower line. The jet level is higher, 


and the curve more concave than for the 
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Fig. 6. The four year mean westerly wind profile for 

Charleston—Caribou — full curve, with indications of 

seasonal variation. The points refer to the different sea- 

sons, but it is apparent that the seasonal scatter of pro- 
file shape is slight. 


annual mean. Even so, the similarity between 
the two curves is striking. 

The scatter of monthly values from the 
four-year mean profile is most pronounced 
below 850 mb. This is perhaps not unexpected; 
the wind-field at low levels is most likely to 
be influenced by surface boundary conditions 
which may have nothing to do with the 
general westerly flow-invariance of the wind- 
speed profile. That it does not vary appreciably 
with longitude over the U. S. can be in- 
ferred from the invariance of its counterpart, 
the momentum profile. 


6. The annual variation of geostrophic wind 


We have remarked that the vertical profile 
of geostrophic mean wind is a quasi-invariant 
of the atmospheric circulation over the U.S.; 
it does not change appreciably with season. 
Hence the seasonal variation of wind at any 
one pressure-level may be taken as typical 


\ —+ + + + + 5 
fs J F M a M J J a 5 ° 
Trace 


o | 


Fig. 7. The seasonal variation of mean westerly wind 
between Charleston—Caribou. 
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of the variation at all levels, except, possibly, 
near the surface. 

In fig. 7 we show the mean annual varia- 
tion of contour-height difference between 
Charleston and Caribou, 1946—49 at the 
200 mb. level. 

The geostrophic Westerly speed of the 
mean circulation reaches a maximum in 
January and a minimum in July. The curve is 
approximately sinoidal in shape. 

The curve shown is based on four years of 
data; the periodic curve in any one year may 
depart markedly from this pattern. Thus the 
maximum in a given year may well occur in 
December or February, or even March; the 
minimum might be displaced to June or 
August. The strength of the mean monthly 
circulation in any month may be displaced a 
month or two from the seasonal mean. In 
the most variable month, the departure from 
the seasonal mean will be 25% or more 
above or below normal in 30% of months. 
Nevertheless, in spite of large ‘casual’ varia- 
tions in the strength of the circulation, there 
is a strong underlying seasonal rhythm. The 
October plateau, showing up in fig. 7, may 
be a real effect, or may be rte simply to the 
inclusion of a run of abnormally sluggish 
Octobers. Such points can only be settled 
by a longer series of data than that used here. 


7. The kinetic energy of middle latitudes 


The mean westerly wind-speed at a given 
pressure level between two stations on the 
same meridian separated by a latitude distance 


Ay is 
gAZ 


ric 2 0 sin 9 Ay 


where AZ is the contour-height difference, 
and 9 is a mean latitude. The total momentum 
over the latitude distance A y of a meridional 
slab of unit thickness is, then 


Po 


Ay f @Udz = ~~~ f AZ dp 
oO 


OCT 


Now we have seen that À Z as a function of 
pressure (height) is a quasi-invariant of the 
mean circulation 


A= (4 Z)m f(b), F (Pm) = 
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We may thus write for the total momentum 


of the slab 


(AZ) 7 
U =: arora ioe 


and for the total momentum of an entire 
latitude belt, assuming the functional form of 
f(p) to be the same in all longitudes, 


Po 
fa (AZ) cot p J fdp 
Q + 
We see, therefore, that the total westerly 
momentum in a given latitude range is pro- 
portional to the difference in the mean con- 
tour heights at the level of maximum geostro- 
phic wind-speed. The integral 


Po 


J fdp 


has the dimensions of pressure. We term it 
pr the pressure-thickness. It represents the 
weight, in pressure units of atmosphere, which, 
possessing everywhere the maximum wind- 
speed, has the same total momentum as the 
actual atmosphere. 

From the vertical profile of geostrophic 
speed, already determined, we find, 


pr = 540 mb. 


Since f is a quasi-invariant of the westerly 
circulation, so is pr. Hence a convenient 
overall index of the strength of the westerly 
circulation is the mean contour height differ- 
ence at a level of approximately 200 mb. 

This formulation is in every way equiva- 
lent to that in terms of the zonal index, for 
which we have the equivalent thickness as 
the quasi-invariant. Whether the pressure 
thickness or the equivalent thickness is used to 
characterize the mean vertical profile depends 
entirely on whether we are dealing with 
contours or fixed level charts. 

The total kinetic energy in a vertical longi- 
tudinal slab of unit thickness between two 
latitudes is 
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Hence 

PUS 
sdZm_ (5) 
We may thus define a kinetic pressure-thickness, 


Po 


Pr KR dp 


Since f is quasi-invariant, so is px. From the 
4-year mean vertical profile between Char- 
leston and Caribou, we find 


px = 385 mb. 


We see from the expression (s) that the 
total atmospheric geostrophic kinetic energy 
between two latitudes is proportional to the 
square of the contour height difference at 
about 200 mb between the two latitudes. 


8. Other longitudes 


It is important to know whether the rela- 
tionships found apply to the N. American 
continent only, or generally characterize the 
entire Westerly belt. Yeu (1950) has published 
a winter cross-section along 76° E, using actual 
rawin data rather than geostrophic speeds. 
The general structure is strikingly similar to 
that found by Hess for Eastern America. Two 
jets occur around latitudes 28° and 42°. Both 
these positions are further south than the 
corresponding jets over America, and jet 
speeds are about 20% higher. In latitude 40° 
the speed of jet level is some 6—8 times that 
at 900 mb, in agreement with what is found 
over N. America. From the appearance of 
the cross-section it would seem that wind-speed 
is a nearly linear function of decreasing pres- 
sure, as over N. America. À profile published 
by CHAUDHURY (1950) shows approximately 
the same characteristics. 

From this it might appear that the structure 
of the westerlies is approximately the same in 
all longitudes. At the moment there are no 
cross-sections available for Western Europe, 
a circumstance which might call for some 
surprise, considering the wealth of data 
available.‘ 


1 Since completing this paper the author has drawn 
a geostrophic cross-section down the Greenwich me- 
ridian (Met. Mag. Dec. 1951). This cross-section is in 
substantial agreement with the conclusions advanced 
here. 
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Indications point to the structure of the 
mean circulation near the zero meridian differ- 
ing appreciably from that found over N. 
America, or Asia. Fig. 8 shows the mean 
contour difference (normalized) between Lark- 
hill (latitude 52°) and Lerwick (60°) for the 
two year period 1948—49. While in a general 
way the vertical profile for the meridian of 
Greenwich resembles that found over N. 
America in corresponding latitudes, there are 
perceptible differences. The relation between 
geostrophic speed and decreasing pressure 
does not show such a close approach to linear- 
ity. Moreover, the jet height is closer to the 
300 mb surface than to the 200 mb surface 
The jet speed for Larkhill-Lerwick is only 
about three times the 950 mb speed, the 
corresponding ratio over similar latitudes ni 
N. America being 6 or more. 


Perhaps the most striking feature of the 
westerly circulation near the Greenwich meri- 
dian is its great variability, as constrasted with 
the relative steadiness of the Aow over N. 
America. Monthly mean normalized contour 
height differences are plotted in fig. 8. They 
show a very wide scatter about the long-term 
average profile shape, with no significant 
seasonal variation. In fig. 6 we indicated the 
monthly ’scatter of profile shape over the 
eastern U. S., and comparing the two figures 
it will be at once apparent that the scatter is 
altogether greater in the case of the British 
figures. The two dotted lines in fig. 8 
indicate the envelopes ‚of extreme scatter 
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Fig. 8. The mean annual vertical wind profile for the 

Greenwich meridian (Larkhill—Lerwick). The points 

indicate the monthly scatter, the dotted lines, the extre- 
mes within which the profile-shape may lie. 
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Fig. 9. S. Hemispheric mean vertical west wind profiles. 
Circles summer, stars winter. Upper set of curves, Mac- 
quarie Is.-Laverton, lower set Auckland—Taieri. 


about the mean vertical profile of geostrophic 
wind. It is evident that in any ‘extreme’ 
month the profile shape may bear very little 
resemblance to its long-term mean. We may 
contrast the profile steadiness over N. America 
and over U.K. by saying that for 4 months of 
the year the mean monthly wind at, say, 
soo mb will depart by more than 6% % of 
the value predicted from the long-term mean 
profile-shape over N. America, but will de- 
part by 25% or more from the long-term 
value in the region of the U.K. In the U. S. it 
is possible to take the long-term mean profile 
as a good approximation to the actual monthly 
profile, but over U.K. this is no longer 
possible; the profile-shape is not quasi-inva- 
riant. 

These facts demand further factual investiga- 
tion. Why is the stability of mean westerly 
flow so strikingly different in different regions 
of the earth? 

RIEHL (1949) has conjectured that the ob- 
served stability in the structure of the Wester- 
lies is due to lee effects of high geographical 
barriers. The westerly flow is notably uni- 
form in geographical barriers. The westerly 
flow is notably uniform in lee of the Andes 
(Ru) and of the Himalayas (Jen). The quasi- 
invariance of the westerly flow-pattern over 
castern N. America may similarly be attri- 
buted to the effects of the Rockies. This ex- 
planation would attribute the instability of 
the flow near the Greenwich meridian to the 
absence of a high barrier to westward, with a 
resulting absence of a quasi-stationary lee flow- 
pattern. 

If this explanation were the correct one 
vertical profiles in the Southern Hemisphere, 
except in lee of the Andes, may be expected to 
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display a structure similar to that found for 
the British Stations. In fig. 9 four southern 
hemisphere profiles are plotted. The top pair 
relate to summer and winter for the station- 
pair, Laverton (38° S, 145° E) and Macquarie 
Island (54°S, 159°E) (Lozwe and Rapox, 
1950). The summer profile is indicated by the 
circled points. These profiles are more re- 
miniscent of the U. K. than the U.S. profiles. 
The jet height is nearer the 300 mb level in 
summer, and the 200 mb level in winter. 
The ratio of jet speed to 950 mb speed is 
over 3 in summer and under 3 in winter, 
in agreement with the value of this ratio ob- 
served for Larkhill-Lerwick. One feature not 
observed in Northern Hemisphere data is 
that the geostrophic speed at jet level is higher 
in summer than in winter. This feature is 
reproduced also in the profiles for Auckland 
(37° S, 175° E) and. Taten (46° See 
(Hurcuincs 1950). The jet height for these 
two stations is closer to the 200 mb level 
both in summer and in winter. It will be borne 
in mind that the Southern profiles are based 
on very limited data — only one month of 
Macquarie Island data, and two months 
from Taieri. However, the profiles shows 
the same characteristics, although the data are 
limited, and refer to different times. There is a 
strong presumption for supposing the mean 
southern profiles to resemble closely those 
found in corresponding regions of the Northern 
Hemisphere. 

On the Riehl hypothesis one would attri- 
bute the N. American type of vertical wind 
profile to the presence of an orographic 
barrier, and the U.K. type of profile to the 
absence of any barrier to westward. The 
facts presented here are certainly not incon- 
sistent with such a presumption. 

On the other hand, an alternative conjec- 
ture equally fits the facts. It might be postu- 
lated that the pattern of flow observed over 
N. America is typical of a continental region, 
while the U.K. type of profile is equally 
typical of an oceanic area. 

Over a continent there is a considerable 
surface frictional drag, with a resultant con- 
version of westerly angular-momentum. I 
this loss of A. M. is to be made good by high- 
level advection, it is reasonable to anticipate 
strong mean westerlies aloft, and a marked 
gradient of wind with height. 
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On the other hand, over an oceanic region 
the A. M. depletion at the surface is less for a 
given pressure gradient, and a smaller turn- 
round of westerly momentum at height will be 
required to make good the losses. Hence a 
smaller increase in westerly wind with height 
will suffice to maintain a steady state. 

Of great interest is the apparent unsteadiness 
of the oceanic circulation, as manifested by 
the variability of the monthly mean profiles 
near the Greenwich meridian. A more inten- 
sive study of oceanic profiles is necessary if 
we are to establish the nature of this instability, 
and its implications for the global balance. 

The unsteadiness of the mean profile over 
Oceanic regions may be a reflection of the 
actual variability of wind, associated with the 
frequency of migrating disturbances. BROOKS 
and CO-WORKERS (1950) have published world 
maps of standard vector wind deviation for 


207 


levels between 700 mb and 130 mb. In these 
charts the global pattern is clear, a maximum 
variability over the oceans in middle latitudes, 
and minima over the continents. 

Minima are found, at all levels throughout 
the year, to be associated with the Himalayas, 
Andes and Rockies. The African minimum 
shows a decided preference for the highest 
ground from Tanganyika to Abyssinia. A 
trough connected with the European Alps is 
also in evidence up to 500 mb. 

This might be taken as evidence favouring 
Riehl’s hypothesis, although it may be difficult 
to show why the presence of perturbations 
in the mean flow should affect the stability of 
the flow itself. Nevertheless, there seems to 
be a distinct association between wind variabil- 
ity and the year to year variability in mean 
flow. 
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Abstract 


From solutions of the linearized vorticity equation applied to systems with a variable basic 
zonal current, it is shown that the effect of damped disturbances is to feed energy into the zonal 
flow and transfer angular momentum from zones of easterlies into zones of westerlies, thereby 
strengthening the zonal circulation of the atmosphere. Amplifying disturbances withdraw 
energy, from the zonal flow and transfer angular momentum from zones of westerlies into 
zones of easterlies, thereby weakening the strong zonal circulation and transforming it into cellu- 
lar vortex motions. Thus the migratory cyclones and anticyclones are necessary elements for 
developing and maintaining the zonal circulation, while, on the other hand, they also act 
occasionally as brakes on the general circulation, keeping it from growing indefinitely. 

The neutral disturbance corresponding to a typical zonal wind profile in winter is investigated 
by numerical integration. The stream line pattern is found to be rather similar to a typical how 
pattern in the upper troposphere during a period of little change, indicating that this flow 


pattern is associated with the nearly neutral or slightly damped oscillations. 


I. Introduction 


When the motions of the atmosphere have 
been smoothed by taking averages along the 
latitude circles and over longer periods to 
eliminate the transient irregularities, we obtain 
systems of zonal currents, the so called pre- 
vailing westerlies and casterlies, which repre- 
sent the more prominent overall pattern of 
the general circulation. This zonal wind system 
is interrupted from time to time by the pres- 
ence of disturbances, both in the westerlies 
and easterlies, the number of which is much 
greater at lower levels than at upper levels. 
One of the fundamental problems that has 
always confronted meteorologists is the ex- 
planation of the nature of the physical pro- 
cesses that take place in the general circulation 
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and also in the secondary circulations associated 
with it. It is evident that a complete answer 
to this problem can be achieved only through 
many different approaches which take into 
account physical processes in the atmosphere, 
among which those maintaining the global 
balance of energy and of momentum for the 
zonal wind systems are fundamental. 

The concept of the balance of angular 
momentum in the general circulation was 
used first by JErrreys (1926), who pointed out 
that the maintenance of the prevailing winds 
in a steady state against surface friction necessi- 
tates a transfer of westerly angular momen- 
tum from the zone of surface easterlies into 
the zone of westerlies. However, the prevailing 
zonal winds are by no means steady. They go 
through cycles of decay and rebuilding, in 
addition to the mean seasonal changes in 
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intensity and location. Statistical investiga- 
tions have shown that the general circulation 
in middle latitudes tends to oscillate between 
extremes of strong and weak zonal motions, 
the time interval between consecutive maxima 
being of the order of magnitude of 4—6 
weeks (ALLEN, 1940; WILLETT, 1947) Thus 
the study of the general circulation should not 
be restricted to the problem of maintaining 
the steady zonal winds against ground fric- 
tion, but should also include the mechanisms 
for their growth and decay. However, it is 
most probable that the process that maintains 
the zonal winds is also responsible for their 
growth and decay and therefore the two 
problems can be studied simultaneously. But 
in view of the complexities of the natural 
process, we cannot expect to obtain a complete 
solution of the problem from any single set 
of equations; the physical processes them- 
selves probably work through many different 
stages, so that a reasonable solution can only 
be achieved through many different steps. 
The purpose of this paper is to investigate 
some essential mechanisms pertaining to the 
overall atmospheric circulations, that are 
directly connected with the development of 
the zonal winds. 

In the past, the mechanism of the general 
circulation which maintains the energy and 
momentum balances, has by and large been 
looked upon as consisting of a large-scale 
convective process, with ascending motions 
in lower latitudes and descending motions in 
higher latitudes, the development of a large 
zonal velocity being attributed to the action 
of the deflecting force of the earth’s rotation. 
Recently the importance of the mean vertical 
motion for the maintenance of the zonal wind 
systems has been questioned by STARR (1948 a, 
1949), who suggests that the horizontal energy 
and momentum transfer can be accounted 
for by the horizontal components of the gen- 
eral circulation and of the lesser circulations, 
or by the eddies with quasi-vertical axes. 
This idea has also been suggested by Rossby 
in his theory of horizontal mixing (1947, 
1949). The investigations of WIDGER (1949), 
Mintz (1949) and Wuite (1950) have sub- 
stantially confirmed this point of view; there- 
fore the possibility of the existence of mean 
meridional circulations will be disregarded in 


this study. 
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In the atmosphere the zonal wind systems 
are most clearly defined above about $ kms. 
Above this level, especially around the level 
of the tropopause where the maximum of the 
meridional momentum transfer occurs, the 
transfer of heat (WHITE, 1950) is relatively of 
little importance and the motion of the air is 
controlled largely by the equation of vorticity 
about the vertical and the continuity equation. 
Below this level, the heat supply and heat 
transfer must play an important role and the 
instability of this layer, which can be inferred 
from the presence of the large number of 
disturbances there, is probably associated with 
these factors, although the exact nature of 
this instability is still not very clear to us. 
This type of instability has been discussed by 
CHARNEY (1947) Eapy (1949) Jaw (1946) and 
Berson (1949). On the other hand, in the 
upper troposphere and lower stratosphere, the 
atmosphere generally moves with a persistent 
flow pattern, indicating that it is in a relatively 
stable state, while very unstable conditions 
happen only occasionally (for example, the 
formation of cut-off lows and highs). Since in 
this study we are mainly concerned with the 
processes which govern the maintenance and 
development of the zonal wind system, it is 
sufficient if we restrict the investigation to 
those processes which are prominent in the 
upper part of the troposphere, assuming the 
existence and not concerning ourselves with 
the origin of disturbances in this field. We 
shall thus divide the atmosphere roughly into 
two layers and study the motions within the 
upper layer only. However, this partition into 
two layers is taken only as a convenient means 
of analysis rather than a real barrier for transfer 
from one region to the other. Energy can be 
received from the lower layer through irreg- 
ular disturbances, and can also be given by 
the upper to the lower layer on some occasions. 
Thus there is a mechanism in replenishing the 
kinetic energy of the upper system, even 
though we treat it as barotropic. In this way, 
we consider the lower layer merely as a 
source or sink of disturbances for the upper 
layer, and the divergence and vertical motion 
as a means of communication between these 
two layers. 

In order to separate the divergence effect 
from the overall processes, we shall assume 
that it takes place intermittently within some 
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specific time intervals, so that after some 
disturbances have been produced in the upper 
layer, the subsequent motions are horizontal 
nondivergent and therefore governed by the 
simple vorticity equation without divergence. 
In this way, the most important problems 
concerning the dynamics of the upper at- 
mosphere can be discussed by the solutions 
of the simple vorticity equation. It can be 
shown by integrating the equations with 
respect to the mass distribution along the 
vertical that the mean motion of the whole 
atmosphere can also be approximately describ- 
ed by this two-dimensional barotropic model. 
Therefore the results obtained may also be 
considered to represent the mean motion. 

Since there is no other energy source in this 
simplified model, the total kinetic energy 
must remain constant during the time under 
consideration. Thus any change of the kinetic 
energy of the mean zonal flow must come 
from the kinetic energy of the disturbances, 
and the study amounts to an investigation of 
the interactions between the lesser circulations 
and the zonal wind systems. In the following, 
the method of harmonic perturbations (both 
in x and in time f) will be adopted. The prob- 
lem is therefore reduced to the discussion of 
the stability of the zonal wind profiles with 
respect to the disturbances. Here the term 
“Stability” is defined according to the reac- 
tions of the disturbances to the zonal wind 
systems in the horizontal plane. In this respect 
three kinds of disturbances can be distinguished: 
the disturbances which withdraw kinetic 
energy from the mean zonal flow are said to 
be amplified, those which supply kinetic energy 
to the mean zonal flow are said to be damped 
and those which neither give nor take energy 
from the mean zonal flow are neutral. Since 
this study is not aimed at explaining the origin 
of the disturbances (although it may help to 
explain the development of some large disturb- 
ances), these disturbances are taken to have 
been added to the upper layer from below. 
It then becomes physically feasible to discuss 
the different effects produced by damped as 
well as amplified disturbances. We shall de- 
monstrate later that one of the effects is to 
transfer westerly angular momentum from the 
belts of slower zonal flow into belts of stronger 
zonal flow, therefore to strengthen the pre- 
vailing zonal winds and to overcome the fric- 
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tional effects, although the disturbances them- 
selves carry no resultant momentum. On the 
other hand, the amplified disturbances take 
kinetic energy from the main flow and transfer 
zonal momentum away from the zone of 
strong westerlies. Thus it appears that the zonal 
mean motion is not merely a statistical result 
of these irregular disturbances, as is frequently 
stated in descriptive studies of the general 
circulation. Rather the general circulation is 
built up and regulated by the disturbances 
(migratory cyclones and anticyclones). Its 
growth and decay are closely connected with 
the activity and development of the disturb- 
ances and these in return are connected with 
the stability of the zonal flow. 


2. The rates of change of the mean zonal 
velocity and kinetic energy 


The equations of motion over a spherical 
earth for horizontal and frictionless flow may 
be written 


Ou Ay A op I neg 
ot o a cos DA 2 a cos Doi 
(ra) 
dv rae I Ov}? 
ot Né 0 28 2 28 (1b) 
where Z =f +Ë = 20 sin ® + 
I dv du cos PT . 
PRET [5 — | is the vertical 


component of the absolute vorticity, hereafter 
referred to simply as the absolute vorticity, & 
the relative vorticity, u and v the zonal and 
meridional components of the horizontal veloc- 
ity, Un? = u? + v?, @ the density, p the pres- 
sure, a the mean radius of the earth, and A 
and ® are the longitude and latitude res- 
pectively. 

As only large-scale motions connected with 
the general circulation will be discussed, all 
quantities may be averaged along latitude 
circles. Thus the velocity components and 
vorticity can be written as u = % +4’, v = 
D+v', and Z = Z +’. Here the bars denote 
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average values defined as i = (27) / udA, and 


(6) 
u', v’, © are the departures from the average 
values, which can also be taken as representing 
a perturbation. The motion is assumed to be 
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purely horizontal and nondivergent within 
the time interval under consideration. This 


2% 
implies / vd} = 2 xd =o. In this respect the 
O 


motion is similar to a geostrophic velocity 
field. Since the effect of the variation of the 
density o in x-direction is small, o will be 
taken as constant unless otherwise stated. The 
above condition then means that there is no 
net mass transport across the latitude circle at 
that level. If the motion at this level is repre- 
sentative of the mean atmospheric flow, this 
is equivalent to the condition 
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When equation (ra) is integrated along the 
latitude circle and use is made of the non- 
divergent assumption, we find that the rate 
of change of the mean zonal velocity is given by 
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= 2 za: Ou cos? ® represents the total north- 
ward transfer of angular momentum across the 
latitude circle. When the convergence of 
meridians is neglected, linear momentum can 
be used and t= — ouv. Thus the time rate of 
change of the mean zonal momentum in the 
narrow zone contained between the latitudes 
® and ® + 6@ is proportional to the net in- 
ward transport of angular momentum into 
this zone, namely Jt /a cos? BA, or alterna- 
tively to the northward transport of vorticity. 

In discussing the development of the general 
circulation it is useful to consider also the time 
rate of change of the kinetic energy of the 
mean zonal flow. This may be dealt with 
quite generally as follows. Multiplying the 
general momentum equation in the zonal 
direction (obtained by combining the first 
equation of motion and the continuity equa- 
tion, both in three-dimensional forms) by the 
mean zonal velocity # and integrating over 
the whole volume » of air contained between 


STI 


the two latitude surfaces ® = Os, and D = M,, 


we get 
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dz 
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Here it is assumed that the variation of o with 
longitude and with time are both small and 
have been neglected. Equation (3) shows that 
the kinetic energy of the mean zonal flow can 
change by a conversion from the kinetic 
energy of the eddy motions, an inflow of 
this energy through the boundaries, and 
possibly by the presence of a mean meridional 
cell. The pressure force does not occur in this 
equation, therefore this energy cannot be 
produced directly by the pressure force. 

The significance of this equation for the 
study of the maintenance and development 
of the mean zonal winds has been discussed 
by the writer elsewhere (Kuo, 1951). It has 
been shown that only the horizontal eddy 
stress u’v’ and the mean meridional motion can 
serve as sources of kinetic energy of the mean 
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in the last member of this equation is of the 
right order of magnitude and sign to account 
for the total frictional loss. We shall there- 
fore restrict our investigations merely to this 
effect, which can be discussed adequately by 
using our simplified model. 

Since the production of the total kinetic 
energy is also of paramount importance in 
understanding the physical processes, we may 
also write down the equation for the time 


zonal flow, and that the term ou’ v’ cos D 
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rate of change of the total kinetic energy of 
the horizontal motion: 


2 2 I Te 
eee = —V. | Pm Se = onary | “Es 


> d 
op, = nz 6 nt) (4) 


where V- denotes the horizontal divergence. 
The significances of the different terms have 
been discussed by Starr (1948 b). It may be 
noted that the pressure force is certainly a 
source for the total kinetic energy although 
not for the mean zonal kinetic energy, hence 
it may be concluded that its effect is to pro- 
duce eddy and meridional motions. 


Suppose we now integrate equation (4) 
over the whole region & bounded by two 
latitude surfaces along which v = 0, the earth’s 
surface where w =o and the top of the at- 
mosphere where 9v,? = o and the motion 
is assumed to be nondivergent, then all the 
terms on the right disappear. The equation 
then expresses the constancy of the total 
kinetic energy for the entire region. Suppose 
now the complete motion is divided into a 

2n 
mean zonal flow u = (2x) * / udAand a per- 
Oo 
turbation w, v’ and denote the kinetic energy 
of the horizontal perturbation in Z by £” then 


we find by making use of (2) and (4) that 
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where the integration along z has been omitted 


and where À denotes the initial value of 
n/a cos y. The same equation can be obtained 
from (3). It should be pointed out that equation 
(s) is applicable only when the entire region 
is considered. For any particular region or a 
small belt, the transport of kinetic energy 
through the boundaries must be taken into 
consideration and therefore we must use 


| equation (4). 
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3. The perturbation equations and 
boundary conditions 


In the upper troposphere and lower stra- 
tosphere, the large-scale motions are to a large 
extent horizontal, nondivergent and baro- 
tropic, we therefore assume that this motion 
may be described by the vorticity equation 
and the equation of continuity. Because of 
the difficulty of finding general solution of 
non-linear partial differential equation, the 
vorticity equation will be linearised and it 
then becomes 


d 
dt 


where to is the basic current depending upon 
® only, Zo is the absolute vorticity of the 
basic current, » is the kinematic viscosity. 
The suffix zero occuring in Ao and zo denotes 
initial values of À = u/a cos ® and Z. This 
equation will be used to investigate the 
changes from the initial state. The frictional 
term on the right is generally small and can be 
neglected except as discussed below. Since the 
motion is assumed nondivergent, the conti- 
nuity equation takes the form 


‘ee: du’ d(v'cos P)] _ 
VE a cos D [5 * I® ] a 
(7) 


and a stream function Y for the perturbation 
velocities may be introduced, defined by 
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where the dots indicate time derivatives. In 
terms of this stream function, the relative 
vorticity ¢’ of the perturbation is given by 


0? y , d? l'A 
0% | dp 
== log (tan ® + sec D). 
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where y 


The perturbation stream function can be 
represented by the sum of a number of elemen- 
tary independent harmonics of the form 


YY — Y (y) ein (A—ut) (9) 
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where n is the wave number around the globe 
and u may be complex, u = u,+ iu;; My re- 
presents the angular phase velocity and ny; 
the amplifying factor of the particular har- 
monic disturbance. When (9) is substituted 
in the linearised vorticity equation (6) we find 


(lo — u) (W — n? 2) — dy y = 


= — irn [Pr — 2m" + ns] (ro) 
where the primes denote differentiation with 
respect to y. This equation must be satisfied 
by all the component harmonic perturbations. 
Since » is extremely small in the atmosphere, 
the terms on the right-hand side of (10) can 
generally be neglected. In this case we get 
the nonviscous differential equation 


(11) 


However, it should be pointed out that if the 
frictional term is neglected entirely in the 
vorticity equation (6), then no distinction can 
be drawn between amplified and damped 
disturbances, since either the harmonic factor 
einQ—Ht) or e-in@—at) will give the same dif- 
ferential equation (11). Therefore we must 
take (11) as the limit of (10) with » extremely 
small, with the necessary corrections for its 
solutions in the vicinity of the singular points 
where u,= Ao. 

The general solution of (11) for each har- 
monic contains two arbitrary constants. No 
specific information can be obtained without 
imposing some restrictions on their values. 
The necessary restrictions are the boundary 
conditions that must be satisfied by each 
harmonic. In the atmosphere, the only natural 
horizontal boundary conditions are those 
needed to make the perturbations finite at 
every point, therefore X and d¥/dy must be 
equal to zero at the two poles. However, 
since the motion over the southern hemisphere 
is generally unknown to us, our discussion is 
necessarily restricted to the northern hemi- 
sphere. It therefore appears best to introduce 
the equator as a line of symmetry, and disregard 
all the seasonal variations. Thus we may use 
the following different set of boundary con- 


ditions under different circumstances: (a) ¥=0 
at the pole, which must be satisfied, and 
either d¥//dy = 0 (symmetric) or Y= 0 (anti- 
symmetric) at the equator. Each of these condi- 
tions will give a zero momentum transfer 
(t = 0) at the equator. In addition to the study 
of the hemisphere as a whole, we may consider 
a perturbation confined to a smaller belt. 
For this case, we use (a) Y = 0 along two lati- 
tude circles. Although these boundaries are 
rather arbitrary, it might be of interest for 
example to examine the belt within the 
subtropics which separates the easterlies from 
the westerlies. It may be noted that 7, defined 
in (2), vanishes at the boundaries. 

By applying one boundary condition to 
the general solution of the particular har- 
monic, one of the two arbitrary constants can 
be eliminated, while the other boundary condi- 
tion enables us to determine the permissible 
value of u for any given n or conversely. 
Therefore the total perturbation stream func- 
tion is given by 

Ot A7) = Lgl; (Yen 4" (128) 
where the arbitrary constants c; are deter- 
mined by the initial disturbances 
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Thus, the harmonics that are present at any 
moment correspond to those present at the 
initial moment and the solution is unique. 
When every component disturbance is known, 
then (12a) can be used to forecast the velocity 
field. But at the present, we can only discuss 
the general properties of the individual har- 
monic components as related to the profile 
of the basic zonal currents. 


4. The solutions of the vorticity equation 
and the nature of the disturbances 


The solution of equation (11) has been dis- 
cussed by the author in a previous study (Kuo, 
1949) for the motion over an infinite plane or a 
cylindrical earth, in which the coriolis para- 
meter was taken as a linear function of ® and 
the profile of the basic current symmetrical. 
Though the results obtained there are still 
generally valid for the motion over spherical 
earth and for asymmetrical profiles, certain 
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modifications are unavoidable and conse- 
quently we shall now discuss the problem 
further. 

Equation (11) shows that if the phase ve- 
locity y lies within the range of values of Ao, 
the point y=y, where u = Ao is a singular point 
of the equation.! The solutions show certain 
peculiarities around this point and require 
careful examination. These questions do not 
exist for those disturbances whose phase 
velocity is distinct from the velocity of the 
basic current. It is therefore necessary to divide 
the disturbances according to their velocities of 
propagation into three groups: (a) those 
propagating faster than the maximum velocity 


max. 4o of the basic current, (b) those travelling 


slower than the minimum basic current min. Ao, 
and (c) those propagating with a velocity 
within these two limits. 

The disturbances in the groups (a) and (b) 


will be discussed very briefly. It can be shown, 


= h—dho/dy> (where 


dy 
h=2(w + Ao) sec h? y + 2 tanh yddo/dy corres- 
ponding to f for motion in a plane) and employ- 
ing Sturm’s oscillation theorem, that under the 
conditions existing in the atmosphere (h > o 
generally) there can be no neutral disturbance 
in (a). That is to say, in general, no neutral 
disturbance that is governed by the simple 
vorticity equation will travel faster than the 
maximum velocity of the basic flow in middle 
latitude westerlies and equatorial easterlies. 
On the other hand, there always exist neutral 
disturbances in group (b), since h(i — Ao)~* is 
then negative (see Kuo, 1949). 


by writing 


Now let us examine the possibility of having 
amplified or damped disturbances in these two 
groups. If u is complex, Y must also be com- 
plex, Y=, +i¥;. Thus, when equation 
(11) is multiplied by the conjugate function 
w* — Y,— iV, we get by separating the real 
and imaginary parts: 


d d¥, 
dy EE dy 


+4 u nz 
lee 


(13a) 


1 Suffix S will always refer to the value of a function 
at the singular point y =y,. 
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dope dpe oad se NRE 
dZo 
. (m 7) VAE 13b) 
a 


where W is the Wronskian of X, and W,, 


49, Ww, a; 


W = ÿ; dy dy > and 2 — (Ao => 


— pop" Ba = ts [do = lg 
Equation (13b), when integrated from one 
boundary to the other, determines the phase 
velocity. When (13a) is integrated from one 


boundary to the other we find 


W: 


her 


which shows that à necessary condition for 


dZo 


having a complex ju, or g: + 0, is to have 7 


change its sign somewhere within the region, 
since the boundary conditions require both Wı 
and W;to be zero. This condition holds both 
for the amplified and damped disturbances in 
the groups (a) and (b), since there is no 
singular point for these perturbations. 

In the atmosphere, the absolute vorticity Zo 
is largely determined by the coriolis para- 
meter and therefore increases with latitude 
over most of the hemisphere. From (13 b) it 
can be seen that under this condition there is 
no disturbance at all in group (a) and only 
neutral disturbances in group (b). Rossby’s 
long waves are examples of these neutral 
disturbances. 

Group (c). This group includes nearly all 
the wave disturbances actually observed in the 
atmosphere, which travel with phase velocities 
within the range of velocities of the basic 
current, min. Ao < 4, < max. Ao. For these 
disturbances, the points y, where #— À are 
singular points of (11). (For the jet stream 
profiles, there are generally two singular 
points for every value of u, one on each side 
of the jet). Only one solution of the equation 
(11) is regular at this point, while the second 
involves a term with the factor log (y — y.) 
(Kuo, 1949), and therefore is many-valued. 
Since the motion cannot be many-valued, the 
change of the value of y in passing from the 
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negative to the positive side of y, must be 
determined properly, by finding the path of 
continuation around y,. This is obtained by 
considering the solutions of (11) as the limits 
of the asymptotic solutions of (10) with a 
very small coefficient of viscosity. Is has been 
shown by many investigators (e.g., LIN, 1945; 
Wasow, 1947; KUO 1949) that this conti- 
nuation must be taken along a path below the 
singular point y(sı) where dAo/dy is positive 
and above the singular point y(s:) where 
dAo/dy is negative. Therefore,” the solutions 
should be taken as 


Wi=nt+a2n?+ ou 
Pa=1+ban+ ….. 
W=1+by?4 


(15a) 


for 7-0 


+G, Wr log 7 
...+G Pı (log |n|—ai) for A 


(1sb) 
where 


TA = GAS 
n= y—y,and G, = ( ) ( ) 2 42. 
LEE dy : dy k 


However, according to Lin (1945) and Wasow 
(1947), the set of solutions are not valid within 
a small region around the singular point for 
the damped disturbances. This means that 
although the neutral and amplified disturb- 
ances can be discussed by employing the 
inviscid vorticity equation (11), some other 
effects need be considered within the small 
region around the singular points in discussing 
the damped disturbances. This point becomes 
clear in the case when the absolute vorticity 
Zo is a monotone increasing function of y: 

With this in mind, we can use the solutions 
(15). According to the two expressions of Y2 
on the two sides of y,, there is a change of 


the value of = from »—e to y. +e, 


given by 


(7) - (7) =. 

dy Jy=y,+e dy /y=y—e 

He ASS 
‚Ge 69 


dy 
There is, therefore, a rapid change in u 


ern 


aE MAGE 
across a singular point if ns is not zero 


and & very small. 


Now we may proceed to discuss the sta- 
bility of the zonal winds for the disturbances 
in group (c). Since the inviscid vorticity 
equation (11) 1s valid for all real y when yu; > 0, a 
necessary condition for the existence of ampli- 
fying disturbances, as can be seen from (14), 
is that dZo/dy should change its sign somewhere 
in the region, the same as for the disturbances 
in the groups (a) and (b). Unlike the disturb- 
ances in the groups (a) and (b), this condition 
is not necessary for damped disturbances, 
because equation (11) and its solution are not 
valid for some parts of the real y axis. 

The same condition is also necessary for the 
existence of neutral disturbances in this group. 
This has been proved in the previous study 
by the writer (Kuo, 1949) by considering 
the jump of the Wronskian across the singular 
point, given by 


Wy, Fe) — W(y,—@) = 


dZo\ | dis | * 

=y ( ) Bs 2? (17) 
Since for neutral disturbances both Y, and Y; 
must be solutions of equation (11), the varia- 
tion of W is zero at all points except the singular 
points where there are jumps given by (17). 
As W is zero at the two boundaries, the sum 
of these jumps must be zero. For the symmetric 


jet profiles the values of the two jumps are 
the same; therefore each must be zero. Thus 


either (2) or Y, must be zero. The impossi- 
bility of the latter has been proved by the 
writer (Kuo, 1949) by the method of com- 
parison, which can also be used for other 
profiles. It can also be proved by employing 
the expression (13b). Thus dZo/dy must change 
its sign and go to zero in the region. Also 
dZofdy must be zero at the singular point 
for symmetrical jets; that is to say, this neu- 
tral disturbance travels with the velocity of 


the basic current Ao(c) at the critical point? (the 


point where dig = 0). This is also true for the 


dy 


unsymmetrical jets when the value of Ao is 
the same at the two critical points. 


2 Suffix c will always refer to the value of a function 
at the critical point y = y(c). 
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For unsymmetrical jets in which the values 


of À at the two critical points ye) and y(c2) 
are not the same, the condition that the total 


i fW ish requi 422 and 
jump o must vanish requires ion 


dZo pe À 

iy to have opposite signs, which assumes 
$2 

the vanishing of dZo/dy somewhere in the 

region. This neutral disturbance travels with a 


phase velocity between the two values of Ao 
at the two critical points, as has been proved 
by Foote and Lin (1950). 


It has thus been shown that a necessary 
condition for the existence of amplifying 
disturbances is that dZo/dy should equal zero 
somewhere in the region and change its sign. 
The same condition is also necessary for the 
existence of neutral disturbances in group (c) 
and of damped disturbances in the two groups 
(a) and (b), but it is not necessary for damped 
disturbances in group (c). 


5. Equations of angular momentum transfer, 
energy exchange and the stability of zonal 
flow 


The relation between the meridional trans- 
port of angular momentum and the stability 
of the zonal flow has been discussed by the 
writer in a previous investigation (Kuo, 
1949). While the results obtained there are 
correct for the particular profile studied, they 
need some modifications and _ elaborations. 
The same problem has also been discussed by 
Macuta (1949) and HASSANEIN (1949). These 
authors show that even with a constant basic 
current, some perturbations are able to pro- 
duce a constant meridional transfer of momen- 
tum. However, these perturbations do not 
satisfy the boundary conditions of this paper, 
and furthermore, a constant momentum 
transfer without accumulation does not help 
the momentum budget, as can be seen from 
equation (2). Because of the importance of 
this momentum transfer and the energy 
exchange between the basic zonal current 
and the perturbations, they will be discussed 
in some detail in the following two sections. 

The perturbation velocity components u’ 
and v’ can be expressed in terms of the stream 
function, 
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W = — É£ cos nË — = sın né | e"Hit 
dy dy 


and v = — n[#, sin n£ + W; cos né] ei, 
where £=A—~y,t is the phase angle. From 
these expressions and by making use of 
equation (13a), the average northward an- 
gular momentum transfer across unit length 
of the latitude circle is found to be given by 


ou'v’ cos? D = — T= 


I ; |7 |2 dZo 
= — onu;e?"Hi ; 
= A en 
yr 


dy (18) 


This transfer can also be expressed in terms 
of the tilting of the troughs and ridges, defined 
as the curves along which the meridional 
velocity is zero. The location of the trough 
and ridge curves is determined by tan n &= 
— W;/, and therefore the angle « which they 


make with the northerly direction is given by 


d(né à 
2 tan &« = Yr sec? nË | hes Ms 


ady 
y 
— | 4 | aZi 
— nf Er Fr d) (19) 
Yı 


Since | ¥ i is proportional to 92, we see that 
the angular momentum transfer given by 
equation (18) is proportional to the product 
of the mean kinetic energy of the perturbation 
and the tilting of the troughs and ridges. 

The vorticity transport can either be ob- 
tained in the same way, or simply by differ- 
net (18) with respect to y. Thus we 
lave 


|v 


i LS ad d 2nu 
v er er 
ki dZo 
- — - y 20 
| À “4 ll |: dy | | ( ) 


This equation shows that the vorticity transfer 
produced by the damped disturbances is in 
the direction of increasing Zo while that pro- 
duced by the amplifying disturbances is in the 
direction of decreasing Zo, therefore it is 
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not always in the direction of the absolute 
vorticity gradient. When the zonal flow is 
stable, it actually is transported in the direc- 
tion opposite to the gradient of vorticity, 
therefore has the effect of increasing this 
gradient; only a very unstable flow will re- 
sult in a thorough mixing and a reduction 
of the existing vorticity gradient. 

Substituting equation (18) into (2), we get 
the time rate of change of the mean zonal 
velocity: 


Ri 
dy Ve 


2 
2 


ou _ fui F 47e 
dt 2 


anust 


(1) 


It also follows that the effect of the damped 
disturbances is to produce an increase of u in 
the region where dZo/dy is positive and a 
decrease where dZo/dy is negative, thus their 
effect is to intensify the zonal flow. The 
effect of the amplifying disturbances is in the 
opposite direction, that is to flatten the jet- 
stream of westerlies and easterlies. 

When the expression (18) for r is substituted 
nto the energy equation (5), it becomes 


2; 
OE’ dr 
Da | à Fo = 
Di ; 


Da 
a : QUE |2 adZo 
= —onuje*™it | Ao — 
2 © : i | Ao a. u | dd 


1 


d® 


(22) 


From (22) it may be seen that if the pertur- 
bation kinetic energy is to increase (u; > o), 
then, besides the condititon that dZo/dy should 
change its sign somewhere in the region, we 
must also have positive dZo/dy associated with 


higher values of À and negative dZo/dy with 


lower values of Ao. This may be considered 
as a second necessary condition for the occur- 
rence of amplified disturbances and is equi- 
valent to the condition given by FJÖRTOFT 
(1950). For the profile of a jet-stream of 
westerlies, this condition is automatically satis- 
fied when dZo/dy changes its sign within the 
region. It is not necessarily satisfied however, 
for example, an easterly current around the 
equator can satisfy these requirements simul- 
taneously only if the current is sufficiently 
strong. 
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It can now be shown that for the jet stream 
profiles, these two conditions are sufficient for 
the occurrence of amplified disturbances. 
That is to say, when they are fulfilled, equation 
(tr) and the given boundary conditions ac- 
tually permit solutions with a positive w;. 

The existence of amplified disturbances for a 
symmetric jet stream velocity profile has been 
proved in a previous study (Kuo, 1949) by 
applying the Sturm oscillation theorem. This 
is done by first establishing the existence of a 
neutral disturbance with u,= Ao(c). For the 
jet-type velocity profile in which the current 
velocities at the two critical points cr and c2 
are the same, and dZo/dy is positive in the 
central part and negative in the outer parts of 
the region, the ratio [A—A(c)|-'dZo/dy is 
always positive and finite. It follows that a 
neutral disturbance with = Ao(c) exists, since 
the solutions of the equation 


dZo 
aye; Ald | i dy | 
dy? wen iE en) AR (c) 


become oscillatory for sufficiently small values 
of n2. When we investigate the eigen value 
u =u+iu; in the neighborhood of this 
eigen-value 4 =Ao(c) we find that (Kuo, 1949) 
if u, > Ao(c) then u; < 0, and if u, < Ao(c), then 
u; > 0. Hence the disturbances which are 
travelling with a phase velocity greater than 


the critical values Ao(c) will be damped and 
those which are travelling with a phase velo- 


city less than the critical value Ao(c) will be 
amplified. 

The existence of neutral and amplified 
disturbances for the profiles with Ao(cx) and 
Ne (c2) differing slightly from each other can be 
discussed by the argument used by Foote 
and Lin (1950). From the above discussions 
one can summarize the necessary and sufficient 
conditions for the existence of amplified and 


neutral disturbances (with w,=— +) as follows: 


If dZoldy equals zero and changes its sign 
dZo 
somewhere in the region and if ue > o where 


ho is large and "< 0 where do is, small, then 
amplified and neutral disturbances may exist. 
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A necessary and sufficient condition for the 

À : AZo: 

existence of damped disturbances is that y's 


not identically zero in the entire region. 

In applying the equations (18) — (22) to 
the damped disturbances propagating with the 
current velocity of some point, it must be 
remembered however that the solution does 
not hold in the small region around the 
singular point. Since this interval is generally 
small, we can neglect the large fluctuations 
within it and just compute the change of 7 
across this interval from (17). Then equation 
(22) can also be used to discuss the damping 
effect. The dissipation of kinetic energy 
through friction within the narrow viscous 
strip can only increase this damping effect, 
and not reverse it. 


6. Momentum and vorticity transfer and 
energy exchange for different velocity 
profiles 


We shall next consider the meridional mo- 
mentum and vorticity transfer and energy 
exchange between the perturbations and the 
mean flow for some typical profiles. 


(a) Jet-stream profile in which the values of 2. 


at the two critical points are the same. 
Assume that there are two critical points 


0 


(poins where = 0), y(cx) and y(c2), one 
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Fig. 1. (a) The mean zonal current and the mean absolute 
vorticity distriburions. The dashed curve represents the 
mean zonal winds produced by damped disturbances at 
some later time. (b) Schematic orientation of the troughs 
and ridges (T) and the meridional distribution of the 
stress T= — ou’v’ cos?® for damped disturbances; the 
arrows indicate direction of vorticity transport. (c) 
Same as (b) for amplified disturbances. 
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to each side of the jet stream of westerlies. The 
profile must be such that Zo decreases with 
latitude in the outer parts of the region and 
increases in the central part, as is shown in 
fig. 1. For this profile, a neutral disturbance 


exists which propagates with the velocity Ao(¢), 
which has its troughs directed from north to 
south, which gives zero momentum transfer 
at all points, and which causes no energy 
exchange between the basic flow and the 
perturbation. 


Now we discuss the damped disturbance 
which moves with a velocity only slightly 


greater than the critical velocity Ao(c), and 
having singular points. The eddy stress 1 
(southward angular momentum transport) 
produced by this damped disturbance (u; < 0) 
decreases slowly from zero at y= y: to a 
negative value up to the point y(sr) — &, 
according to equation (18); then it decreases 
rapidly across the singular point, correspond- 
ing to the change of W given by (17), up to 
y(sr) + €, since dZo/dy is positive at y(sı). 
Although the inviscid theory fails to be valid 
within the small belt 2 ¢ around y(sı) we may 
approximate to the distribution of 7 withi 

this small interval by the dotted curve, as is 
indicated in fig. tb. As the change of + across 
this interval is known, its values at all other 
points where the inviscid theory holds can be 
computed from (18). Thus it increases from a 
negative value at y(s:) + € to a positive value 
at y(s2) — e, where it has another rapid decrease 
up to y(s2)+ &, and then decreases to zero as 
the other boundary is approached. The vortic- 
ity transfer is northward in the central region 
between y(cı) and y(c2) and southward in the 
outer regions, as denoted by the small arrows. 
The direction of the trough and ridge lines 
for these damped disturbances is represented 
by the curve T. It is seen that the troughs and 
ridges are oriented with their convex-side 
toward the cast in the vicinity of the jet 
stream, and momentum is transferred from the 
outer parts into the central part, (from regions 
where the zonal current is small into the re- 
gions where it is great), therefore having the 
effect of intensifying the zonal currents. It is 


: A 
also seen that t and the wind shear — have 


dy 


Opposite signs in most parts of the region; 
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therefore the damped disturbances supply 
their energy into the main zonal flow. The 
dashed curve in fig. 1a represents the zonal 
wind distribution produced by the damped 
disturbances at some later time. 

For the amplifying disturbances, 7 is given 
by (18) at every point. It increases from zero 
at y=yı to a positive value, has a higher 
rate of increase near the singular point, and 
attains its maximum value at y(c:), then 
decreases to a negative minimum value at 
y(c2), and increases to zero as the other bound- 
ary is approached, as is shown in fig. Ic. 
Thus positive t is associated mainly with 
positive dAo/dy and negative t with negative 
dio/dy, showing that the kinetic energy of 
the perturbation will increase with time. The 
vorticity transfer is in the opposite direction 
as produced by the damped disturbances. 
The trough lines and ridge lines are oriented 
with their convex side toward the west. 
There is a flow of angular momentum from 
the central part into the outer parts, with the 
effect of flattening the profile of the zonal 
motion. 

For the disturbances without singular points 


(u, + ne the general features of the distribu- 
tion of t and the orientation of the trough 
and ridge-lines are the same as have been 
discussed above. 

The results obtained above agree with 
those obtained in a previous study (Kuo, 1949) 
for symmetric profiles and over a plane or a 
cylindrical earth. When the asymmetry is not 
great, the general conclusions obtained also 
hold for the profile with different values of Ao 
at the two critical points. For this kind of 
profile, a constant momentum transport can be 
produced between the two critical points 
by a neutral disturbance as has been discussed 
by Foore and Lin (1950). However, no 
other effect on the mean zonal flow will be 
produced by this neutral disturbance except a 


gradual reduction of this difference of Ao 
between the two critical points. 


(b) Jet is limited to middle latitudes while the 
vorticity distribution at higher and lower latit- 
udes is dominated by the coriolis parameter. 


In the atmosphere, even when the westerlies 
are strong, and the absolute vorticity has a 
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maximum to the north and a minimum to 
the south of the jet, the negative gradient of 
vorticity generally does not extend to the 
equator, but Zo attains a second weaker maxi- 
mum at a lower latitude, showing the predom- 
inance of the effect of earth’s rotation. This 
type of vorticity distribution, such as is re- 
presented in fig. 2, occurs much oftener than 
the case (a) and therefore deserves separate 
discussion. 


In this case, there and four critical points in 
the zonal wind profile where dZo/dy vanishes, 
the two inner ones having higher values of 


Ao. If Ao is the same at these two inner critical 
points, y(c:) and y(c3), a neutral disturbance 
will travel with this velocity. Another neutral 
disturbance corresponding to the two outer 
critical points may also exist. 


For damped disturbances, the distribution 
of t and the direction of the troughs and 
ridges are represented by the curves in fig. 
2b. It starts from zero and increases to a maxi- 
mum at y(cr), then decreases from y(cı) to 
y(c2). Assuming that this damped disturbance 


is travelling with a velocity greater than Ao(c2), 


so that (=) is positive, then 7 will have a 
dy /sı 
sudden decrease across this singular point, 
and reach a large negative value at y(sr) +e. It 
increases from there to the point y(s2) —¢, 
where it reaches a positive value, and then 
decreases rapidly across the singular point 
y(s2). It attains a small negative value as y(c4) is 
approached and increases to zero near the 


dZo . 
pole. It is seen that 7 —— is mainly negative, 


dy 
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Fig. 2. Same as Fig. 1 for a weaker zonal wind profile. 
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therefore the total kinetic energy of the 
perturbation decreases with time, being fed 
into the zonal flow, at least in part. Also, angu- 
lar momentum is transported into the central 
part of the zone of westerlies, strengthening 
the westerlies and also developing an easterly 
current north of the equator, thus gradually 
increasing the degree of instability. For this 
case, the troughs and ridges are tilted with 
their convex side toward the east in the zone 
of westerlies, but the direction is reversed in 
the equatorial belt and also near the pole. 


For the amplifying disturbances, the distri- 
bution of r is represented by the curve in 


fig. ac. It is seen that 7 dAo/dy is mainly positive 
and has large values in the central part of the 
profile, so that the perturbation is withdrawing 
kinetic energy from the main flow. But in 


the outer parts, 1 dAo/dy is negative, so that 
the perturbation is also supplying kinetic 
energy into the main flow, strengthening the 
zonal flow in those parts, as can be seen from 
equation (22). Thus this case is relatively 
stable as compared with the case (a). 


It should be pointed out that these dis- 
cussions of the distribution of 7 with respect 
to the zonal wind profile apply only to the 
individual situations. Because the jet-stream 
of westerlies fluctuate both in their latitudinal 
position and strength, and because different 
effects are produced by different kinds of 
disturbances, the distribution of the average 
values of t for a longer period may differ 
considerably from that given by the corre- 
sponding mean profile. 


-20 m sec"! -10 


o b c 


Fig. 3. Distribution of t and vorticity transfer in the 
equatorial easterlies. (b) for amplified disturbances. 
(c) for damped disturbances. 
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(c) Waves in the high evel easterlies in 
equatorial belt 


Although there is no natural boundary 
between the high level westerlies in middle 
and higher latitudes and the high level easter- 
lies in equatorial latitudes, different kinds of 
perturbations exist in these two regions. The 
disturbances in the westerlies move rather 
rapidly towards the cast and are numerous, 
while in the easterlies there are relatively few 
large disturbances, and these generally move 
slowly towards the west or remain stationary. 
Therefore we may take the equatorial easter- 
lies as a separate region. 


The case when the easterly current is more 
concentrated in the equatorial belt and flat 
in the outer regions, with dZo/dy negative 
around the equator and positive in the outer 
parts will be discussed. According to the dis- 
cussion in section $, both damped and ampli- 
fied disturbances can exist in this case. An 
examination of equation (18) shows that for 
amplified disturbances, 7 is positive north of 
the equator and negative to the south, with 
the same sign as the shear of the basic current, 
as indicated in fig. 3b. This gives an increase 
of the perturbation kinetic energy. The 
trough and ridge lines for these amplifying 
disturbances are oriented with their convex 
side toward the east, contrary to that in the 
zone of westerlies. This tilting gives an equator- 
ward flow of westerly momentum and will 
tend to flatten and possibly destroy the zonal 
flow pattern. 


For damped disturbances with phase velocity 


within the range of À in the belt, the distribu- 
tion of 7 is represented by the curve in fig. 3c, 
with the proper rapid change across the 
singular point. It is seen that 7 is negative to 
the north and positive to the south of the 
equator, so that the perturbation kinetic 
energy will decrease. The troughs and ridges 
are tilted with their convex side toward the 
west, and momentum is transferred awav 
from the equator, with the result that the 
easterly current is strengthened near the 
equator, increasing the degree of instability. 


Besides these disturbances with u, within 


the range of Ao in the belt, we might also have 
disturbances moving toward the west at a 
speed faster than the maximum easterly current. 
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The general features of the distribution of r are 
the same as those described above. 


(d) When dZo/dy is positive in the entire re- 
gion 

When the prevailing zonal winds are weak, 
the value of Zo is determined primarily by the 
coriolis parameter, therefore may increase 
with latitude over the entire hemisphere. 
According to the preceding discussions, this 
case is stable and cannot be disturbed spon- 
taneously. The same conclusion may also be 
reached by applying the non-linearized vor- 
ticity equation dZo/dt =o (FjÖRTOFT, 1950; 
STARR, 1950; KUO, 1950). However, the 
linear theory also tells us that all possible 
neutral disturbances travel with phase velocity 


u < min. Ao. Thus if any disturbance is added to 
the system and is moving with a velocity 


within the range of Ao, which is generally the 
case because added vorticity moves with the 
current, the motion cannot follow the law of 
conservation of absolute vorticity at every 
point, for some other effects necessarily enter 
into the problem (either viscosity or diver- 
gence), at least in some small regions. If we 
allow viscosity into the small region around 
the singular points, then damping effect will 
take place. 


Although viscosity is important for these 
damped disturbances, its effect will be small 
outside the small strips around the singular 
points. Therefore the distribution of the stress T 
is still given by equation (18) in these regions, 
and increases with latitude. Across the small 
strips around the singular points, t has rapid 
decreases, attains a negative value at y(sz) +e, 
and goes to zero as the boundary y= Ya is 
approached. Corresponding to these negative 
jumps of r, the mean zonal current decreases 
within these strips. The mean vorticity 
decreases around y(sr) and increases around 
y(s2), indicating a tendency of creating a mini- 
mum and a maximum of Z along these lati- 
tudes respectively. 


One physical explanation of this process is 
that high concentrations of vorticity are 
always present in the small-scale disturbances, 
which under the influence of the northward 
increase of Z, will result in a transport of 
cyclonic vorticity toward higher latitudes 


5—104318 
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(Kuo, 1950). The accumulative action of 
eddy viscosity near the singular points in 
spreading out the existing vorticity concentra- 
tions helps to build up extreme values of 
vorticity along these latitudes. 


7. The neutral disturbances in the upper 
troposphere in winter 


In order to compare the general results 
obtained above with actual flow patterns in 
the atmosphere, a solution of the perturbation 
equation for some normal zonal wind profile 
is desirable. Here the study will be restricted 
to the neutral disturbances which probably 
represent the flow pattern observed on the 
upper air weather charts during a period of 
relatively little change. 

The average basic current studied in this 
section is given in fig. 4. It is taken to represent 
a mean profile of the layer from 900 mb to 
100 mb, and resembles a typical profile at the 
500 mb level across a section of strong wester- 
lies. In order to simplify the integration, the 
profile has been greatly smoothed and adjusted 
to make Ao have the same value at the two 
inner critical points. 

For this neutral disturbance, the phase 


velocity is given by the current velocity at the 


inner critical points. The corresponding wave 
number and the eigen function can be ob- 
tained by Taylor’s expansion. Thus, if Y and 
its derivatives at the starting point y= yrare 
known, the value of Yand its first derivative 
at the nearby point y = y: + 7 are given by 


(ys + 9) = F (ys) + ny) + 


where the primes denote differentiation with 
respect to y and 


OW (yz) = Wie (ys) = WY (yr — 9) 
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Fig. 4. The mean zonal wind profile and absolute vor- 


ticity distribution in winter, with wpwy = 40 m. sect, 


Since Y has to satisfy the vorticity equation at 
every point, the second and higher derivatives 
have to be computed from (11) by successive 
differentiation and by use of the values of the 
lower derivatives already obtained. Thus 


pr = [en] 7 


= (m — F),, (y) 
Yu = (nt — F)  — Pp 
Div = m — F)W"— 2 Fw — FW 


Il 


and so on, and Y can be integrated in a step-by- 
step manner. 

To start the integration, the solutions corre- 
sponding to the distribution of F near the 
equator are obtained. In our case, F is equal to 
— 23 near the equator. Thus, for y = 0, (11) be- 
comes Y" = mW, where m? = n? + 23, and the 
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solution that satisfies the boundary condition 
P(o) =ois 
= A sinh my (24) 
The first and higher derivatives for the starting 
point can be computed from this solution and 
the values of and its derivatives at y =0 + 
computed from the Taylor expansions. The 
value of à Wiv at the starting point can either 
be neglected or computed from the solution 
given above. If the condition u’= 0 = —9W/dy 
is used for y = 0, we get another solution. Since 
the difference is small, it will not be discussed. 
Since the other boundary condition applies 
at y = © (the pole) in our coordinate system, 
it is necessary to replace it by the condition 


P'+nY =0 (25) 


for a very large value of y, i.e., X should de- 
crease exponentially as the pole is approached. 
For the profile we are studying, the function F 
can be represented by F =—2 (y—o.8)-? = 
= — 2x7? for y > 1.3 (north of 65°N), to a 
high degree of approximation. Thus, for y> 1.3 
equation (11) becomes 


Jaa (m +2 Vo 


i (26) 


The solution of (26) which satisfies the condi- 
tion (25) is given by 


P=C (" oa =) ey (27) 


and is valid for y > 1.3. 


When the function Ÿ with its derivatives 
obtained from (24) has been integrated to y= 
= 1.3, the ratio Y/Y should agree with that 
given by (27). The value of n that fulfils this 
requirement is the eigenvalue of the problem, 
and the corresponding ¥ the eigenfunction. 
It is found that the neutral oscillation corre- 
sponding to the two inner critical points of 
the profile in fig. 4 has the eigenvalue n = 9.5. 
The phase velocity of this neutral disturbance 
is 4, = 0.05 ©, about 12 deg. lat per day. The 
soo mb and 300 mb maps during the winter 
months January and February show that the 
wave number is generally from 8 to to, and 
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the phase velocity is about 10—15 degrees of 
longitude per day. The eigenfunction Y and 
its first derivative are given in fig. 5. It is 
seen that the perturbations are much stronger 
on the north side of the jet than on the south 
side; this agrees with the statistical result 
obtained by RIEHL, YEH and La SEUR (1950) 
that the maximum value of 72 occurs generally 
to the north of the jet. 

In the perturbation stream function, the 
arbitrary constant has been chosen to give a 
maximum perturbation velocity u’ = 12 m 
sect while the maximum velocity of the 
basic current is 40 m sec=!. The perturbation 
kinetic energy is less than 10 per cent of that 
of the basic flow. 

Since only an integral wave number is 
permissible, we shall take = 9 instead of 9.5. 
The total stream function is plotted in fig. 6. 
It is seen that closed cyclonic and anticyclonic 
centers are produced to the north and south 
of the jet stream of westerlies as a result of the 


-500 


Fig. 5. The eigen function W and its derivative for the 
neutral disturbance. 


MH nt 


Fig. 6. Composite contour lines in meters. 
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perturbation. However, the anticyclonic cen- 
ters are very weak, much weaker than the 
cyclones to the north. In order to produce 
the large semi-stationary anticyclonic centers in 
the subtropics, it is necessary to superpose 
another slowly moving disturbance, largely 
determined by the two outer critical points 
in the profile. Since this generally does not 
effect the flow pattern in middle and higher 
latitudes, it will not be discussed. 


8. Summary 


From the discussions above, it is seen that 
different disturbances have quite different 
effects on the zonal flow, according to 
whether the basic flow is stable or unstable for 
the particular disturbance. In the upper 
atmosphere under normal conditions, the abso- 
lute vorticity is largely determined by the 
coriolis parameter, increasing toward the 
north in most regions. This situation is stable 
for most. disturbances. Thus, when disturb- 


. ances are added to the upper layer, most of 


them will undergo a damping process, trans- 
port zonal momentum from regions of 
weaker to regions of stronger zonal motion, 
and supply kinetic energy to the zonal flow. 
Parts of these zonal momentum and energy 
being transported downwards to overcome 
ground friction. The remainder will have the 
effect of strengthening and sharpening the 
zonal westerlies and easterlies, and thereby 
also increase the degree of instability. 

When the process of building up of the 
zonal wind systems through the damping 
disturbances has gone on for a sufficiently 
long period of time, the currents become 
strong and sharp, and maxima and minima of 
absolute vorticity are produced, the zonal 
wind distribution becomes unstable for the 
more effective disturbances, with longer wave 
lengths and smaller phase velocities. (From the 
vorticity distribution in fig. 5 and the solution 
obtained in the previous section, it may be 
inferred that the wave number of maximum 
instability is about 4 or 5 during the periods 
when the jet stream is strong in the winter 
months.) The amplitude of these disturbances 
will then increase to form closed cold lows to 
the south and warm highs to the north of the 
jet. As these amplified disturbances have the 
effect of transferring angular momentum from 
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regions of higher to regions of lower zonal 
velocity, and withdraw kinetic energy from 
the zonal flow, the intensity of the mean 
zonal flow is greatly reduced by these disturb- 
ances, and the zonal wind profile is flattened. 

After this breakdown, the absolute vorticity 
distribution is again dominated by the coriolis 
parameter, so that the situation is stable for 
other disturbances, and the damping effect 
operates once again, to rebuild the prevailing 
zonal winds. Thus according to this study, 
the general circulation will go through cycles 
of decay and rebuilding, due to the effect of 
these disturbances. It may be noted that 
their kinetic energy, at least in part, must be 
destroyed through the convergence effect, or 
transferred to the lower layer, so that the 
energy of the upper system will not continually 
increase. 


HSIAO-LAN KUO 


From the foregoing discussions, it may be 
concluded that the damped disturbances are really 
necessary parts in building up the prevailing zonal 
winds and overcoming the loss of kinetic energy 
and angular momentum through friction in the 
lower atmosphere. On the other hand, although 
large amplified disturbances occur only occasion- 
ally, they act in such a way as to prevent the pre- 
vailing winds from growing indefinitely, and 
therefore act as brakes for the general circulation. 
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Abstract 


Using a simple atmospheric model, the difference between the solutions of the primitive 
and the quasi-hydrostatic equations is demonstrated. The results show that the quasi-hydrostatic 
processes are largely determined by the initial field of pressure wherever the meteorological 
*‘noise’’=phenomena (in the form of short-period oscillations about a trend of the meteorological 

elements) depend on the initial field of the non-geostrophic windcomponents. This trend 
can be computed with sufficient accurracy also by a numerical integration of the primitive 
equations, which are mathematically simpler to handle than the quasi-hydrostatic equations. 

In this case the initial wind fields can be computed from the observed pressure field by using 


the geostrophic approximation. 


To find the trend the results have to be smoothed. 


Die hydro- und thermodynamischen Diffe- 
rentialgleichungen erfassen bei ihrer Anwen- 
dung auf atmosphärische Vorgänge nicht nur 
die relativ langsam ablaufenden sogenannten 
quasistatischen Prozesse, sondern zwangsläufig 
auch gewisse kurzperiodischen Abläufe, bei 
denen Ausbreitungsgeschwindigkeiten in der 
Größenordnung der Schallgeschwindigkeit auf- 
treten und die man deswegen als meteorolo- 
gischen Lärm bezeichnet. 

Im allgemeinen wirken sich diese kurz- 
periodischen Vorgänge mit relativ geringen 
Amplituden nur in Form einer Unruhe im 
zeitlichen Gang der meteorologischen Zu- 
standsvariabeln aus und werden in den synop- 
tischen Karten weitgehend weggeglättet, so 
dass man sich bei einer rechnerischen Vorher- 
sage mit der Erfassung der quasistatisch ver- 
laufenden Anderungsprozesse begnügen könnte. 
Beim numerischen Integrationsverfahren wird 
aber gerade der zeitliche Gang, also der par- 
tielle zeitliche Differentialquotient, dessen Be- 
rechnung auch ohne vollständige Integration 


aus den Differentialgleichungen verhältnis- 
mässig einfach ist, prognostisch zur Ermittlung 
der im nächsten Zeitintervall eintretenden 
Anderung der betrachteten Zustandsgrösse 
ausgewertet. 

Benutzt man für den numerischen Rechen- 
prozess die hydro- und thermodynamischen 
Grundgleichungen in unveränderter Form, so 
wird man bei Berechung der Tendenzen 
ungeglättete Momentanwerte erhalten. Aus 
den entsprechenden Differenzenquotienten las- 
sen sich demzufolge nur Änderungen für die 
betreffende Variable für ein Zeitintervall ent- 
nehmen, das klein ist gegenüber der Schwin- 
gungsdauer der kurzperiodischen Lärmvor- 
gänge, was einen Aufwand an Rechenarbeit 
erfordert, der eine rationelle Auswertung des 
numerischen Integrationsverfahren erschwert. 

Die Anzahl der erforderlichen Zeitschritte 
kann wesentlich reduziert werden, wenn diese 
für die Rechnung unbequemen, meteorolo- 
gisch uninteressanten kurzdauernden Vor- 
gänge aus den Differentialgleichungen elimi- 
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niert werden, was allerdings auch eine grund- 
legende Änderung in der Gestalt der Differen- 
tialgleichungen zur Folge hat. 

Im folgenden soll an Hand eines einfachen 
atmosphärischen Modells versucht werden, in 
den Lösungsfunktionen die kurzperiodischen 
Anteile von den quasistatischen zu separieren 
und abzuschätzen, für welche durch Druck- 
und Windfelder definierten Anfangszustände 
eine Elimination der meteorologischen Lärm- 
vorgänge berechtigt erscheint. 


I. Ausgangsgleichungen für das zugrunde- 
liegende atmosphärische Modell 


Als Ausgangsgleichungen sollen die be- 
reits von J. G. CHARNEY (1949) diskutierten 
linearisierten Gleichungen einer homogenen 
Atmosphäre zugrundegelegt werden. Bei Gül- 
tigkeit der hydrostatischen Grundgleichung 
und bei Annahme eines durch 


U=— © H, = const 
V=0 


definierten Grundzustandes schreiben sich diese: 


Du = — gh, + lv (1) 
IDPS — lu (2) 

Bewegungsgleichungen 
Dh = = v — Hu, (3) 


Kontinuitätsgleichung, 


wenn We = vz = Wy = Vy = h, = 0 vorausge- 
setzt, die geringe Abhängigkeit der Parameter 
l, H von y also als unerheblich betrachtet 
wird. In den Gleichungen bedeuten H Höhe 
der homogenen Atmosphäre, U, V Geschwin- 
digkeitskomponenten des Grundzustandes; 
h, u, v die entsprechenden Grössen der Stö- 
à 7 
dt ah ox”? 
Das Gleichungssystem lässt sich in dieser 
einfachen Form allgemein integrieren. Wegen 
der vorgenommenen Vereinfachungen werden 
zwar die gefundenen Lösungen tatsächliche 
Wettervorgänge nur sehr weitgehend idealisiert 
wiedergeben. Zur Demonstration der charak- 
teristischen Unterschiede in den Dimensionen 


rungsfelder; D = = 2 © sing. 
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und in der Abhängigkeit von den Anfangs- 
feldern zwischen kurz- und langperiodischen 
Prozessen dürfte jedoch das Modell ausreichend 
sein, so daß dieses zumindest einen Einblick 
in die Wirkungsweise des atmosphärischen 
Mechanismus vermittelt und einen Hinweis 
geben dürfte, auf welche Weise und unter 
welchen Voraussetzungen eine derartige Eli- 
mination auch in den nichtlinearen dreidi- 
mensionalen Diff. Gl. System möglich und 
berechtigt ist. 

Das System (1—3) läßt sich auf eine Diff. 
Gl. der Variablen h reduzieren. Aus (1, 2) er- 
hält man 


(D? + P)u = — gDh, 
0 à 
bzw. wegen D Er = D 
(D? + lux = —gDhyx (4) 
(Di + Bo = gi (s) 


Führt man diese Gl. in (3) ein, so resultiert 
eine homogene lineare Diff. Gl. III. Ordnung 
in h allein: 


d 


D (D:— gH =) h + Ph, =0; "10 


die im Falle ! > 0 in die Diff. Gl. fiir Schwere- 


= 


wellen (D: — çH =) h = o übergeht. Das 


Ox? 


System (4—6) ist dem System (1—3) völlig 
äquivalent und soll allgemein integriert wer- 
den, wobei insbesondere der EinfluB einer 
Variation der beliebig vorgebbaren Wind- 
felder u und v auf die zeitliche Entwicklung 
des durch die Größe h charakterisierten Druck- 
feldes untersucht werden soll. 

Wir beschränken uns dabei auf das Anfangs- 
wertproblem, bei dem zu einem bestimmten 
Zeitpunkt t= o die Funktionen h= h°, u = wo, 
v=v° für alle x-Werte vorgegeben sind, 
wodurch h (x, t), u (x, f) und v (x, £) eindeutig. ' 
bestimmt sind. 


II. Berechnung der Druckfunktion h (x, t) 
Der Ansatz h = hoAei &—® führt nach Ein- 

setzen in Gl. (6) auf die Frequenzengleichung 
— r (U— 03 + gHy?(U — à) — Pc=o 


AVY 


lle A 


1 ove 


ie ee LO dre. 


— tt NT 


Ré: ds), wee AL 
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UP 


bzw. auf A3—m2A + 
yp 


Ant 
wenn zur Abkürzung gesetzt wird 

m? = gH + = und | 

A =U—c (7a) 


Für U? < gH ~ 105 [m? sec”2] existieren drei 
reelle Wurzeln (Ar, Az, As), wovon Ax 


l , : 
; Az, Az nahe bei + m liegen. 


m2y2 


Aı läßt sich als Grenzwert eines Ketten- 
bruches darstellen mit dem Bildungsgesetz 


nahe bei 


UP 
y2 UP 


2 Ar ie) 
2 À 2 ? ; 
m I,n 


Aynntı = EN 
wre my 


Der Grenzwert Ar= lim A;, werde mit d 
n — 00 
abgekürzt. Die beiden weiteren Wurzeln A2, 


Az ergeben sich aus 


A+ +d + (d?—m?) = 0, 


Ax = ds um 
my? 
a=—1 + m3 -—+m; 
M=+ \ [ne — 3 
Be m2 Pret 
> 4 2 
(8) 
Cee U— NE, k =I, 2, 3. (9) 


Mit den drei Werten für c erhält man als all- 
gemeineren Ansatz für h 


3 
h(x, t) = ho > Aye? — 4) (ro) 

k=ı 
Der Lösungsansatz (10) erfaßt alle Lösungen mit 
Anfangsfeldern in h°, u°, v° in Form harmo- 
nischer Funktionen in x der gleichen Wellen- 
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2n 
y 

Anfangsfunktionen ergibt sich leicht durch 
Superposition von particulären Lösungen der 
Form (10) mit verschiedener Frequenz » zu 
Fourierreihen (bei periodischen Anfangsfel- 
dern). Man erhält bei periodischen Anfangs- 
funktionen für h°, u°, »° mit der Intervallänge 


länge L= . Eine Anpassung an beliebige 


1 Bus 27 
L* wegen L = —; y = — n für 
5 n Ex 
+ 00 3 
Mei) Ni hove), Let, 
n=—oo k=r 
L* 
Beer OT 1; 2 2. (11) 
2% 


worin noch die 44, aus den Anfangsbedin- 
gungen zu bestimmen wären. Da die Berech- 
nungsgleichungen für Ay sich für alle » 
entsprechen, genügt es, diese bei einer fest- 
gehaltenen Frequenz » aufzustellen. Z. Zt. 
t=o muß gelten: 


h(x, 0) = he 


und wegen (1—3) 
Barone: = Pe 2 


Dike RUE 


uo + H(gh2, — 1 vz) 
(12) 


Denkt man sich die Anfangsfunktionen für 
ho, u°, v° nach Fourier zerlegt, so möge für ein 
bestimmtes v gelten: 

h? = hove"* ho, Uo, Vo = komplexe Am- 
plituden der betrachteten 


Uy = Une“  Partialwelle z- Zt-t= o mit 


Vy = Ve”? 


L= oa (13) 
Setzt man (ro) und (13) in (12) ein, so erhält 


man zur Bestimmung der A4 folgendes lineares 
Gleichungssystem: 


ones I 
k=1 


t Zur vereinfachten Schreibweise soll künftig der 
Index » wieder fortgelassen werden. 
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z uo Ul Vo 

28 A Jel = F 

= ho gr ho 

= 

3 2 

ArAr = gH acl ae ee 

- gv? h v ho 
k=1 


(14) 


Mit den Bezeichnungen (8) für A; erhält man 
anstelle von (11) für eine bestimmte Frequenz 
(siehe Fussnote 1): 


hls) = hol Aveo) + [[A: + 
+ A3) cos v»Mt + (Az — À) i sin » Mt] - 


v[e—(u+5) | h (11a) 


und anstelle von (14) 


Ar+ (A, Ae A3) = I 


Ai (+4) + M(4;=4;)=F 


LIN 


2 


d'A: + (m: +) (A As) dM (4, A.) = C 


(14a) 


Für die Amplituden der Funktion h in Gl. 
(t1a) erhält man daraus: 


(m +2) —G—a(r +4) 
Ax = : = 
M2— 2 @ 
4 
A, + A; = I — Ax 
I d 3 


Zwecks bequemerer Variation des Anfangs- 
zustandes in u und v soll angesetzt werden 


’ Lo) ee 
9 = Uo EX = Oj ha ho e”* 
I 
' Odile 
vo = vo ee — (1 — e)i = ho er (16) 
I 


, gy ‘, 
Da vi =i ho e*, ug — 0 die Komponenten 
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des zu h° = ho e®* zugehörigen geostrophischen 
Windes bezeichnen, so stellen die dimen- 
sionslosen komplexen Zahlen 6 und & ein Maß 
dar für die Auslenkung des Anfangswindfeldes 
gegenüber dem geostrophischen. Die Beträge 
von 6, e geben das Verhältnis der Amplituden 
der nichtgeostrophischen Komponenten zur 
geostrophischen wieder, die Winkel von &, 6 
deren Phasenverschiebung gegenüber der geo- 
strophischen Windkomponente. 

Mit diesen Ansätzen (16) für u und v° 
erhält man anstelle von (14) 


3 

> A = ii 
k=1I : 

3 Öiv 

DE AA | — €) U — The ex | = Je 
k=1 

: oil 

Di ArA- LeHe == pe u|= G (14b) 
k=1 


Mit (14) und (15) ist die Funktion h(x,t) in 
(1x1 bzw. 11a) bestimmt. Bevor die auftretenden 
Amplituden A, numerisch berechnet werden, 
sollen die weiteren Funktionen u und v er- 
mittelt werden. 


Ill. Berechnung der horizontalen Wind- 
felder u und v 


Die Funktionen u, v gewinnt man nach 
Vorliegen der Lösung h aus Gln. (4 u. 5), 
in welchen die rechten Seiten als inhomogene 
Terme betrachtet werden können. Einsetzen 
von h(x, f) aus (10) in (4) und (5) ergibt: 


(D2 + Eu = — ghoi2v: - 


3 
oy ; 
. N ApAp RN und 


ET 


3% Li 
3 > Ay + ei (x — cpt) 
— 


(D? + b)v = ghoily 


k=1 
mit den Lösungen: 
Ar : 
D gho ae el? (~— cpt) + 
SAN iz 
k=ı 9? : 


+ (ou cos It + B,is in It) ei” x —Un 


7s. idées 


Le ed 
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el” A) + 


LRU 


+ (&, cos It + Byisin Ite” @-un 


Die Lösungen der homogenen Differential- 
gleichungen (4,5) verschwinden? jedoch, da 
= Pu=Xy= B,=0, wie man erkennt, wenn 
man die Lösungsfunktionen in Gln. (1—3) ein- 
setzt und die Frequenzengleichung (7) beach- 
tet, so daß man für u und v schließlich erhält: 


EIGENE IN er (e— cpt) — 


3 
v : 
+ en EN Bee 


k=I 


= + he {Biere ae [(B: un B;) cos vMt + 
al rot d 
+ i(B: — B;) sin mel (ei) (17) 


3 2 
v= a hoi > RB Ar er Eck) — 
l er De 


8 
Vv À 
= © hoi Y Cer) = 


k=I 


-7 huil C: ei” at) u [(C2 + C3) cos vMt + 
ed 
+ i(C2 — C3) sin »Mf] fs ( >’ ] (18) 


worin bedeutet: 


a l AkAk 
DR ae: 
k 3 k k D, 
i Ar 
(© iS Be 


Unter Benutzung der Frequenzengleichung in 


“ der Form: 


2 Im Falle » =o tritt allerdings die inertiale Fre- 
quenz | auf. Man erhält für die Anfangsfelder : 

h (x, 0) = u(x, 0) = 0; v(x, 0) = vo = const die 
Lösung: 


À UV: 
u = Vo sin It; v = vo cos lt; h = vo — sin It. 
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ee] 
3 n+2 
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> gH U\gHv x 


ergibt sich für 


y? DAT 
Ko= | GEN. es 
yp I Uo 
nee ace ai 
RE 


und damit für die dimensionslosen Koeflizien- 
ten: 


Br == 


B: u B; =di — By, C2 G— I—e—Cy 


BB, = = =. e+ re Sap, 


CCE ete u 


Eine Vorstellung über die Größenverhält- 
nisse der in den Funktionen h, u, v auftretenden 
Amplituden Ay, B,, C, in Abhängigkeit von 


der Wellenlänge L = == vermittelt die Ta- 


belle 1, in welcher 1 = 10-4 [sec], U=15 
[m - sec-1], gH=8 - 104 [m2 sec~?] angenom- 
men wurde. 


IV. Elimination des meteorologischen 
Lärmes 


Die Abschätzung der Größenverhältnisse 
ergibt, daB in dem meteorologisch interessie- 
renden Spektrum mit ausreichender Näherung 
gesetzt werden kann: 
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19 950°0 + 3 zooo‘o + zoo00'o — 


19 ıLz‘o + 8 €£00'0 + LEo00'o — 


19 064'0 + 3 1L00‘0 + 111000 — 


10 S+9‘o + 3 £ooo‘o + 691000 — 


19 8+L‘o + 3 z£oo‘o + L6r100‘0 — 


10 S18‘o + 3 1000°0 + 167000 — 


19 098‘O + 3 7700‘0 — 007000 — 


; 
0. re 


—9 


IQ 2-01: 100 + 39500 — #_g0‘0 — 


19 z_OI 
19 :-01 
10 = 01 
19 2-01 
19 z-01 


19 z_0I 


DR" een 


+ 61‘0 + 3 1lr‘o = eo LOL 
+ 9S‘o + 3 06+'‘o Tr VAR 
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fg — tg 
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19 Lyıo‘o + 2 6$$‘o —4 or. 
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19 + 3 Torro‘o + £yzo‘o 
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1Q + 3 #0900°0 + ££zo‘o 


rg 
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rn ee Er 
ee À a 
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aaa nmr eee | 
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19 #01 
19 9-01 - VI‘O — 3 Igz10‘0 — g£10'0 
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Br 
m“ 
H 
Ae ae ae 


m? 


H » 
HR öi; B,;—B,r — 
BE inet, 
re: Ca Ge 


Die Wellen h, u, v setzen sich aus zwei 
hinsichtlich der Schwingungsdauer wesentlich 
unterschiedlichen Anteilen zusammen. Die 
Terme mit den Amplituden Ar, Br, Cr mit 

27 


d 2 T; 
einer Schwingungsdauer T; = — = — stellen 
ver’. u 


in reeller Schreibweise harmonische Funktionen 
dar, die ohne Formänderung mit der konstan- 
ten Verlagerungsgeschwindigkeit cr sich in 
positiver x-Richtung fortbewegen. Die Terme 
mit den Amplituden A» + 43; B:+ B3; Ci+C; 
ergeben in dem Koordinatensystem 


X — (u + 2) = Kconst 


stehende Wellen mit einer Schwingungsdauer 
DST IE 
Teure 
Tab. 2 gibt in Abhängigkeit von der 
Wellenlänge die Schwingungszeiten T: und 
T2, 3 wieder: 


Tab. 2. Schwingungsdauer JT; und T;,; in Ab- 
hängigkeit von der Wellenlänge L für U = 15 
[m - sec], 1 = 1074 [sec], gH =8 - 104 [m?- sec ?] 


5 bs Cr = U—d M Tr Tr, 3 
[106 m[ | [m - sec=:] | [m - sec“) | [Std] [Std] 
30 3,89 555 2 140 15 
25 5,03 488 1 380 14 
20 6,61 426 840 13 
15 8,76 370 475 II 
10 11,39 325 244 9 
5 13,90 294 100 5 
3 14,95 283 19 H 
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Die Terme 


ho = hoAxeï «a 


y ; 
UQ = 8” hoBı ey Kan 


I 


Vert i 
Dos Fe ihoCre? &— at) 


charakterisieren also die eingangs als ”quasista- 
tisch verlaufend” bezeichneten Anderungspro- 
zesse, während die Terme 


3 
hr = ho Aye” &— et) = 


k=2 


= ho [(A2 + A;) cos vMt u 


A;)i sin vMt] e” [»-(u+ :) | 


an 


3 
v : 
UL = = ho > Be? RN = 


k=2 


= a ho [(B: + B3) cos »Mt + 


+ (B,— B;)i sin vMi] fe (u+ )] 


3 
12 KS 0 ; 
vL = rn ho > Ce” (epi!) — 


k=2 


— e ho [(C2 + C3) cos vMt + 


+ (C:— C;)i sin »Mi] e” [=—(u+ 4 


die wegen ihrer relativ geringen Schwingungs- 
dauer als ”meteorologischer Lärm” bezeich- 
neten Änderungsprozesse wiedergeben. 

In Fig. 1a, b, c sind die Verläufe von h, u, v 
sowie hg ua, vo für L = 5 - 106 [m] und e = 
= |e| = 1071, 6=|6| = 10° als Funktion von 
t an der Stelle x = o dargestellt. 

‘Es ergibt sich die für die zeitliche Entwick- 
lung atmosphärischer Felder allgemein cha- 
rakteristische Tatsache, daß der Verlauf der 
kurzdauernden Schwankungen um einen zeit- 
lich zügigen Gang in den meteorologischen 
Elementen wesentlich bestimmt wird durch 
die — im Modell durch & und 6 definierten — 
ageostrophischen Windfelder, die man in der 
Ausgangssituation vorfindet, während die — 
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durch ho, to, va wiedergegebenen — zü- 
gigen Verläufe weitgchend durch die Druck- 
verteilung zur Ausgangszeit festgelegt werden. 

So ergibt sich mit den Näherungen (19a) bei 
einer Auslenkung des Anfangswindfeldes v° 
gegenüber dem geostrophischen Windfelde v,° 
um € = [ele”# für ha eine reelle Amplitude 


von 
/ 
Veen 
m m? 


und eine Phasenverschiebung »x gegenüber 
ho mit e— o, die sich bestimmt aus 


| Ar] = 


sinva =| Aı |”: - = le | sin yp. , 


worin y. den Winkel der im allgemeinen kom- 
plexen Zahl & bedeutet. Bei einer Beschrän- 
kung von |e] < 10° ergibt sich, daß die 
Partialwelle in ho eines durch & = |e | ee 
charakterisierten  Anfangsfeldes gegenüber 
einem durch ¢= o definierten Anfangsfelde in 
v° eine Amplituden- und Phasenänderung er- 
fährt, deren Betrag beschränkt bleibt durch 


[AA:| =e] £5, ee 


m? 


Die gleichen Abschätzungen ergeben sich auch 
für die Amplituden Br, Cr der Funktionen tio, 
VQ. 

Die Tabelle 3 gibt diese maximalen Beträge 
für |e] = 1071 in Abhängigkeit von der Wellen- 
länge wieder. 


Tab. 3. Maximale Amplitudenänderung und 
Phasenverschiebung, hervorgerufen durch eine 
Auslenkung des Anfangswindfeldes v° um = 10-1 


ee gegenüber dem geostrophischen Windfeld. 


Maximale Phasen- 
verschiebungıo 


Maximale relative 
Amplitudenänderung 


]10° ml | A Ar| | vo | 
30 20,5% 15e 
25 3,4 % 2,00 
20 4,4% 2,50 
15 58 % ener 
10 7,6 % 4,4° 

5 9,3 % 5,4° 
I 10,0 % Sn 


Der EinfluB der ageostrophischen Kompo- 


nente u bzw. à im Anfangswindfeld auf Am- 


K. HINKELMANN 


plitude und Phase der Funktionen he, 4a, va 
ist unerheblich. 

Dagegen wächst das Amplitudenverhältnis 
zwischen kurz- und langperiodischen Schwin- 
gungen der Druckfunktion h und der ageostro- 
phischen Komponente u — für kurze Wellen 
näherungsweise ausdrückbar durch 

VE 
TE bzw. [Bs 


bei festgehaltenem 6 — der Einfluß von & ist 
bei kurzen Wellen unerheblich — mit abneh- 
mender Wellenlänge stark an, so daß bereits 
bei Wellenlängen in der Größenordnung von 
einigen 1000 km, wie Fig. 1a, b veranschau- 
lichen, nicht mehr lediglich von einer durch 
den meteorologischen Lärm bewirkten Un- 
ruhe im zeitlichen Gang der meteorologischen 
Elemente gesprochen werden kann. 


Es muß jedoch beachtet werden, daß sowohl 
6 als auch e nicht unabhängig von der Wellen- 
länge angenommen werden können, da sicher- 
lich allgemein 


lim 0.= 0 und um a = 0. 
Lo Lt 


Das ergibt sich schon daraus, daß die Lösungen 
(11, 17, 18) mit 6 > 0, e > 1 in die Lösungen 
der Differentialgleichungen fiir reine Gravita- 
tionswellen, bzw. bei Annahme adiabatischer 
Zustandsänderungen in die Lösungen der 
Diff. Gl. für Schallwellen konvergieren, die 
ja für schr kleine Wellenlängen die tatsäch- 
lichen physikalischen Vorgänge richtig be- 
schreiben. 

Weiterhin dürfte durch Superposition meh- 
rerer Wellenfunktionen verschiedener Wellen- 
länge und Phase eine Glättung der kurzperiodi- 
schen Lärmschwankungen bewirkt werden. 


Da den kurzdauernden Schwankungen we- 
gen der Ungenauigkeit in den Meßdaten für 1, 
v, z.Zt. t=0 keine meteorologische Bedeu- 
tung beigemessen werden kann, wird man 
sich bei der rechnerischen Prognose mit der 
Ermittlung des zeitlichen zügigen Grund- 
verlaufs — im Modell charakterisiert durch 
die Funktionen he, ua, va — begnügen 
müssen. 

Dieser zügige Grundverlauf soll abkürzend 
als Trend” bezeichnet werden, analog der 
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Ausdrucksweise in der statistischen Konjunk- 
turforschung. 

Ignoriert man also diese kurzperiodischen 
Vorgänge hr, ur, vr und sieht man von 
dem geringen Einfluß ab, den die Anfangsfel- 
der der ageostrophischen horizontalen Wind- 
komponenten auf den Verlauf der Funktionen 
he, ua, vo ausüben, so bestimmt sich der 
Trend der Variabeln h, u, v in dem betrach- 
teten Modell lediglich durch die Vorgabe der 
Druckfunktion h° zu einem bestimmten Zeit- 
punkt, und zwar erhält man mit diesen Ver- 
nachlässigungen für 


h = ho er & am) 


Ze as dho eitx—at) 


— £? iho el” Rat), 


ae (20) 


SI 


wenn unter den gestrichenen Größen die oben 
definierten ausgeglichenen Variabeln verstan- 
den werden. Für ein beliebiges Ausgangs- 
feld h° ergibt sich durch Summation über »: 


CO 
er 
n=—c 
oo 
= , Ly* 
Tes os pd bon CO Ch, pm 
P nr 2 T 
oo 
p= 7 Diskrete (202) 
n=—oco 


In den Lésungen der linearisierten Gleichun- 
gen ergibt sich diese Darstellung einfach durch 
Wegstreichen der den meteorologischen Larm 
zum Ausdruck bringenden Funktionen hr, 
ur, V_ sowie der die ageostrophischen Kom- 
ponenten (e, 0) enthaltenden Terme in den 
Amplituden Ax, Br, er 


Da die Funktionen (20a) für alle Anfangs- . 


lagen, bei denen die Felder der horizontalen 
Windgeschwindigkeiten nur unwesentlich aus 
dem geostrophischen Windfeld ausgelenkt 
sind, die Lösungen der Differentialgleichungen 
(1—3) mit für synoptische Belange ausreichen- 
der Genauigkeit beschreiben dürften, so ist es 
naheliegend, die Differentialgleichungen von 
vornherein so zu modifizieren, daß deren 
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Fig. ra. Lärm- und Trendkurven der Funktion h 


h/ ho 


----- h] ho 
ern hQlho 1 


L = $:106m, e=6= 1071 


U = 15 msec 1, gH = 8- 104 m? sec”? 
= i1oO-4 sec! 


Fig. rb. Lärm- und Trendkurven der Funktion u 


Bl vos NL xs 109 ate 0 = 10 
----- u] |vog| U = 15 msec, gH = 8 - 104 m? sec”? 
FY uQl|vgo| | = 1074 sec! 


Lärm- und Trendkurven der Funktion v 


Tugs ce 


v| |vog| L = 5+108 m e = 6 = 1071 
----- 0/ |vog| U= 15 msec, gH = 8-104 m? sec”? 
le vQ/|vog| 1 = 104sec71 


2h ip, 513 ) 


Je, [std] 


oh 
Fig. ıd. Lärm- und Trendkurven der Funktion h = = 


ll len LR RE or — lo“ 
----- hf ho U=15 msec, gH = 8-104 m? sec”? 
& (ho)ı/hg 1 = 104 sec 
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Lösungen nur die Trendkurven liefern, wie sie 
im benutzten Modell durch die Näherungen 
(20a) beschrieben werden, ein Rechenver- 
fahren, das auf die direkte Messung oder 
Berechnung von Windfeldern der Ausgangs- 
wetterlage verzichtet. 

Die zu den Lösungen (20a) zugeordneten 
Differentialgleichungen anstelle der Gln. 
(13) ergeben sich, wenn cr = U — d durch 


U (: — a) angenähert wird, zu 
v2m? 


u Seo 2 Dh, (21) 

v = : hy (22) 

Dh = uy v — Hu (23) 
g x 


‘bzw. in auf die Variable h reduzierter Form 
Di — 2 h, —0 (24) 
XX gH 


Die Mitnahme der ageostrophischen Wind- 
komponente # in dem Vektor 0 ist in diesen 
Gleichungen nur erforderlich, wenn die zeit- 
liche Entwicklung der meteorologischen Va- 
riabeln berechnet werden soll. Wie die Tab. ı 
ausweist, macht die ageostrophische Kompo- 
nente u nur wenige Prozent der geostrophi- 
schen Komponente v aus, so daBes bei Betrach- 
tung des Windfeldes selbst ausreicht, aus der 
vorliegenden Lösung für h bzw. p den geo- 
strophischen Wind durch Gradientbildung zu 
errechnen und den tatsächlichen Wind durch 
diesen zu approximieren. 


V. Numerische Integration 


Wegen der erheblichen mathematischen 
Schwierigkeiten, die der Lösung der allge- 
meinen nicht vereinfachten nichtlinearen hy- 
dro- und thermodynamischen Gleichungen 
entgegenstehen, wird man gezwungen sein, 
wenigstens näherungsweise eine Lösung durch 
eine numerische schrittweise Integration her- 
beizuführen. 

Man begnügt sich zunächst mit der Errech- 
nung der Tendenz der meteorologischen 
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Variabeln a : 
ot 
logische Variable bezeichnet, ersetzt dann nähe- 
rungsweise den zeitlichen Differentialquotien- 
ten durch den entsprechenden Differenzen- 


wenn y eine beliebige meteoro- 


quotienten = und ermittelt durch Annahme 


eines festen Wertes für At die während der 
Zeit At eintretende Änderung Ay, die man zu 
dem zur Zeit t = o vorgefundenen Wert von 
y addiert. 

Durch fortgesetzte Anwendung dieses Re- 
chenverfahrens erhält man eine diskrete Folge 
von Feldern y (t= 0; At; 2At; 3At..), die bei 
hinreichend klein gewählten Intervallen At die 
Lösung der Differentialgleichung hinreichend 
genau zu approximieren vermögen. 

Es bieten sich dabei zwei hinsichtlich des 
Grades der Rechenschwierigkeit und des Um- 
fanges der auszuführenden Integrationsschritte 
unterschiedliche Rechenverfahren an, je nach- 
dem ob die zeitliche Anderung der Momen- 


tanwerte 2 oder die zeitliche Änderung des 


at 
Ja 
Trends = berechnet wird. 
[4 


Es besteht also die Möglichkeit, die Diff. Gin. 
in den Momentanwerten oder in den zeitlich 
ausgeglichenen Werten anzuschreiben. 

1. Für das gewählte Modell einer homo- 
genen Atmosphäre führt die erste Berech- 
nungsmethode auf das folgende Rechenschema, 
das sich unmittelbar durch Auflösung der Gin. 
(1—3) nach den lokalen zeitlichen Differential- 
quotienten ergibt: 


„= -—Uu,— gh, + lv w& vo (25a) 

v = — Ur, — hu » =. (25b) 

3 > IU 

h, = — Uh, + — v — Hu, & Ah (25c) 
£ At 


Diese Gleichungen erfassen sowohl den aus- 
geglichenen zügigen Grundverlauf (Trend) als 
auch die kurzperiodischen Oszillationen (Lärm). 
Bei der Wahl des Zeitinkrementes At ist des- 
halb zu berücksichtigen, daß die ermittelten 
Differentialquotienten nur während einer rela- 
tiv kurzen Zeitspanne als konstant betrachtet 
werden können. 
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end As als Funktion der 


ot ot 

Zeit dargestellt. Beide Werte unterscheiden 
lah 
ot 
halb auch bei der prognostischen Verwer- 
tung der entsprechenden Differenzenquotienten 
in gleicher Weise unterschiedliche Zeitinter- 
valle verwendet werden müssen. 

Zu einer sinnvollen Festsetzung des In- 
krementes At gelangt man, wenn At < T 
gewählt wird. 

Bei Verwendung der Berechnungsgleich- 
ungen (25) führt dies zu 


In Fig. 1d ist 


sich in ihrer Größe erheblich, , wes- 


si 
| Ot 


fe 
i al Bes Ga 
= PET NUE 
Stellt Lin die kürzeste Wellenlänge dar, auf 
deren Erfassung gerade noch Wert gelegt 
wird, bzw. deren Ermittlung aus dem Beo- 
bachtungsmaterial des Stationsnetzes gerade 
noch möglich ist, so ergibt sich für 
(oa 


At < M: 
Lin wird für großräumige (zirkumpolare) 
Darstellungen in der Grössenordnung 10° [m] 
liegen, so daß — vgl. Tab. 2.— At < 1 [Std] zu 
wählen wäre. Der Vorteil dieses Berechnungs- 
schemas liegt darin, daß zur Ermittlung der 
Tendenzen nur elementare Rechenprozesse 
(vier Spezies, räumliche Differentiationen I. 
Ordnung) erforderlich sind, die ohne sonder- 
liche Schwierigkeiten durch automatische Re- 
chenmaschinen ausgeführt werden können. 

Die zwangsläufige Mitberechnung der kurz- 
dauernden Lärmprozesse — ausgezogene Kur- 
ven der Fig. 1 — stellt dabei lediglich ein 
notwendiges Übel dar, da den ermittelten 
Funktionen hr, ur, vr, die näherungsweise 
auch bei einer numerischen Integration der 
Gin. (25) erfaßt werden, so lange keinerlei 
physikalische Realität zukommt, als auch die 
Abweichung vom geostrophischen Windfeld 
im Ausgangsfelde meBtechnisch nicht erfaßt 
werden kann. 

Um nicht eine übertriebene Genauigkeit 
vorzutäuschen, wird es notwendig sein, nach 
Ermittlung des zeitlichen Verlaufs der meteoro- 
logischen Variabeln diese völlig unsicheren 
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Bestandteile durch Glättung mit den üblichen 
Methoden der Statistik wieder zu eliminieren. 

Bei der praktischen Berechnung wird man 
also willkürliche, sich nur im Rahmen der 
MeBfehler haltende Werte für e, 6 annehmen 
können, ohne daß die Genauigkeit der mitt- 


leren Verläufe — h, 7, v — merklich beein- 
trächtigt wird. Naheliegenderweise wird man 
als Ausgangsfeld für die horizontalen Wind- 
geschwindigkeiten das geostrophische Wind- 
feld zugrunde legen, ohne daß jedoch diese 
berechneten Winddaten für den Verlauf der 
Trendkurven oder -felder kritisch sind. 

Der Nachteil des Rechenverfahrens besteht 
darin, daß nur relativ kurze Zeitintervalle 
(einige Minuten) verwendet werden können. 
Dieser Nachteil wird jedoch weitgehend aus- 
geglichen durch die Einfachheit der erforder- 
lichen Rechenprozesse. 

2. Bei Verwendung der hydro- und thermo- 
dynamischen Diff. Gln. mit ausgeglichenen 
Größen führt das Berechnungsschema bei 
Anwendung auf das betrachtete homogene 
Modell auf die Gln. (21—23) bzw. Gl. (24), 


aus welcher man die Tendenz des Trends in h 


durch Auflösung nach h, mit der formalen 
Lösung 


Ban 


Der Vorteil eines solchen Rechenverfahrens, 
in dem die Elimination der kurzperiodischen 
Schwankungen bereits in der differentiellen 
Form der Gleichungen, nicht erst in den 


. Lösungsfunktionen vorgenommen wird, be- 


steht darin, daß wesentlich größere Zeitinter- 
valle At bei Ersatz des Differentialquotienten 
durch den Differenzenquotienten verwendet 
werden können. 

Mit den gleichen Abschätzungen wie unter 
1) ergibt sich wegen der größeren Schwin- 
gungsdauer Tı für 


er 
APR 19 Std. 
I 


Die Anzahl der zur Üburbrückung des 
gleichen Vorhersagezeitaums notwendigen In- 
tegrationsschritte nach Verfahern 2) reduziert 
sich also gegenüber dem Verfahren I) um 
rund das 20-fache. 
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Andrerseits erfährt jedoch das Berechnungs- 
schema insofern eine für den Rechenaufwand 
wesentlich ins Gewicht fallende Komplikation, 
als bereits zur Ermittlung der Tendenz die 
Lösung einer Differentialgleichung — im 
Modell Gl. (24) — erforderlich wird. 

Die sinngemäßen nichtlinearisierten und 
nicht spezialisierten Diff. Gln., die sich nach 
Elimination des Lärmes ergeben, lauten ent- 
sprechend Gln. (21—23): 


LE CR Vi Li Vh, oder bei Ver- | 


l [2 dt 
wendung der Druckfunktion p. | 
LU Da Pate A Tr end à 

Se CRETE Tgp | 
+ dv ‘ V3 und ) 
(26a) 

pce V3:00 

ot Sad 


$,V = 2 dimensionale Vektoren 
3 dimensionale Vektoren 


(26b) 


sl 
= 
| 


zu denen noch die dritte Komponente der 
Bewegungsgleichung, die in Form der sta- 
tischen Grundgleichung eingeführt werden 
könnte, sowie eine thermodynamische Diff. Gl. 
tritt. 

Die Elimination der Lärmvorgänge im all- 
gemeinen Fall führt demnach darauf hinaus, 
in dem Diff. Gl. System die Bewegungsglei- 
chungen für die horizontalen Komponenten zu 
ersetzen durch die durch H. Pmiprs (1939) 
aufgestellte sogenannte quasigeostrophische 
Windgleichung (26a), die nach Linearisierung 
unmittelbar auf die Form der Gln. (25a, 25b) 
führt. 

Es kann zur Zeit nicht abgeschätzt werden, 
welches Verfahren für eine praktische rechne- 
rische Wettervorhersage rationeller sein wird, 
da dies wesentlich von der technischen Leis- 
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tungsfähigkeit der zur Verwendung gelangen- 
den Rechenmaschinen abhängen wird. Jedoch 
dürfte auch das im allgemeinen als wenig 
aussichtsreich betrachtete Verfahren der nu- 
merischen Integration des ursprünglichen 
Gleichungssystems unter Mitnahme meteoro- 
logischer Lärmvorgänge bei Einsatz sehr 
schnell rechnender Automaten erfolgverspre- 
chend sein. 

Bei einer Behandlung der Gleichungen im 
zweidimensionalen Rahmen wird die Auflö- 
sung einer immerhin sehr komplizierten Diff. 
Gl. erspart, im dreidimensionalen Falle — bei 
Ersatz der vertikalen Komponente der Bewe- 
gungsgleichung durch die statische Grundglei- 
chung, der eine Elimination der in vertika- 
ler Richtung verlaufenden kurzperiodischen 
Schallvorgänge bewirkt, — wird man ledig- 
lich auf die Notwendigkeit einer Integration 
über die vertikale Koordinate, also auf die 
Lösung einer gewöhnlichen Diff. Gl. geführt, 
in welcher die horizontalen Koordinaten als 
Parameter auftreten. Ein für eine numerische 
Integration geeignetes Gleichungssystem unter 
Einschluss meteorologischer Lärmvorgänge in 
der Horizontalen wurde von A. ELIAssEN (1949) 
angegeben. 

Die Unkenntnis in den Ausgangsfeldern der 
Windgeschwindigkeit ist dabei, zumindest für 
höhere Breiten, kein ernsthaftes Hindernis, 
da diese Felder ohne Bedenken durch das 
durch den Kräftegleichgewichtszustand defi- 
nierte geostrophische Windfeld ersetzt werden 
können. 
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Le bilan de la rotationnelle 


absolue dans l'atmosphère 


Par J. VAN MIEGHEM, Université de Bruxelles 


(Manuscrit reçu le 27 Juillet 1951) 


Abstract 


The equation of balance of the absolute vorticity is derived from V. Bjerknes Circulation 
Theorem, formulated in an arbitrary moving frame of general coordinates. The flux and the 
rate of production of absolute vorticity are defined. The case of the component of absolute 
vorticity normal to the equiscalar surfaces of a meteorological quantity is considered. From the 
expression of the horizontal convective flux of the vertical component of the relative vorticity, 
a relation between the horizontal meridional flows of relative vorticity and angular momentum 


is deduced. 


Rapportons l’atmosphère à un système de 
coordonnées curvilignes quelconques (x!, x?, 
x3), animé d’un mouvement variable d’un 
point à un autre et d’un instant à un autre. 
Désignons respectivement par Vj, v; et Wi, 
(i=1, 2, 3), les composantes covariantes en 
variables x’ de la vitesse V de l’air par rapport 
au système de référence inertial (système géo- 
centrique fixe), de la vitesse v de l'air par 
rapport au référentiel mobile (x", x?, x3), et 
de la vitesse W du référentiel mobile par 
rapport au référentiel fixe, (V; = W; + v)). 
Soient yi; = y; (i, j = 1, 2, 3) les coefficients de 
la métrique en variables x’ et yÿ = yf les 
mineurs algébriques normés des éléments y; 
du déterminant y = ||y;; ||. Les yj, yŸ et y 
sont des fonctions des x’ et aussi du temps # 
lorsque le référentiel (x', x2, x3) se déforme 
tout en se mouvant; par contre, lorsque ce 
référentiel se meut à la manière d’un corps 
solide, les yi, yy et y sont indépendants du 


temps f. Rappelons que = eo = div W (VAN 
Vi 


MIEGHEM et VANDENPLAS, 1950). 
Ensuite, introduisons les composantes con- 
travariantes & de la rotationnelle absolue 


rot V, on sait que 


6—104318 


OV eV 
Ox! i 


Vath eeu yi = Vs (1) 
où (i, j, k) représente une permutation tour- 
nante de (1, 2, 3) (loc. cit.). Rappelons que 
ER: Vy est la projection orthogonale du 
vecteur rot V sur la normale à la surface de 
coordonnées x3 = Cste au point considéré. 
Notons encore l'identité évidente 


dV y & 
Oxk 


en x! et f. 


Ole 102,3) (1°) 


Enfin, soient p la pression atmosphérique, 
a le volume spécifique de l'air et F; les compo- 
santes covariantes de la force de frottement F. 
Posons 


op do ER 

Dai Oxi D Dy = NY Stet 
IF; IF; NEE" 2 
aa i ae (2) 


où, comme dans (1), (i, j, k) est une permuta- 
tion tournante quelconque déduite de (1, 2, 3). 
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Alors, le théorème de la circulation de V. 
Bjerknes s'écrit, (loc. cit.) 


d En 00 
soy ae ie (3) 


où En, sy et fy sont les projections orthogonales 
des vecteurs rotV, s et rotF sur la demi- 
normale N à l'élément do de la surface o prise 
à l'instant t, (dt = o). 

En remarquant que 


Endo = EN, ôo = yy Éôxiôxi = Vij Ox! Sx, 
(4) 


où (i, j, k) représente une permutation tour- 
nante de (1, 2, 3) et où N, sont les composantes 
covariantes du vecteur-unité N porté par la 
normale N, on établit aisément ae cit.) que 


d(Endo) _ /(dV; mw! 
dt =(5, em 


\ 


it Ox! dxf, (yf, = 1, 2,3) (5) 
Oxi a es 2 


En transformant (sn + fn) do comme Endo a 
été transformé en (4) et en substituant (5) 
dans (3), on trouve la forme rotationnelle des 
équations de la dynamique atmosphérique en 
variables x’, à savoir: 


I dVy & dvi 
Vy dt er 


Ove 3 
Ek a Ei aa sk + fk; 


(i, k= 1, 2, 3) (6) 


d’où, finalement, grace A (1’/), l’équation du 
bilan de la composante & de la rotation- 
nelle absolue rotV, 


1 avy & Lipa ER kei 
ET re Vy (ei — veg] = 
ao Ef, (i,k =r, 23) (6a) 


ou encore, 


a avy & 


vy at + div (tv — vk rot V) = sk + fk, 


(6b) 


De ce bilan, il résulte que le flux, dans le 
référentiel mobile, de la composante & de la 
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rotationnelle absolue rotV est défini par les 
composantes contravariantes 


Cki = EF yi — yk Ei (7) 


et que le taux de production de &, dans ce 
référentiel, s’exprime à l’aide de sk + fr, c’est- 
à-dire à l’aide des solénoïdes isobares-isostères 
et des tubes de tourbillon de la force de frotte- 
ment F. 

Remarquons que C** = 0 et que, par suite, 
le flux de &* à travers la surface x* = CS 
est nul. De plus, le flux CH, (i + k), se décom- 
pose en un flux convectif E* vi le long du courant 
atmosphérique dans la surface x* = CS et en 
un flux non convectif — v*& qui dépend du 
courant atmosphérique v* à travers la surface 
xt = CH Lorsque la surface x = Ces 


dak 
surface substantielle (5 = vk =0} le flux non 


convectif est nul. 


2. Supposons qu’à chaque instant t, les sur- 
faces équiscalaires d’une grandeur météorolo- 
gique, — la température potentielle © par 
exemple, — constituent une famille de sur- 
faces à un paramètre. Identifions ces surfaces 
avec les surfaces de coordonnées x3 = Cste et, 
leurs trajectoires orthogonales avec les lignes 
de coordonnée x3 = var. On a, alors, y31 = 


Y32= YP SYP S00, YByv3= 1, y = V33 Ye 


avec Vy = pri Y22— (yr2)?, et 
en 8 — aS &3 = for V os Pes 
we Vy33 ER VY33 = TOtz V, vz = VY 33 v3, 


(8) 
où € et vz ne sont autres que les projections 
orthogonales de rotV et de v sur la normale z 
à la surface x3 = Cste passant par le point 
considéré. Cette normale est tangente à la 
ligne de coordonnée x3— var. En particulier, 


lorsque x3 = 0, ona v3 = et Vy35 =|O| 


où |\7O| est la valeur algébrique de l’ascen- 
dant de la température potentielle @. 


En substituant (8) dans (6a), on obtient le 
bilan de £, 


SONT zen 2 [Vyx (Ev — vz €] 
Vyx gt Vx axe 
= og oe RE (0) 


bin ach Ad 


nés élu ia: el a lal 


ver 


ee a OR 


7 À bain dt du = = TN aimstttés 


LE BILAN DE LA ROTATIONNELLE ABSOLUE DANS L’ATMOSPHERE 


ou encore, 


1 2Vyxt 1 _2(Vyxév) 
ie ot Mr a xt gi 
\Yx VY» x 

_ vz a) 05 — 
=e Æ 7 ANS = rk 
VYx Oz Ox EVE 


(9a) 


où dz = Vy,, 0x3 et où s; et f. sont les nombres 
de solénoides isobares-isostères et de tubes- 
unités de tourbillon de la force de frottement 
9 

par unité d’aire, dans les surfaces x3 = Cste. La 
somme s; + f. représente le taux de produc- 
tion de € dans x3 = Cste. 

Lorsque, dans l'équation (9a), on néglige les 


ENYA! ae on 


E 
termes en —— „— ) 
| > 33? 


VYx 
obtient une équation établie par S. PETTERSSEN 
(1950). Cette remarque fixe le degré de va- 
lidité de l'équation de Petterssen. Notons que 


I OVyy : 
EA LE. exprime la déformation subie par 
/ dt 
Vx 
1.0 
les surfaces x3 = Cste et que — ave re- 
VP ee 


présente la courbure moyenne de ces surfaces. 

3. Proposons-nous d’étudier le flux de la 
composante € de la rotationnelle dans le cas de 
coordonnées curvilignes triorthogonales (xt, 
x, x3). Dans ce cas, on a: yy = py! =O pour 
is jet yy; = I respectivement pour i = I, 
2, 3. Désignons par x, y, z les axes tangents aux 
lignes de coordonnée x! et par vy, vy, vz les 
composantes coordonnées de la vitesse relative 
v de lair dans le triédre trirectangle xyz, 
supposé dextrogyre. Les axes x, y, z sont 
normaux respectivement aux surfaces x1=— 


Cste, x? = Cste, x3 = Cste. On a donc vy = 
Vy.v! = fy: Win De Vy,, v? = 12: VV 42 
Vz = V7, 03 = Vs : VY35 Ox = Vy,, 0x1, 
dy NN Y2 SEA = VY33 9X3, 


et 


ARE 2 à 
oa re : 
Voy, Nam ox? 
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où wz désigne la composante coordonnée 
suivant l’axe z de la rotation instantanée du 
trièdre xyz par rapport au référentiel inertial. 
Nous avons vu que le flux de ¢ se produit dans 
les surfaces x3 = Cste et que la partie convec- 
tive de ce flux s’exprime, dans le triedre xyz, 
à l’aide des composantes coordonnées (Cv, 
&v,, 0). 

Si nous portons uniquement notre attention 
sur la partie &, de & liée au mouvement de 
air par rapport au référentiel mobile; on 
voit que 


Vase 


avecn( = er We. 
Quant au flux convectif de £, dans le référen- 
tiel mobile, il est défini par les composantes 


coordonnées (Cy vx, Sy ¥), 0) dans xyz. Re- 


marquons que 


ees. : 
| Cy Vx = “rar (ya Vx Vy) — Vy div, v + 
LR M 
2 yn ay (Yale — (x)? ]} 
— I d 
Lx Due Yır dy pe ee Vy) Hr 1% div, Vv == 
a 1144 op at 
| ie 7. ae. { Y22 [ (vy) (vx) ek 
ou 
à Id ı 
div, v a We I (Vy22 vx) ate 
I I ı 7 
ee V ne 
Var dy ( oe } 


On voit ainsi que, dans les surfaces x3 = Cste, 
le flux convectif de la composante normale 
Cy = rot, v de la rotationnelle relative rotv, à 
travers les lignes x!= var et x?=var, dépend 
en premier lieu de la corrélation, le long de 
ces lignes, entre les composantes v, et v, de la 
vitesse relative de l’air et de la corrélation, le 
long des mêmes lignes, entre chacune de ces 
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composantes et la divergence du courant dans 
les surfaces x3 = Cste. 

En particulier, lorsque le système de réfé- 
rence (x, x2, x3) est fixe par rapport à la terre 
et lorsque les coordonnées x1, x2, x3 coin- 
cident avec les coordonnées sphériques A, , r, 
On a: Pir = 1? COS? Y, Y22 = 1, Y33 = I. Le flux 
méridien ®, vers le Nord, à travers le parallèle 
c, de la composante verticale &4 = rot; v de 
la rotationnelle relative rotv par rapport à la 
terre, est alors donné par 


D = p Ex v, 0x = 


I 
1 cos? p dy ¢ 


+ v, div, Vv Ox, 


FCO OD Vy Vy ONE 
(10) 


où nous avons posé 6x == r cos pod et fi = 
ite Wie 

r oy. 

La formule (10) montre que le flux ® dépend 
1°) de la corrélation, le long du parallèle c 
considéré, entre la composante zonale vy du 
vent et la divergence horizontale du courant, 
et 2°) de la variation en latitude de la corréla- 
tion, le long du même parallèle, entre les 
composantes zonale (v,) et meridienne (vy) 
du vent v. En particulier, lorsque les creux ct 
les dorsales bariques sont orientés du SW 
vers le NE, ces deux corrélations sont positives; 
de plus, si la déviation vers l'E, par rapport 
aux méridiens, des creux et des dorsales, 
diminue lorsque la latitude augmente, la 
corrélation entre v, et vy, le long des parallèles, 
decroit vers le Nord. Dans ce cas, les termes 
du second membre de (ro) sont positifs et il 
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en est de meme du flux ® qui est donc dirigé 
vers le Nord. 

Enfin, notons que le flux méridien vers le 
Nord, à travers le parallèle c, du moment 
cinétique relatif de l'air par rapport à l'axe de 
la terre est donné par la formule de Jeffreys 


WY = G1 cos? 0 Vx Vy OX, 
fr 


x Peas : [4 > > 5 is 
où o désigne la masse spécifique de l'air; par 
conséquent, si l’on fait abstraction de la variation 
latitudinale de o, on a 


2 ae + À vx div, v Ox, 


os or cos?y doy 2 


(10a) 


où © désigne une valeur moyenne de o au 
niveau r considéré. La relation (10a) fait 
apparaître le lien entre le flux méridien ® de 
la composante verticale de la rotationnelle re- 
lative rotv de l’air et le flux méridien Ÿ du 
moment cinétique relatif de l'air par rapport 
4 l’axe de la terre. 

Le taux de production, par unité de surface 
horizontale, de la composante verticale rot; v 
de la rotationnelle relative rotv, est égal à 
$:+ fx. Remarquons que 5, est positif ou 
négatif suivant que la rotation qui amène 
l’ascendant horizontal Y,& de volume spéci- 
fique « de lair sur le gradient horizontal — 
Wap de la pression atmosphérique p est cyclo- 
nique ou anticyclonique. Comme # = rot, F, 
on en déduit que f. est positif ou négatif 
suivant que la composante verticale rot: F de 
la rotationnelle de F est cyclonique ou anti- 
cyclonique, c'est-à-dire, d’après une remarque 
de S. PETTERSSEN (1950), suivant que la cir- 
culation horizontale de l’air au point considéré 
est anticyclonique ou cyclonique. 
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CORRESPONDENCE 


Dean Sik: 

The August 1951 issue of Tellus contains a very 
interesting article by Dr. F. A. Berson (1951) on 
some aspects of vortex motion on the earth. In the 
course of presenting a review of several existing 
theories on this general subject Dr. Berson refers 
to a result obtained by Rossby for the force acting 
on a solitary vortex in a resting environment (ba- 
rotropic and non-divergent) on the earth. In parti- 
cular, equation (1) of Dr. Berson’s paper gives the 


I : 
value - R?wß as the force per unit mass acting on 
4 


a simple vortex. This value is obviously obtained 
from Rossby’s equation 


R 
Rp (ro wr, 


which occurs as equation (3) in Rossby’s first paper 
on this subject (1948). The present writer would 
like to point out that this is merely the net Coriolis 
force acting on the vortex and that in the non-di- 
vergent model there will exist a net pressure force 
on the vortex even initially when the vortex has no 
translation and the environment it at rest. That this 
it so is indicated by Rossby’s later result (1949) as 
well as by an analysis carried out by Long (1951). 
(Both of the latter indicate that the initial total 
force on the vortex is approximately one-half that 
given by Rossby’s first result.) ; 


Dr. Berson does not apply this formula in his 
paper, but it may be well to point out this distinc- 
tion at this time. 


REE EE RoE ING ES 


BERSON, F. A., 1951: A synoptic-aerological study of in- 
teraction between vortices. Tellus, 3 172—190. 

Lone, R. R., 1951: A theoreticai and experimental study 
of the motion and stability of certain atmospheric 
vortices. J. Meteor., 8, 207—221 (see pp. 213—215). 

Rosspy, C. G., 1948: On displacements and intensity 
changes of atmospheric vortices. J. Mar. Res., 7, 175— 
187. 


Rossby, C. G., 1949: On a mechanism for the release of 
potential energy in the atmosphere. J. Meteor., 6, 
163—180 (see pp. 178—197). 


Sincerely, 


Norman A. Phillips 
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August issue uf Tellus 1951, 


F. A. Berson on vortex interaction. 


Due to regrettable mistake a wrong drawing has been substituted in 
Fig. 1. The downwards pointing arrow should point upward and to its 
right the minus sign be replaced by a plus sign. — In the text to Fig. 
8a, p. 187, second line from the bottom read »corresponding broken 
symbols» instead of »primed numbers», 
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